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Abstract 

We present calculations of the structure functions !F 2 {x,Q‘^) and 1Fl{x,Q‘^), concentrating on 
small X. After discussing the standard expansion of the structure functions in powers of a.s{Q‘^) we 
consider a leading-order expansion in ln(l/x) and finally an expansion which is leading order in 
both ln(l/x) and as{Q‘^), and which we argue is the only really correct expansion scheme. Ordering 
the calculation in a renormalization-scheme-consistent manner, there is no factorization scheme 
dependence, as there should not be in calculations of physical quantities. The calculational method 
naturally leads to the “physical anomalous dimensions” of Catani, but imposes stronger constraints 
than just the use of these effective anomalous dimensions. In particular, a relationship between 
the small-x forms of the inputs T 2 {x, Qq) and !Fl{x, Qq) is predicted. Analysis of a wide range of 
data for T 2 {x, Q^) is performed, and a very good global fit obtained, particularly for data at small 
X. The fit allows a prediction for 1Fl{x, Q^) to be produced, which is smaller than those produced 
by the usual NLO-in-as(<5^) fits to 1 F 2 {x,Q‘^) and different in shape. 
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1. Introduction. 


The recent measurements of T 2 {x.,Q‘^) at HERA have provided data on a structure function 
at far lower values of x than any previous experiments, and show that there is a marked rise 
in at very small x down to rather low values of [1][2]. Indeed, the most recent 

measurements have demonstrated that the rise persists for values of as low as 1.5 GeV^. 

The qualitative result of a steep rise at small x initially led to surprise in many quarters. The 
main reason for this was that standard methods used to fit the data were based on the solution of 
the Altarelli-Parisi evolution equation [3] (along with convolution with coefficient functions) at the 
two-loop level, using flat input parton distributions at starting scales of Qo r\j 4GeV^ (e.g. [4]). 
This method followed the reasoning that steep behaviour can only come about from perturbative 
physics (the Donnachie-Landshoff pomeron used to describe soft physics has behaviour [5], 

and we will take steep to mean any powerlike behaviour steeper than this), and that a starting 
scale for perturbative evolution should be high enough for as{Q^) to be fairly small (^ 0.3) and 
to avoid any significant corrections from “higher twist” {Aqq^/Q'^) corrections. This procedure 
results in an effectively steep^ behaviour at small x [6], but only after a long evolution length, and 
therefore at values of » 4GeV^. 

Thus, the data led to a degree of optimism amongst those advocating an alternative description 
of small-x structure functions, i.e. using the BFKL equation [7]. This equation provides the 
unintegrated gluon Green’s function which includes the leading power of ln(l/x) for any power of 
as (where ag is taken to be fixed). It was traditionally solved analytically in the asymptotic limit 
X ^ 0, or numerically for finite x, and predicted a powerlike behaviour of for the gluon 

distribution function, where A = 41n2Q;s and oig = (3/7r)Q!s, i.e. A ~ 0.5 if Ug ~ 0.2. This was 
assumed to lead to T 2 {x,Q‘^) behaving like x~^ and could be claimed to be in rough qualitative 
agreement with the data, even if A was somewhat high. It could also be seen as some justification 
for choosing powerlike inputs (with A ~ 0.2 — 0.3) for the parton distributions (e.g. [8]), which 
could then enable a good fit to the data using the Altarelli-Parisi equation. 

However, it was also convincingly demonstrated that it was possible to generate the observed 
steep behaviour by being a little less conservative concerning the region in which perturbative 
evolution could be applied. Gliick, Reya and Vogt had in fact predicted a sharp rise in F 2 {x,Q‘^) 
at small x, even for ~ IGeV^, by using two-loop evolution from roughly valence-like parton 
distributions at a starting scale of Qq = 0.34GeV^ [9]. This attracted criticism not only on the 
grounds that higher twist corrections should be important in the region of evolution, but also 
because ag{Q‘^) was as high as ~ 0.5 at the lower end of the range, and therefore perturbation 

^ We use the term effectively steep since the slope at any x and is a function of these variables rather 
than a constant, i.e. 7-2(x,Q^) ~ x~^ where A ~ 0.7^ ln~2(0.1/x), and is smaller than any 
power of X as X ^ 0. 
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theory itself should be questionable. Starting at a higher scale, Qq = IGeV^, Ball and Forte were 
able to fit the small-x data using their double asymptotic scaling (DAS) formula [10], which is a 
simple, but very accurate, approximation to the solution of the one-loop evolution equation with 
flat^ inputs and which is valid in the region of small x {x ^ 0.01). They also showed that two-loop 
evolution with flat inputs starting at ~ 2GeV^ could fit the data available at that time very 
well [11], and had a similar shape to the DAS result (as indeed it must, since both fit the data). 

This state of affairs clearly left scope for argument about the real underlying physics describing 
the small-x behaviour of Those who used Altarelli-Parisi evolution from small scales 

could be accused hrstly of working in regions where perturbation theory was questionable and, 
perhaps more importantly, of ignoring terms of higher order in (but also higher order in ln(l/x)), 
which seemed from the BFKL equation to have very important effects. Conversely, those who used 
the BFKL approach could be accused of ignoring all but the leading-ln(l/j;) terms (and hence 
ignoring the large-rc data) and also of working in a less well-defined theoretical framework than the 
renormalization group approach based on the factorization of collinear singularities [12]. Starting 
from an input for the parton distributions with A ~ 0.25 at values of Qq ~ 4GeV^ was taking 
the best of both worlds. However, this lacked a real justification for the choice of input, which 
was signihcantly steeper than that expected from non-perturbative physics, but rather smaller 
than that from the BFKL equation, and also ignored potentially important ln(l/x) terms in the 
evolution. 

A significant step forward in the investigation of small-x structure functions was the devel¬ 
opment of the /cT-factorization theorem [13] [14]. This is the prescription for the way in which 
an off-shell photon-gluon scattering amplitude can be convoluted with the unintegrated gluon 
Green’s function calculated using the BFKL equation to provide the small-x structure functions 
themselves (once convoluted with a bare, off-shell gluon density). Hence, it enables one to find 
effective moment-space coefficient functions (or in the case of and massless quarks, 

where dJ- 2 {x,Q‘^)/dlnQ‘^ is calculated, a mixture of coefficient function and anomalous dimen¬ 
sions) within the BFKL framework. Numerical calculations performed using this method were able 
to match the available data in a qualitative manner [15], as did similar calculations [16] using a 
modification of the BFKL equation, i.e. the GGFM equation [17]. 

However, the /cr-factorization formula was also shown to be very important if one insisted on 
working within the rigorous framework of the traditional renormalization group approach [18]. In 
this approach the infrared poles in the calculated coefficient functions, or photon-parton scattering 
amplitudes, are removed order by order in ag and absorbed into the essentially nonperturbative 
parton distribution functions. In order to guarantee the independence of physical quantities on the 

^ Parton inputs behaving like x~^, as opposed to are referred to as flat, i.e. the parton density 

rescaled by x is flat. 


2 



factorization scale the latter then evolve in ln(/x|p) according to the renormalization group equa¬ 
tions governed by calculable perturbative anomalous dimensions. By showing how /^T’-factorization 
fits within the collinear factorization framework, Catani and Hautmann were able to calculate all 
the renormalization group anomalous dimensions to lowest nontrivial order in a* for each power of 
ln(l/x), and similarly for a number of coefficient functions. 

In order to explain this properly we digress for a moment. We may write the Altarelli-Parisi 
splitting functions as 


1 


P{x,as{Q^)) = -Y,a:{Q^)(y^ 


n — 1 


\n^{l/x) 

a-nm -i -1- terms less singular as x ^ 0 , 

m\ 


n=l ^m=0 

or taking the Mellin transform of the splitting function (weighted by x) 

'y{N,asiQ‘^)) = [ x^P{x,as{Q^))dx, 

Jo 

and working in moment space, 

oo n 

7(iV,a,(Q2)) = + 


( 1 . 1 ) 


( 1 . 2 ) 


n=l 


m=l 




(1.3) 


where A(N) is regular as ^ 0. In the normal loop expansion one solves the evolution equations 
order by order in a*. We may however write 


- OO y OO 1 TTi _1 / "I / \ \ 

P{x,as{Q^)) = - a”((5^)( anmcCiQ'^)— -j ——] -h terms less singular 

X \ ml J 

n = 0 ^772 = 1 ' 


as X —> 0, 


(1.4) 


or in moment space, 


OO OO OO 

7(A^, as(g")) = <(Q") dr,ma^{Q‘^)N-^ + ^ a:{Q^)an{N), (1.5) 

n=0 m=l n=l 

where the an{N) are regular as iV ^ 0, but have singularities at negative integer values of N. 
When one calculates the structure function by converting back to x space (i.e. performing the 
inverse Mellin transform), these singularities lead to contributions suppressed by powers of x, and 
so are negligible at small x. Therefore, ignoring these contributions and expanding the an{N) about 
= 0 we get 


7(Ar,a,(Q2)) = ^ a_a-(Q2)^-rn = (1.6) 

n=0 m=l —n n=0 


where the expressions are strictly convergent only for |A^| < 1. Thus, in order to include the leading 
powers in ln(l/x) for the splitting function, we have to take the n = 0 part of the expression (1.6) 
for the anomalous dimension. Next to leading order (NLO) in ln(l/x) is the n = 1 part, and so on. 
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Using the /cr-factorization formula Catani and Hautmann demonstrated that within the renor¬ 
malization group framework 'jgg(N,Q‘^) and were the same renormalization-scheme- 

independent expressions as the effective anomalous dimensions [19] given by the BFKL equation. 

i.e. 



7°^(a.(g2)/lV) = ( 1 - 7 ) 


where 7 °g(as(g^)/A^) is given by the iterative solution of 


1 = 


N 


-xiil) 


( 1 . 8 ) 


and, 

U(7) = 2V’(1) - V’(7) - - 7)- (1-9) 

This solution as a power series in as{Q‘^)/N exists only for |A^| > A(Q;s(g^)). The anomalous 
dimension develops a branch point at = A(Q;s(g^)), but the series in (1.7) is analytic and 
convergent outside the circle |A^| = A(g^). In fact, this anomalous dimension leads to the gluon 
distribution function having an asymptotic powerlike behaviour of so it does not make 

sense to talk about its Mellin transform for it simply does not exist in this region of 

moment space. However, the inverse Mellin transform may be performed by analytic continuation 
into the region 'ReN < A(go) (but |A^| > A(go)) from 'ReN > A(go)) if so desired. 

Catani and Hautmann also derived expressions for 'yjj:{as{Q‘^)/N) and 'yjg{as{Q‘^)/N) in 
certain factorization schemes {'y^f{as{Q‘^)/N) and 'y^g{as{Q‘^)/N) being zero), and also for the 
coefficient functions ^{as{Q^)/N), C(^ ^{as{Q^)/N), C 2 i{asiQ‘^)/N) and C 2 i{asiQ'^)/N) (all 
zeroth-order quantities being zero except which is unity). This facilitated calculations of 

structure functions within the normal renormalization group framework, but including much of 
what is often called the BFKL physics (i.e. the leading-ln(l/rE) terms), and indeed, a number 
of calculations were performed [20]-[23], and in most cases comparisons with data made. These 
calculations used different methods of solution, made rather different assumptions and used different 
ansatze for unknown terms. Consequently different results were obtained. The conclusions which 
could be drawn regarding the inclusion of the leading-ln(l/a;) terms depended very much on which 
of the approaches was taken. However, it seemed that by including these terms it was not possible 
to improve upon the best fits for the small-x data using one- or two-loop evolution from soft inputs 
[20] [22]. Indeed, many ways of including them made the fits significantly worse, and this seemed 
to be universally true if the fits were more global, i.e. constrained by large-x data [20]. Also, 
it seemed that there was a very strong dependence on the factorization scheme used to perform 
the calculations when including the leading-ln(l/x) terms [20] [22] [24] [25], and a number of new 
factorization schemes were invented, e.g. the SDIS scheme [24] and the Qq scheme [25]. 
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The high precision of the most recent HERA data constrains theory far more than previously, 
and has changed the above picture somewhat. The best recent global fits seem to come from those 
intermediate approaches which use NLO perturbation theory with a quite steep (and completely 
unexplained) input for the singlet quark with A ~ 0.2 and a similar form of small-x input for 
the gluon [26] (unless Qq is less than ~ 4GeV^, in which case the gluon must be flatter or even 
valence-like). Fixed order perturbation theory using flat or valence-like inputs and low Qq fails 
at the lowest x values, and for fits to the small-x data alone relatively steep inputs for the singlet 
quark, i.e. A ^ 0.2, seem to be absolutely necessary [27]. Moreover, approaches including the 
leading~ln(l/x) terms now seem to fail [28] [29] in practically all factorization schemes. 

In this paper we will take issue with the above conclusions. In particular we will demonstrate 
that the apparent failure of approaches using the leading-ln(l/x) terms, and certainly the factoriza¬ 
tion scheme dependence, is due to incorrect methods of incorporating these terms. Indeed, Catani 
has already shown how to obtain factorization-scheme-invariant results in the small-x expansion 
by writing evolution equations in terms of the physical quantities, the structure functions and 
“physical anomalous dimensions”, rather than in terms of parton densities and of the usual anoma¬ 
lous dimensions [30]. In this paper we will go further and show that the correct leading-order, 
renormalization-scheme-consistent (RSC) calculation of the structure functions must naturally in¬ 
clude some leading-ln(l/x) terms in the form of these “physical anomalous dimensions” . It also 
provides limited predictive power at small x, giving justification for the input J^ 2 (x,Qq). We will 
discuss this method of calculation, then make detailed comparisons to data, and demonstrate with 
the aid of the new HERA data that it leads to a very good global fit to all J^ 2 {x, Q^) data. Indeed, 
the complete RSC calculation, including leading-ln(l/x) terms, is clearly preferred by the latest 
data, particularly that at small x. 

We note that a very brief presentation of the complete RSC calculation of structure functions 
has already appeared in [31].^ However, this current paper gives a complete, and very detailed 
discussion of the correct calculation of structure functions, as well as examples of commonly en¬ 
countered pitfalls, whereas there is only the barest outline of the full calculation in [31]. Moreover, 
there is a far more comprehensive presentation of the comparison with experimental data, and with 
alternative approaches, in the current paper. 

This paper will be structured as follows; after giving a brief outline of the different possible 
types of scheme dependence in the calculation of structure functions, we also quickly review the 
work of Catani, illustrating that it is indeed very easy to guarantee factorization-scheme-invariant 
results. We then give a very detailed description of the calculation, within moment space, of 
structure functions in various expansion schemes using the normal parton language. The first 

^ As in this shorter presentation, the present paper deals with the heavy quark thresholds in a rather 
naive manner, i.e. the quarks are taken to be massless, with a particular flavour becoming active only 
above a certain Q^. 
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part of this, regarding the normal loop expansion, will be largely a review, but will highlight 
points usually not discussed in detail, particularly the role played by the scale at which the parton 
distributions are input and the form of the inputs. The second part, discussing the leading- 
ln(l/a:) expansion, will present the only correct way to perform this expansion. One finds that by 
remedying the factorization scheme dependence one is forced towards an alternative derivation of 
Catani’s factorization-scheme-independent anomalous dimensions. However, we will also show 
that, in order to make an ordered calculation of quantities, more care is needed than simply 
working with a factorization-scheme-independent set of variables, and demonstrate that there 
is a certain degree of predictive power for the form of the inputs for the structure functions at 
small X. We also explain why the standard solutions using the small-x expansions are strongly 
factorization scheme dependent. To conclude this section we present the argument that there 
is a unique renormalization-scheme-consistent calculation of structure functions, which applies 
to both large and small x. We present this calculation for both the currently academic case of 
the nonsinglet structure functions and for the phenomenologically important case of the singlet 
structure functions. We then discuss how we move from moment space and obtain our x-space 
solutions, and the qualitative form these solutions must take, i.e. our best attempt at predictions. 
After this long theoretical presentation we consider the comparison with experiment. We fit the 
data for J- 2 {x,Q‘^) using the renormalization-scheme-consistent solutions, and compare to global 
fits at NLO using the normal loop expansion. We also comment on alternative fits to the data made 
by other groups, and on determinations of ag- We conclude that the full renormalization-scheme- 
consistent calculation gives the best global fit to structure function data, particularly at small 
X. Finally we investigate the phenomenological consequences for Tl{x,Q‘^), and also preliminary 
indications for the charm structure function. Good measurements of either (but preferably both) 
of these quantities, particularly at small x, would help determine whether the approach developed 
in this paper, and hence the inclusion of leading-ln(l/x) terms in the calculation of structure 
functions, is indeed correct. 

2. Scheme and Scale Choices. 

For simplicity we work in moment-space for much of this paper, i.e. define the moment-space 
structure functions by the Mellin transform, i.e. 

F{N,Q^)= [ x^-^T{x,Q‘^)dx. (2.1) 

Jo 

The moment-space coefficient function is defined similarly but, as with the definition of the anoma¬ 
lous dimension, we define the moment space expression for the parton distribution as the Mellin 
transform of a rescaled parton density i.e 

f{N,Q^)= [ x^i{x,Q^)dx. (2.2) 

Jo 
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Let us consider the most general moment-space expression for a structure function, i.e. the sum of 
the products of the expressions for hard scattering with a certain parton (the coefficient functions) 
with the corresponding, intrinsically nonperturbative parton distributions. 

F{N, Q^) = C°'{N, //|’), (2-3) 

a 

/ii? is a factorization scale separating the ultraviolet physics from the soft infrared physics, and on 
which the left-hand side of (2.3) is independent, is the renormalization scale, i.e. the appropriate 
scale at which to define the coupling constant, and on which F{N, Q^) is again independent. 

Despite being intrinsically nonperturbative, the parton distributions evolve according to the 
perturbative renormalization group equation 

dfaiasijJ^R^^fiR/l^p) ^ ^7„b(a^(//2^),//y/i|)/b(a«(^|j),^y^|,). (2.4) 

In principle there are many choices to be made when performing a perturbative calculation of 
structure functions to a finite order. However, two of the choices left open in the above expressions 
are common to all perturbative calculations in quantum field theory; the choice of renormalization 
scheme and subsequently the choice of renormalization scale. When one removes ultraviolet di¬ 
vergences from perturbative calculations, there is always a freedom in what type of regularization 
is used and/or how much of the finite part of a calculation is removed at the same time as the 
divergent parts. Although the all orders calculation of a physical quantity is independent of the 
convention, the renormalization scheme used, the perturbative expansion of the quantity is not. 
Hence, neither is the definition of the expansion parameter, the running coupling constant. All 
couplings in quantum field theory satisfy a renormalization group equation, 

E = -«“(>■")). (2-5) 

where /x is the scale at which the coupling is defined, and the convention of the minus sign is 
introduced in order to make the /3-function for QCD positive. The solution to this equation 
depends on the coefficients bn and on a scale A. This scale, and the value of the coefficients beyond 
re = 1 are renormalization scheme dependent {bo = (11 —2AIj/3)/47r and bi = (102 —38Alj/3)/167r^ 
in all renormalization schemes which use dimensional regularization, where Nf is the number of 
quark flavours). Hence, a choice of renormalization scheme amounts to a choice of which expansion 
parameter is to be used, and consequently the form of the perturbative expansion. One would 
hope the choice would be such as to make the series converge as quickly as possible, but this is 
difficult to guarantee. Conventionally the MS renormalization scheme is used. Also the choice of 
the appropriate scale to be used in the coupling for a given process must be made. Sometimes 
this is relatively clear, but in deep inelastic scattering it is not entirely obvious. Traditionally, the 
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simple choice is taken, and indeed, the data seem to favour this choice [32], Hence, we 

will assume this to be the case from now on. 

Having settled for a particular choice of renormalization scheme and renormalization scale the 
ambiguities due to ultraviolet regularization have been dealt with. Our fundamental equations 
become 

F(iV,Q2) = (2.6) 

a 

and 

( 2 - 7 ) 

where «« is defined in a particular renormalization scheme. The remaining ambiguities are due to 
the particular problems in calculating quantities in QCD, i.e separating the physical quantity into 
the perturbative coefficient function and the intrinsically nonperturbative parton distribution. To 
begin with we have the freedom of choosing the factorization scale fip- As with renormalization 
scheme dependence this does not affect the all-orders calculation, but does affect the form of the 
perturbative expansion, since it affects the scale at which the coupling in the coefficient function 
is evaluated. One might imagine that it is desirable to choose to be large in order to make 

the expansion parameter asifij?) as small as possible, and indeed is nearly always chosen to be 
equal to the hard scattering scale Q^. We shall also make this simple choice. 

This leaves us with our defining equations^ 

FiN, g2) = ^ C“(iV, as{Q^))UN, Q2), (2.8) 

a 

and 

= (2.9) 

In principle there are now choices to be made for the scale at which to begin the evolution of 
the parton densities, Qq, and for the inputs for the partons at this scale. This question will be 
dealt with in detail later in this paper, so we will leave it for the moment. This still leaves us 
one more freedom in our calculation, i.e. how we choose to remove the infrared divergences from 
the bare coefficient functions and hence how we define our parton distributions. Starting from 
any particular choice for the definition of parton distributions one may always choose a new set of 
parton distributions by an invertible transformation 

fa{N, Q^) = Y, Uab{N, a,)MN, Q^), (2.10) 

b 

^ Because we consider only massless quarks, once we set = njp the only scale dependence in the 
coefficient functions comes from as{Q^). 



where Uab{N, as) has a power series expansion in as such that Uab{N, a^) = Sab + 0{as)- The struc¬ 
ture functions will clearly be unchanged as long as the coefficient functions obey the transformation 
rule 

C-{N,as) = {U^)JiN,as)C\N,as). ( 2 . 11 ) 

By substituting fa{N,Q‘^) = Ylib^abi^i^s)fb{N,Q‘^) into (2.9) we easily find that the new par- 
ton densities evolve according to the standard evolution equations but with the new anomalous 
dimensions 

%b{N, as)=Y.Y. ocs)^od{N, as)U^j^{N, a,) + ^ a,). (2.12) 

The matrix U must obey a number of conditions in order that physical requirements on the parton 
distributions are maintained, e.g. flavour and charge conjugation invariance, fermion number con¬ 
servation and longitudinal momentum conservation (see for example the second of [18]). However, 
none of these needs to be satisfied simply in order to keep the structure functions unchanged. 

The transformation defined above is called a change of factorization scheme and is constructed 
precisely so that the physical structure functions are left invariant. However, it is important to 
realize that, unlike the changes in renormalization scheme, a change in factorization scheme leaves 
the expression for the structure functions unchanged not only to all orders, but order by order in 
as, and calculations performed carefully at a given well-defined order in one scheme will lead to 
precisely the same results for the structure functions as those in another scheme. Also, the coupling 
constant to be used depends only on the ultraviolet renormalization.^ We will illustrate this point 
in §4, where we will calculate structure functions to a well-defined order in a variety of expansion 
schemes and demonstrate factorization scheme independence of our expressions. However, we will 
first illustrate Catani’s recent proposal for the construction of factorization-scheme-independent 
structure functions. 


3. Evolution Equations for Structure Functions. 

It is, as Catani noticed [30], very simple to obtain factorization-scheme-independent expres¬ 
sions for the structure functions, or more precisely, factorization-scheme-independent effective 
anomalous dimensions governing the evolution of the structure functions (reflecting the fact that 
we cannot make a concrete prediction of the structure functions at a given scale using QCD, but 
only how they change with Q^). In order to obtain these effective anomalous dimensions all one 
has to do is eliminate the parton densities from the equations (2.8) and (2.9). 

^ It is perfectly possible (if somewhat perverse) to choose the MS scheme to remove ultraviolet diver¬ 
gences, but the MS scheme to remove infrared divergences from the bare coefficient functions. In this case 
it is the MS scheme coupling constant that must be used. 
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In order to demonstrate this, let us now be a little more careful in our definition of the relevant 
structure functions. There are two independent structure functions F 2 {N,Q‘^) and Fl{N,Q‘^). In 
general we may write 

F,iN, Q2) = Q^) + Ff^^iN, Q^) {i = 2,L), (3.1) 

where the singlet and nonsinglet structure functions are defined by 

Ff(N, Q2) = d{N, Q^) + Cf{N, Q^), (3.2) 

and 

Nf 

Fl^^iN,Q^) = CdiN,a,)^e]d^iN,Q^), (3.3) 

i=i 

where Nf is the number of active quark flavours, f^{N,Q‘^) and fd{N,Q‘^) are the singlet and 
nonsinglet quark distribution functions respectively, and g{N,Q‘^) is the gluon distribution. In 
any factorization scheme obeying the requirements of flavour and charge conjugation invariance 
the renormalization group equations for the singlet and nonsinglet sectors will be decoupled. The 
equations for the nonsinglet distributions will be ordinary differential equations, 

\Iq 2 =lNs{N,as)fdiN,Q^), (3.4) 

while those for the singlet sector are coupled 

d (rff{N,a,) ^dN,as)\(d{N,Q^)\ 

(ilnQ 2 g{N,Q^) J \jgf{N,as) jgg{N,as) ) V g{N,Q^) )' 


Let us first consider the simple case of the nonsinglet structure function F^^{N,Q^). Multi¬ 


plying both sides of (3.4) by Yl7=i 6? can clearly write 


dlnQ^y Cf ^(AT, a,) ) ~ Ci^^(iV, a,) 

which becomes the factorization-scheme-independent equation 


(3.6) 


^ r2,iV5(iV, as)F2^^{N, Q^), (3.7) 

where r 2 ,iV 5 (A^, a^) = jNsiNjas) + dln{C 2 ^{N,as))/dlnQ'^. Therefore, we have an effective 
anomalous dimension governing the evolution of each of the nonsinglet structure functions, and 
clearly r)^‘^(AI, Og) must be a factorization-scheme-independent quantity (and is in principle mea¬ 
surable). The solution to this equation is trivial; 


Ff^(iV,Q2) 


Fi^^iN,Ql)exp 



^2,ns{N, as)d\n 


(3.8) 
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The situation for the singlet structure functions is more complicated. As we see from (3.5) the 
evolution equations for the singlet quark density and gluon density are coupled. However, using 
(3.2) for i = 2, L we may solve for the parton densities in terms of the two structure functions 
and the coefficient functions. Substituting these into (3.5) we then obtain the coupled evolution 
equations 

d fFi{N,Q^)\ _ ft22{N,a,) t2L{N,a,)\ f Fi{N,Q^)\ 

dlnQ^\F^{N,Q^)) [rL2{N,as) tLL{N,as) J \Fi{N,Q^) J ' 


The expressions for the physical anomalous dimensions, f 22 (A^, a*), i' 2 L{N,Cis), etc., are straight¬ 
forward to derive in terms of anomalous dimensions and coefficient functions in any particular 
factorization scheme using the above procedure, but result in rather cumbersome expressions. It is 
simplest first to define a factorization scheme such that Tf (A", Q^) = f^{N, Q^), i.e. (A, ag) = 1, 
= 0 (this is generally known as a DIS type scheme [33]®). In terms of the coefficient 
functions and anomalous dimensions in this type of scheme we have 


f22{N,ag) = -fff{N,ag) - “''| 7/g(A, ag), 


i'2L{d^, Otg 
fL,2(A^, Os 


7/9 (^, 


ag 


Cl{N,agy 

= CUN, ag)jgfiN, ag) - (A, ag)jgg{N, a. 


dC[(N, 


ag 


Cf(N a _ (CUN.a.))\ 

dlnQ 2 Cl{N,ag) ’ 


dlnQ2 


+ C[{N, a^)7//(A, 


ag 


tLL{N,ag) = 7gg(A, a*) + 


dln(C£(A,a,)) , Ci{N,a 


din 


-k 


Ci{N,ag 


-lfg{N, ag 


(3.10) 

Thus, we have a factorization-schemeAnvariant set of anomalous dimensions governing the 
evolution of the structure functions. Before going any further let us remark that we believe there is 
a (purely technical) problem with the above expression. As is well known (and as we will discuss in 
the next section), Fl{N,Q‘^) starts at an order of ag higher than F 2 {N,Q'^). Because of this there 
is an intrinsic asymmetry in the above definitions, with r 2 L(A, a^) beginning at zeroth order in 
otg, t' 22 {N, ag) and T ll{N, Og) beginning at first order, and T l 2 {N, a*) beginning at second order. 
The fact that one of the physical anomalous dimensions has a part at zeroth order in ag seems 
against the spirit of the perturbative approach. In practice, the result is that, if one solves the 
evolution equations order by order in ag the solutions are very different from those obtained from 
an order-by-order solution using the parton densities, which we know work very well for all but 
possibly the lowest values of x yet probed. (We will discuss this point rather more in §4.) 


® We call it a DIS “type” scheme because satisfying the above requirement still leaves freedom in how 
we may define the gluon density, and thus we are still considering a family of schemes. 
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A trivial modification of Catani’s approach is therefore to accept that Fl{N,Q‘^) contains an 
extra power of ag, and to define the new structure function Fl{N,Q‘^) = FL{N,Q‘^)/{as/{2Tr)). 
The longitudinal coefficient functions are likewise changed to C1{N,as) = C'£(A', as)/(«s/(2vr)), 
and the singlet evolution equations become 

d fFi{N,Q^)\ _ fT22{N,ag) T2LiN,ag)\ f Fi{N,Q^)\ 
dlnQ^[F^{N,Q^)) [rL2{N,ag) TLL{N,as) J \Fi{N,Q^) J ^ 


where the r(A^, q;s)’s are defined precisely as in (3.10), but in terms of of C1{N,as) rather than 
C1{N,as). This procedure restores the symmetry between the physical anomalous dimensions, 
and makes the order-by-order-in-ojs calculations essentially the same as when using evolution of 
parton distributions. It is, of course, trivial to obtain the physical Fl{N,Q‘^) from Fl{N,Q‘^)7 
Having made our redefinition of the quantities with which we work, we now have a direct 
relationship between possible calculations using the evolution equations for structure functions 
and the solutions using the parton densities. At present the parton anomalous dimensions and 
coefficient functions are known to order a^. It is easy to see that this allows us to derive each of 
the T’s to order a^. For example, at first order in as they are 


=7;;(iV,a,)- 


0 , 1 ^ at ^ “s) _ 0,i 


Clo,iiN,as) 


as), 




(Xa 


r^’}(N,as) 

"" CU,{N,asy 

Tl[^{N,as) = Cl,yN,as)^^J{N,as) - Cl./N, as)^^^(N, ag) + Cl,yN,as)^^/iiN,as 


{Clo,iiN,ag))^ , 




a. 


0,7 ^r ^ 




as), 


(3.12) 

where the super-subscript n,l denotes the (n + l)-loop quantity. The O(a^) expressions for the 
physical anomalous dimensions are straightforward to derive, but are rather complicated. Similarly, 
from the known expansions of the parton anomalous dimensions and coefficient functions in the form 
a'^Y)m=i-n^fn{as/N)"^, we can calculate rl^{N,ag) and r 2 , 2 (iV,as), T^L{N,as) and 
(where both are trivially zero), and r 2 ^(A", Og) and r 22 (A^, a^). This is the same order as for the 
parton anomalous dimensions, with the longitudinal anomalous dimensions having similar structure 
to the gluon anomalous dimensions and the r 2 a(A^, cks)’s having similar structure to the quark 


^ Similarly, it is also best to work with the rescaled nonsinglet structure function {N,Q^). Doing 
this, the analogous procedure to (3.7) and (3.8) leads to F^’®(A,Qo) beginning at zeroth order in as{Qo), 
as does F^^^{N,Qq). 
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anomalous dimensions. The exact form of these physical anomalous dimensions may be calculated 
using the expressions for the parton anomalous dimensions and coefficient functions in the standard 
DIS scheme in [18]. The resulting expressions are relatively simple, being 


^Idc^s/N) 

^Idc^s/N) 


(2vr; 




n .=0 




ac 


- oo 

-^(i7;,(a./iV) + Y.W>Ja.lN)T), 

^ ' n=0 

7°g(a«/iV), 


(3.13) 


where q and Cf ^ are the one-loop contributions to 'yjg{aa/N), C[^(as/N) and 

Cl Yois/N) respectively. We note that each of these anomalous dimensions is renormalization 
scheme invariant as well as factorization scheme invariant, as we would expect for leading-order 
physical quantities. 

Solving the evolution equations for the structure functions using any subset of the currently 
known physical anomalous dimensions guarantees a result which is factorization scheme indepen¬ 
dent. However, since the evolution equations are coupled, a simple expression for the solutions, such 
as (3.8), is impossible. There is, in fact, considerable freedom in how we may solve the equations. 
We could, for example, simply put all of the anomalous dimensions currently known into (3.11) and 
then find the whole solution. Alternatively, we could solve using just the order anomalous di¬ 
mensions and then try to perturb about this solution in an ordered manner. These two approaches 
would lead to rather different answers, but both would be factorization scheme independent. The 
problem of obtaining a correctly ordered solution for the structure function will be discussed in 
detail in the next section. We will initially use the familiar parton distributions and coefficient 
functions, and show that, even when using this approach, if we solve producing a well-defined 
expansion for the structure functions, we automatically avoid the problem of factorization scheme 
dependence. 


4. Ordered Calculations of Structure Functions. 

There are in principle many different expansion methods one may use when obtaining solutions 
for the structure functions. The standard one is simply solving order by order in as- But there is 
also the expansion in leading powers of ln(l/x) for given powers in (or equivalently in powers 
of N~^ in moment space), as we have already discussed. One can also combine the two expansion 
methods, and indeed, we will later argue that this is the correct thing to do. Nevertheless, we 
will begin by outlining the procedure for making a well-ordered calculation of structure functions 
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using the standard loop expansion. Although this is well known, we feel it is worth presenting it 
pedagogically, and making some points which are not usually highlighted, especially concerning the 
role of the starting scale. We may then discuss the more complicated cases of the expansion in 
leading powers of ln(l/a:) and the combined expansion. 


4.1. Loop Expansion. 


We begin by introducing some new notation. In order to solve the evolution equations for the 
parton densities order by order in Ug and hence obtain expressions for Q^) and Ff{N, Q^) 

we make use of equation (2.5) to rewrite the evolution equation for the nonsinglet parton density 
as 

d f^^(N 0^1 

dag{Q^) = -lNs{N,ag)f^^^{N,Q% (4.1) 

where ^NsiN,a.s) = al{Q‘^)jNsiN,a.s)/f3{as). Similarly 


a2(Q2 


das{Q‘^ 


V 9{N,Q^) ) 


(jfg{N,as)\ (f^{N,Q'^)\ . . 

y%f{N,ag) %g{N,ag)J [ g{N,Q^) J’ ^ > 


with similar definitions for 7 jj(A^, CKg) etc. as for yjvs(iV, Os)• Each of the y’s may now be written 
as 

OO 

7(iV,a,) = ^7-’*(iVX+i, (4.3) 

n=0 

where we also have the analogous definition 


7(iV,a,) = ^7-’*(iVX+i, 

71=0 


(4.4) 


for the normal anomalous dimensions.® Thus, at zeroth order the 7 (A^)’s only differ from the 7 (A^)’s 
by the normalization factor bo. Beyond this order the differences are more complicated. 

Using our definition, the evolution equations may be solved order by order in ag- When 
doing this it is necessary to choose a starting scale Qq for the perturbative evolution of the parton 
distributions (or equivalently a starting value of the coupling ag{Qo)), and specify input parton 
distributions at this scale. Let us discuss the choice of this scale briefly. Qq must clearly be chosen 
to be large enough that perturbative evolution should be reliable and also such that higher twist 
corrections should be very small. Traditionally evolution only takes place up from this starting 
scale. This has been both for the simple reason of convenience, and also because some form of the 
inputs has been expected at low starting scales. Until a few years ago the requirements described 
above led to a choice of Qo ~ 4GeV^. In the past couple of years this value has often been chosen 

® We use the superscript I in order to denote that this is the standard as expansion of the anomalous 
dimensions rather than the expansion in (1.6), i.e. the leading-ln(l/x) expansion, which will be used more 
frequently in this paper. 
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to be rather lower, due to the apparent success of evolution from lower starting scales, and also 
because much of the interesting small-x data is now at < 4GeV^. These choices have been 
accompanied by guesses for the form of the inputs at low starting scales, e.g valence-like [9], or flat 
[4] [10]. However, with the quality of the most recent data, these guesses for the form of the inputs 
no longer seem to lead to good quantitative agreement with data. 

We make no assumptions of the above sort about the value of Qo- We do require it to be high 
enough to be in the perturbative regime and to avoid higher twists, but acknowledge that there is 
no reason why Qg cannot be chosen to be quite large, and evolution away from the starting scale 
performed both up and down in Q^. Taking this open-minded approach we then simply assume 
that only perturbative effects can lead to deviations from soft behaviour of the structure functions, 
and also demand that the form of our well-ordered expressions for the structure functions is as 
insensitive as possible to this choice, thus making the choice of Qo open as possible. We will 
see the consequences of this unusual approach to the input scale and the inputs for the structure 
functions as we progress. 

We begin by solving for the nonsinglet parton distributions, which are an easily understandable 
model. In this case the solution is particularly simple. Integrating both sides of (4.1) we obtain 








'-k=0 


a 


AQ 


a 






exp 


• oo 

E 

•n=l 




n 


(4.5) 


Perhaps unconventionally, we explicitly express the input as a power series in 

0(s{Qq)- There are two reasons why this is necessary. Firstly, changes in the starting scale Qo 
to Q!s((5o)“dependent changes in the expression for the evolution term, which must be compensated 
for by Q;s((5o)“dependent changes in the starting distribution in order to leave the whole expression 
for the parton distributions unchanged, as required.® Let us examine this briefly by looking at the 
change of the lowest-order piece of (4.5), i.e. 





(4.6) 


under a change in starting scale, Qq —> (1 + 5)Qq, where 5 is some constant. Each of the terms in 
our expressions for the structure functions would now be written in the forms shown above, but 
as functions of (1 + (I)(5o. We may regain expressions in terms of Qq by expanding the coupling 
constant as((l + ^)Qo) ™ the form 


a,((l + <5)Q®) = asiQl) “ 56o«"(Qo) + 0(«s(Qo))- (4-7) 

® Since the structure function is obtained by multiplying the parton distributions by the Qg-independent 
coefficient functions, the parton distributions must be Qo'independent to make the structure functions so. 
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Under this change in the input coupling constant, the evolution term in (4.6) undergoes a change 


Wm) 



+ higher order in as{QQ)- 


(4.8) 


This change in the parton distribution due to the variation in the leading term can be countered, 
up to higher orders, by a change in the order Us^Qq) input of the form 


(iV, Ql) = QD- (4.9) 

Because changes in the evolution term due to a change in Qq begin at first order in asiQo), 
are therefore absorbed by terms in the input at first order and beyond, the zeroth-order input 
must be insensitive to such changes and is QoUndependent: Qo) = /^,o(-^)- 

difficult to see that higher-order changes in the parton distributions due to changes in Qq can all 
be accounted for by changes in the higher-order inputs and that in order to satisfy this constraint 
alone, it is consistent for the higher-order inputs to be equal to functions of N dependent only on 
the anomalous dimensions. Therefore, these higher-order inputs are in some sense perturbative, 
consisting of perturbative parts multiplying the fundamentally nonperturbative As we 

will soon discuss, there are other constraints to be satisfied, e.g. factorization scheme independence 
of the input for the structure function, and this slightly complicates the above picture, but does 
not change the main conclusions. 

Thus, we see that it is necessary to express the input as a power series in OsiQ^) in order to 
make the parton distribution Qg-independent, but that only one intrinsically nonperturbative input 
which is Qg-independent is needed. Usually in analyses of structure functions the parton inputs are 
taken to be a single as((5o)Andependent function which is implicitly allowed to be Qg-dependent. 
Phenomenologically, this is normally much the same, but we stress the formally correct expression 
for the input here since it is rather important when constructing properly ordered solutions, and 
leads to some predictive power, especially in the small x limit. 

Another reason for explicitly writing the input as a power series in asiQ^) is that it makes little 
sense to demand that the starting distribution should not have a perturbative expansion, unless 
one insists that there is something special about a particular factorization scheme. If there is not, 
then any transformation from a scheme in which the starting parton distribution is purely zeroth 
order in as(Qg) will lead to a starting distribution with a power series expansion in as{Qo), but 
again, where there seems to be some underlying nonperturbative input which is unaltered 

by the change in factorization scheme. 

Accepting that the parton inputs should be a power series in agiQ^), and substituting our 
solution for the parton distribution into (3.3), we obtain the general expression for the nonsinglet 
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structure functions 


F,^^iN,Q^) = 


Nf 


m=l 


j=l k=0 
0,1 


k2 + E 


a,{Qp\^-sW 

OisiQ^ 


exp 


E 

•n=l 


n 


r/siN) 

(4.10) 


It is clear that this may be written as 

'OsiQ, 


fP^{n,qP = 


2\ \ 4^s(^) “ EL 


OsiQ"^ 


E E Fp±{N,Ql)ar^{Ql)aP{Qp, (4.11) 


n=0 m=0 


and that, once a choice of renormalization scheme and starting scale have been made, each of the 
Fnm{Ff) must be invariant quantities under changes of factorization scheme in order to guarantee 
the scheme independence of the whole structure function. However, it is also clear from (4.10) that 
each of the FP^{N) is potentially a function of the starting distribution, the anomalous dimension 
and the coefficient function. Any well-ordered calculation of the structure function should include 
all complete terms in (4.11) up to a given order in n and m, and no partial terms. In practice it is 
possible to work to a given order in n including all m <n, i.e. to expand to a given order in powers 
of as{QP plus powers of Q!s(go)) if fh® y’s and Q’s are known to this order. (It is impossible to 
work to a given value of m including higher values of n — m without knowledge of 7 to order n — m, 
i.e. it is impossible to expand just in powers of as(g^).) 

We shall briefly describe how to construct this ordered solution for the structure functions, 
working up from zeroth order, so that at each order the solution is factorization scheme dependent. 
Consider first calculating Fp^{N, Q^) by working to zeroth order in CP^{N,as), ^Ns{N,as) and 
the starting distribution (remembering that this is gg-independent). To this order 




asjQl 

as(Q^ 




and therefore 


Nf 


i=i 




(4.12) 


(4.13) 


Using the one-loop expression for the running coupling, as is appropriate for a lowest-order cal¬ 
culation, each of the quantities in this expression is factorization scheme independent and indeed, 
also renormalization scheme independent. Therefore we have a consistent leading-order (LO) ex¬ 
pression. 

If we calculate FP^{N, g^), again by combining the coefficient functions with the solutions for 
the parton evolution, then working to the order n = 0 in all quantities leads to FPq{N,Q^) = 0 , 
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since the zeroth-order coefficient function is zero. However, looking at the expressions for the 
general solutions, (4.10) and (4.11), one sees that the only contribution for n = 1 in (4.11) comes 
from working to first order in Cf^^{N,as) and to zeroth order in j]ys{^,Ois) and the starting 
distribution. This leads to the LO expression for of 


(iV,Q=^) = asmClUN)'yejf^^^,{N) 


NS 


Nf 


1 = 1 


a 


m 


a 






(4.14) 


Again, using the one-loop expression for the running coupling, every term in this expression is both 
factorization scheme and renormalization scheme independent, giving a well-defined LO expression. 

We now consider the first correction to these expressions. The first-order expression for the 
renormalization group equation is 


««(Q") 


das{Q'^) 


-7^'5(iV,a,s)/^;5(iV,Q^) - a.(Q")7Al5(^)4,o 


-iJ 


iN,Q^ 


(4.15) 


with solution 


= (asiQl) - 


ij 


cNS, 


cNS, 


a 


M 


a. 


■XQ^ 


7^s(^) 


(4.16) 

Multiplying by as{Q‘^)CX^^f^i{N) and adding to we clearly obtain all terms in the expres¬ 
sion (4.11) for {N,Q‘^)) at n = 1. Adding this to (4.13) we obtain the factorization-scheme- 
independent expression for FX^^{N, Q^) up to n = 1 . 

We note that this is not the same as finding the complete solution to the renormalization group 
equation including all terms in the anomalous dimension up to first order in 7 and multiplying the 
solution by the coefficient function up to first order. This procedure would involve the exponenti¬ 
ation of the anomalous dimension, and thus would include incomplete parts of the F,{^(A")’s for 
n > 1 , and would be a factorization-scheme-dependent, and hence physically ambiguous quantity. 

Similarly to F^^{N,Q‘^) we can obtain the NLO factorization-scheme-independent expres¬ 
sion for FIX^{N, Q^). The expression at n = 2 is obtained by adding the product of the first-order 
coefficient function and to the product of the second-order coefficient function and 

/^q(A', Q^). The NLO FXX^{N,Q'^) is the sum of this and (4.14). When working to NLO for 
either structure function we now have expressions which are renormalization scheme dependent. 
This scheme dependence compensates for the renormalization scheme dependence of the two-loop 
coupling constant (which has a renormalization-scheme-dependent value for Aqcd), and it is this 
expression for the coupling that we should use at this level. Doing so guarantees the renormaliza¬ 
tion scheme independence of the structure functions up to corrections of higher order in a^, i.e. 
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It is now simple to see how to construct factorization-scheme-independent structure functions 
order by order. Defining the nth order renormalization group equation by 


OisiQ^) 


dasiQ^) 


m=0 


(4.17) 


it is easy to prove by induction that the solution contains all terms in the full solution with a given 
sum of powers of as{Q'^) and asiQ^) multiplying the everpresent {as{Qo) factor. 
Thus, defining part of (4.11) by 

asiQo 

OisiQ'^ 

is given by 






Fi^^UN)ar^{Ql)aT{Q^) 


(4.18) 


m=0 


Nf n 

Fl!rfiN,Q^) = E"' E (4.19) 

j m=0 

The nth"Order scheme-independent term in the expression for the structure functions is given 
by the product of the mth"Order coefficient functions and the (n — m)th^order solutions to the 
renormalization group equation summed over m. The nth^order scheme-independent structure 
function is then the sum of these terms up to order n. We must however remember that including 
all F^^{N, Q'^) up to order n is working to (n + l)th nontrivial order, and requires the (n + l)-loop 
coupling in order to make the expression renormalization scheme invariant up to higher orders in 
a*. Similarly, including all F^^{N,Q‘^) up to order n is working to nth nontrivial order, and 
requires the n-loop coupling. 


This procedure clearly provides factorization scheme independence and renormalization scheme 
independence for this method of expansion. We can also see how it relates to the discussion of the 
factorization-scheme-invariant evolution equations in terms of the structure functions. In order to 
do this let us consider the solution for the non-singlet structure function F^^{N,Q^) again. We 
may rewrite our general solution (4.10) in the form 



exp 


E 

•n=l 


«(Qo) - <(Q^)) ~n,l 


n 


r/siN) + 


r-ln 




d 

dlnq'^ 


\n(l+Y,C^AN)c^TiQ"))dlnq^ 

m=l 


(4.20) 

This way of writing F^^{N, Q"^) is particularly useful since it separates the solution into the value 
of the structure function at Qq (the term in square brackets), and the ratio of its value at a 
different with this initial value (the rest of the expression). Clearly these two quantities must 
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be separately factorization scheme independent. Also, it is clear that this solution is of exactly the 
same form as (3.8), and we can express it simply in terms of an input for Q^) at Qq and a 

physical anomalous dimension which governs the evolution, i.e. our input and our prediction. Since 
these are both physical quantities they must be separately renormalization scheme independent. 

The input and the evolution will mix with each other if we make a change in the starting 
scale, however, as with the parton distribution. Examining the effects of such a change for the full 
physical quantity gives us information about the form of the input. As we have already seen, at 
lowest order the nonsinglet structure function is just the sum of the charge weighted nonsinglet 
parton distributions. So the lowest-order input is just 

Nf 

= (4.21) 

i=i 

Making the change of starting scale and consequently of as{Qo) already considered, the change in 
the lowest-order evolution is as in (4.8), and this leads to a change in the lowest-order structure 
function which is of higher order, and which can be absorbed by a change in the NLO input for 
the structure function of 

AFi^f{N,Ql) = 6boj°JsiN)Fi!o{N). (4.22) 

In terms of parton distributions and coefficient functions 


Nf 


2,1 


{N, Ql) = ^ e2(/^^_^(iV, Ql) + 


(4.23) 


i=i 


We chose the change in in (4.9) so that the structure function would be independent 

of Qq, and it is clear that that is consistent with (4.22) and (4.23). However, we can now say 
more about the form of /^‘^(iV, Qq). Because it is the leading term in the expression for the input 
which depends on a<j((5g), F^f {N,Qq) must be renormalization scheme independent. However, 
is renormalization scheme dependent, so Q^) must also be renormalization scheme 

dependent in a way such as to cancel this. Hence, must not only have a part like 

ln(( 5 o) 7 Ar 5 (A')/,^,o(A^) ill order to maintain Qg-independence of the structure function, but also a 
part like —C'^i^j(A^)/^o(A^) in order to maintain renormalization scheme independence, i.e. 




0,1 


pNSi 


(4.24) 


where A^s is some unknown scale parameter. So we see that 


Ql) = HQl/AMs)l^jtsi^)F^,i{N) ^ \n[Ql/Aj,s)T^N'siN)F^,i{N). (4.25) 

It is clear that this does not spoil our argument that (A", Qg) consists of perturbatively calcula¬ 
ble quantities multiplying the fundamentally nonperturbative It is also clear that all the 
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higher-order inputs may be chosen to be perturbative functions multiplying this nonperturbative 
input, and therefore that the input for the structure function is a perturbative power series (depend¬ 
ing on the physical anomalous dimension) multiplying the single nonperturbative factor 
which from (4.21) may be interpreted as a fundamentally nonperturbative input for the structure 
function F^(f{N). Hence, demanding invariance of our expression for the structure function under 
changes in Qq leads us to a power series expression for the input, but with only one (for each 
quark) really nonperturbative factor for this input. We also see that if Qq = A^s, the first-order 
perturbative correction to {N,Qq) vanishes. Hence, we might expect A^s to be some scale 
typical of the transition between perturbative and nonperturbative physics, i.e A^s ^ IGeV. 

Separating the expression for the structure function into a definite input and evolution part 
also enables us to view the loop expansion in an alternative manner. We see that when expanding 
out to nth order in the loop expansion we are including all terms where the order of the input part 
for the structure function added to the order of the evolution part of the structure function is less 
than or equal to n. This is clearly the same as including all powers of ag in the expression for the 
structure function up to nth order, but gives us some additional physical interpretation. 

Writing the solution as in (4.20) does, however, also illustrate that demanding factorization 
scheme invariance does not on its own force us into the strictly defined loop expansion. It is clear 
that we could, if we wished, expand the input and evolution term out to different orders in 
still maintaining factorization scheme independence. However, this is not a sensible approach for 
reasons of renormalization scheme dependence. If we were to expand out the input and evolution 
terms to different orders in we should really use ag itself calculated to a different order in each 
case, surely a perverse thing to do. Also, when using the resulting expression to evaluate the 
structure function at some away from we would only have a subset of the terms at some 
given power of Ug (where Ug may represent q;s((5o) or ag{Q'^)). Each of these terms is presumably 
of similar magnitude, which is rather importantly equal to, or greater than, the magnitude of the 
uncertainty due to renormalization scheme dependence. This is even the case even if we simply use 
the product of the input specified to nth order and the evolution calculated to nth order. In this 
case there is no problem in what definition of ag to use: we simply use the n-loop expression for 
the input and evolution and hence the n-loop coupling for both. However, the resulting expression 
contains, for example, a term like ag“^((5o)(cts~^(f5o) “ but none like 

(Qo) ~ Q^) which should be of the same size. One only 

has a complete set of terms of order up to A: = n — 1. Moreover, this type of term is higher 
order than the overall uncertainty in the expression due to the renormalization scheme dependence, 
which is of order a^F^Q {N,Q‘^). Exponentiating the solution for the evolution part, once it is 
found to a given order, is also redundant, since that would introduce only a subset of the terms at 
higher orders in ag for the evolution, and these terms would be renormalization scheme dependent. 
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Thus, all ways of obtaining factorization-scheme-invariant expressions other than the loop 
expansion contain terms which are essentially redundant using this expansion method, i.e. those 
beyond the order where all complete parts of the loop expansion exist. These other terms in the 
expansion are the same order as uncertainties due to renormalization scheme dependence. The 
loop expansion gives the minimum expression which is totally factorization scheme independent, 
and moreover, renormalization scheme consistent (forgetting complications due to powers of 
for the present) to a given order. If solving using the factorization-scheme-invariant equations in 
terms of physical variables it is useful to bear this in mind.^*^ 

The solution for the longitudinal structure function is much the same as for F^q{N, Q^). We 
may write as 
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\n(j2cZi,i{N)aT-\q^))dW 


m=l 


(4.26) 

So as with F^^{N,Q^) the full solution in terms of parton densities, anomalous dimensions and 
coefficient functions is easily separated into its input and evolution parts and therefore directly 
compared with the solution using the effective physical anomalous dimensions (where if we solve 
using the Ts, the solution is expressed in terms of an input and evolution part once we write the 
overall power of as{Q‘^) multiplying F^^^{N,Q‘^) in the form as{Ql) x {as{Q'^)/as{Ql))-) Once 
again, from (4.26) we see that the loop expansion multiplies orders in the input by orders in the 
evolution in a systematic manner. 


The situation for the singlet structure function is similar to that for the nonsinglet, but is 
complicated by the fact that we now have coupled evolution equations for the quark and gluon 
distributions. This makes it impossible to write a closed form for the solution to the renormalization 
group equations in the way we did for the nonsinglet case in (4.20) and (4.26), but the equations 
may be solved order by order in the same way. 

It is also true that, when solving using the parton model, other ways of obtaining factorization- 
scheme-invariant expressions other than the loop expansion would need, in comparison, more complicated 
prescriptions. 

Comparing (4.26) with (4.20) we see that there are no new fundamentally nonperturbative parts in 
the input. Hence, the complete expression may be interpreted as a perturbative part multiplying the 
nonperturbative factor which is F'^q{N). It is only the perturbative factor which is different in the two 
nonsinglet structure functions. 
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The lowest-order solution to (4.2) is 





9o{N,Q^) 


9^{N) 


O^siQl) 

as{Q^) 


j^iN) 

+ 9o (^) 


asjQl) 




(4.27) 


where and ^^\n) are the eigenvalues of the zeroth-order matrix for the 7 ’s, and fQ’^{N) + 

/(f’“(A^) = foiN) and g^iN) + gQ{N) = go{N). Similarly to the nonsinglet case, the lowest- 
order factorization-scheme- and renormalization-scheme-independent expressions for the structure 
functions are 

Fi{N,Q^) = fiiN,Q^), (4.28) 


and 

FiiN,Q^)=asiQ^){Cl,{N)f^iN,Q^) + Cl,iN)goiN,Q^)). (4.29) 

Under a change in starting scale, the change of the input for F^iN, Q^) begins at order as{Qo), and 
that for FI{N,Q‘^) begins at order o^Qq). Therefore, as with the nonsinglet quark distributions, 
the lowest-order inputs for the partons are Qo^independent. Also, if we expect any powerlike 
behaviour to come only from perturbative effects, then these QQ-independent inputs for the quark 
and gluon are analytic for > 0 . 

Again in clear analogy to the nonsinglet case, the first-order expression for the renormalization 
group equations is 

« d //f(iV,g 2 )\ ^ fff(N,Q^) 

Ua.iQ^) [9 i{N,Q^) J Wj{N,as) 7^’‘(N,a,)J [9i(N,Q^) 

(4.30) 

and the solution is 



(4.31) 

and similarly for gi{N, Q^), where Ql) + /f “(iV, Ql) = fi{N, Qg). 
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This is clearly of the same form as (4.16) and, as in the nonsinglet case, the n = 1 term in the 
factorization-scheme-independent expression for F 2 {N,Q^), including all single powers of as{Q‘^) 
or as{Qo) multiplying {a,siQo)/asiQ'^)y+’-^^\ is 


/f (iV,Q2) + asiQy{ClyN)fi{N,Qy + C£,(iV)5o(iV, Q^)). 


(4.32) 


In the same manner the n = 2 term in the factorization-scheme-independent expression for the 
longitudinal structure function is 


«.(Q")(C'i,i(iV)/f(iV,Q") + C'£i(iV)ffi(iV,Q")) + a^(Q2)(c{^2(iV)/o^(iV,Q2) + c'£^2(iV)5o(iV,Q2)). 

(4.33) 

As for the singlet structure function, there is a simple prescription for calculating nth“Order 
factorization-scheme-independent structure functions. Defining the mth“Order solution to the 
renormalization group equation as the solution to 


a.(g2) 


d 

das{Qy 


\gm{N, g2) J 


m 


fc =0 


f7^f(N,ay 7%{N,ay\ 

\7gf{N,ay 7^g{N,as)) 


the nth“Order term in the expression for the structure functions is given by 


(fy.yN,Qy\ 

ygrn-k{N,Qy J ’ 

(4.34) 


= Y, aT{Qy{Cl^{N)ft^{N,Qy + Cl^{N)g^_^{N,Qy,) (4.35) 

m=0 

and the nth~order structure function is the sum of all such terms up to n. The comments concerning 
the order of the renormalization scheme independence and the order of the coupling constant to 
be used for F 2 {N,Q‘^) and Fl{N,Q'^) for the singlet case apply in exactly the same way as for 
the nonsinglet case. Once again, the loop expansion is a well-ordered way to expand the structure 
functions. 


We may now discuss the relationship to the solutions using the evolution equations for the 
structure functions. First, we can explain in more detail why we feel it is appropriate to take out 
the power of as{Q‘^) from the longitudinal structure function when using this approach. Using the 
parton model the lowest-order expression for F 2 consists of two g^-i^dependent factors, 

the inputs for the structure function, multiplying (agiQ^)/as{Q'^)y+’- ■ Similarly, we may think of 
the lowest-order expression for F^{N, Q^) as being two inputs for the structure function multiplying 
the g2-dependent factors {as{Qo)/as{Q‘^)y+’-~^, with the extra power of asiQ"^) coming from the 
coefficient function. So, at lowest non-trivial order the evolution of the two structure functions 
is different. If we were to solve the evolution equations for the structure functions themselves at 
lowest order, then we would do it in the same way as for the parton distributions: obtaining the 
g2-dependence from the eigenvalues of the lowest-order anomalous dimension matrix, and the 
inputs multiplying the two evolution terms from F 2 {N, Qq) and F^{N, Qq) and the eigenvectors of 
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the matrix. Thus, the terms governing evolution would have to be the same for both structure 
functions, which is clearly in contradiction with our lowest-order result using the parton model 
(also, the lowest-order input for Ff {N,Q‘^) would be zero). This leads to the solutions for the 
structure functions being built up in a very different way using this method of solution from that 
using the parton model. 

Hence, we choose not to forget the success of the parton model, and let the longitudinal 
structure function be multiplied by an overall power of agiQ"^), as discussed earlier. Working with 
the r(A^, as)’s in (3.12) and solving the lowest-order evolution equations, one then obtains precisely 
the same lowest-order solution for the structure functions as above. It is not difficult to see that, 
although the T^’^’s are not identical to the the eigenvalues of the anomalous dimension matrix 

are the same. Thus, we obtain the factors of {otsiQo )™ these solutions as we would 
hope, and get the extra factor of {aa{Q‘^)/as{Qo)) in the expression for F^{N, Q^) from the overall 
power of UsiQ'^) multiplying the solution (and simultaneously get the extra power of as{Qo) in the 
input). The eigenvectors are different from those in the parton model, of course, but this takes 
account of the fact that the inputs are now those for the structure functions themselves, rather than, 
in the case of Ff the parton densities weighted by coefficient functions. Working beyond 

this lowest order and finding the nth“Order solutions to the evolution equations for the structure 
functions, in the same way as those for the parton densities, as outlined above, leads to all powers 
of as up to power n multiplying the factors of (as((5o)/Q^s(Q^))^^’^ ^2 powers of 

Us up to power n + 1 multiplying Thus, it produces the 

same result as the nth^order solution using the parton distributions and the loop expansion. 

Unlike the nonsinglet case it is rather difficult to see how to express the solution in terms of 
factorized inputs and evolution parts; the evolution parts will not be simple exponentials, as in 
(4.20). However, the solution for the structure functions can be written as a series of terms each 
of which can be factored into a part dependent only on UsiQ^) and one depending on both UsiQ^) 
and asiQ"^), the latter either vanishing at = QE or being equal to {as{QQ)/as{Q‘^)y+’-^~^^ 
and thus equal to unity at = QE Tor each of these terms the former part may be interpreted 
as an input and the latter part may be interpreted as the evolution. Within the loop expansion we 
then include all terms where the order of the input part plus the order of the evolution part sums 
to less than or equal to some integer n. 

By examining the form of the inputs under a change in starting scale, as with the nonsinglet 
structure functions, we find that the only fundamentally nonperturbative parts of the inputs are 
the zeroth-order parts, F 2 q{N) and Ff ^{N), with all other inputs being in principle expressed as 
perturbative functions of the physical anomalous dimensions multiplying one of these nonperturba¬ 
tive components. If we wished to stay within the language of partons then we would equivalently 
have the two nonperturbative parton inputs /^{N) and go{N), where the link between the two 
pairs of nonperturbative inputs is 

Fi,o{N) = /o^(iV), (4.36) 
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and 


= CliiN)fi{N) + Ci,{N)go{N). (4.37) 

We choose to think of the expressions for the structure functions as the real nonperturbative 
functions, since they are, of course, the physically relevant quantities.This is consistent with the 
reasoning that go (N) is (5o”independent because the invariance of the physical quantity Ff {N, ) 

under a change in starting scale leads to q{N) being Qo”independent, which then, using (4.37), 
leads to go{N) being Qg-independent. 

The extra complexity of the solution in the singlet case, when compared to the nonsinglet case, 
also means it is far from obvious how to express the physical anomalous dimensions in terms of the 
parton anomalous dimensions and coefficient functions simply by comparing the forms of the solu¬ 
tions in the partonic language and purely in terms of structure functions. However, the expressions 
for the physical anomalous dimensions can be found at each order using Catani’s expressions (3.10). 
One could then solve for the structure functions in a manner different from the loop expansion. 
An example would be to calculate the solution to the whole evolution equation for the structure 
functions using the physical anomalous dimensions up to a given order. This would be by definition 
factorization scheme invariant, but would only be renormalization scheme invariant up to the same 
order in the solution as the order of the anomalous dimensions. The rest of the solution would 
contain only a subset of the possible terms of a given form obtained from the full calculation, and 
would therefore still have no real signihcance. Hence, we conclude that the evolution equations in 
terms of structure functions do nothing to alter the strict ordering of the solution using this method 
of expansion, but do make finding this ordered solution somewhat easier. 

This whole discussion of the order-by-order-in-Os expansion scheme is perhaps a little aca¬ 
demic since the errors invoked in performing the calculations without paying strict heed to the 
formally correct procedures are rather small. For example, let us consider factorization scheme de¬ 
pendence. When using this standard method of expansion only two different factorization schemes 
are generally considered, the MS scheme and the DIS scheme. The first of these simply involves 
calculating the coefficient functions using standard techniques and using the MS scheme to remove 
both the the infrared and ultraviolet divergences and hence provide the anomalous dimensions. 
The second makes a change of parton distributions, so that the singlet quark distribution is equal 
to the singlet structure function. It is clear from the above discussion that, if a calculation is made 
in one scheme to a well-defined order (usually next to leading order), and then a transformation 
to the other scheme made correctly, precisely the same result will be obtained. In practice, small 

We could, if we wished, make a change in factorization scheme where the zeroth-order part of the 
matrix U was not the unit matrix, but a function of N. This changes no physical result, but does change 
the zeroth-order parton inputs. 
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differences are sometimes noticed between calculations using different schemes within the loop ex¬ 
pansion to NLO, but these come from well-understood sources. One common source is that the 
starting distributions in both schemes are described by a simple functional form, e.g. 

xfix, Ql) = AiQl)x-^^Qo)^i _ + e(g 2 )^V 2 + (4.38) 

rather than the formally correct expression of a power series in as{Qo) with essentially perturba¬ 
tive coefficients convoluted with a Qo^i^dependent nonperturbative function of x. If the starting 
distribution is of the form above in one scheme, then the starting distribution in the other scheme 
will not be modelled precisely by a function of the same form. However, the error is in general 
small. Alternatively, if the calculations are not done in a well-ordered manner, e.g. simply solving 
the whole evolution equation using the two-loop anomalous dimensions, then differences between 
the calculation done in the two schemes, or between this type of calculation and the correct NLO 
calculation, will be of NNLO. Again this usually results in only small differences.^^ Similarly, small 
differences would be obtained by working in two different renormalization schemes, and of course, 
there is nothing which can be done about this. 

There are a couple of points we wish to make here. Firstly concerning the form of the input in 
(4.38) we note that, if one considers the input to be consistent with the loop expansion, e.g. NLO- 
order evolution should be accompanied by a NLO input, the power of A(gg) should not correspond 
to parton distributions much steeper than flat for the singlet quark or gluon at this order. This is 
because the first”Order-in~Q;s(gQ) input should be accompanied by no more than a single power of 
ln(l/x), this being all that is required to absorb the change in the zeroth-order evolution under a 
change in starting scale at this order in agiQ^)- In practice restricting A(gQ) in this way is rather 
important for the fits to the low x data, and would mean that NLO fits to small-x data would be 
very poor. The only way to avoid this is to let A(go) be an artificial free parameter, in which case, 
if (4.38) describes the singlet quark density, it must be ~ 0.2 — 0.3 for practically any Q^. This 
value is totally unjustified, and the need for this steepness in the input for the quark is a clear sign 
of the limited usefulness of the NLO-in-as calculation at small x. We will see that the situation is 
more more satisfactory when using the correct expansion method. 

Ignoring this problem with the inputs, there is another reason for being concerned about the 
validity of the loop expansion at small x. The main reason for the smallness of the differences 
between inaccurately performed NLO calculations noted above, even at small x, is that the differ¬ 
ences between these calculations do not contain terms which are any more leading in ln(l/x) than 

For a comparison of calculations done at NLO using different methods see [34]. One sees here that 
starting with the same input parton distributions the different calculations agree to better than 5% except 
for very small x and large , where the discrepancy can approach 10%. The differences observed can be 
compensated for almost completely by small changes in inputs. 


27 



the NLO calculation itself. Hence they are genuinely an order of ag down on the NLO calculation, 
with no small-x enhancement. This is also true for calculations done in different renormalization 
schemes (or at different renormalization scales). The real NNLO contribution will be higher order 
in as but will also contain terms at higher order in ln(l/x), and will therefore be potentially large 
at small x. We therefore stress that the relative insensitivity of structure functions to changes in 
renormalization or factorization scheme for calculations which are not carefully ordered is no guar¬ 
antee that genuine higher-order corrections will be small at small x when using the loop expansion. 
Indeed we would naively expect them to be large. This is a point to which we shall return. 

In contrast to the insensitivity when using the loop expansion, when using the leading-ln(I/x) 
expansion very large differences between calculations done in a large number of different factoriza¬ 
tion schemes have been noted. This is an obvious sign that the calculations are not being done 
in a well-ordered manner, and that the ambiguity introduced by lack of care in the calculations is 
greater in this method of expansion than the standard loop expansion. We will now demonstrate 
that this is indeed the case. 

4.2. The Leading— ln(l/x) Expansion: Parton Distributions. 

It should clearly be possible to define a well-ordered, and hence, factorization scheme indepen¬ 
dent expansion in leading powers of ln(I/x), or equivalently, in leading powers of 1/N in moment 
space. We will now demonstrate that this is indeed the case. As in the loop expansion, we will first 
work in terms of the traditional parton distribution functions and coefficient functions, and see how 
this results in expressions containing the physical anomalous dimensions. Doing this enables us to 
see how large factorization scheme dependence can arise when calculating less carefully within this 
expansion scheme. 

First we must set up our notation and qualify the statements made in the introduction con¬ 
cerning the leading-ln(l/x) expansion. We stated that the anomalous dimensions could be written 
in the form (1.6), and this means that using the form (4.1) and (4.2) for the evolution equations 
we may write, 

OO OO OO 

7 (a^, q^) = Y, <iQ") E r^oc 7 iQ")N-^ ^ E cxiiQ")rMQ^)/N). (4.39) 

71=0 m=l —n n.=0 

So, the 7 *^’s only differ from the 7 *^’s by the normalization factor of bo, but at higher orders the 
relationship is more complicated. In particular, in the MS renormalization and factorization scheme 
we may write 

OO 

%giN,as{Q^)) = Y<{Q^)TggMQ^)/N), (4.40) 

n=0 
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where the series expansion for jgg{as/N) is known to all orders (all the coefficients being positive). 
This expression for Jgg{cis/N) is renormalization scheme independent since any change in renormal¬ 
ization scheme can only bring about a change in the coupling at O^a^), and this does not require a 
change in Jgg{cis/N) to keep physical quantities invariant at leading order in this expansion scheme, 
as we will soon see. ^gg{as/N) is of course renormalization scheme dependent, since it must change 
to absorb part of the effect of the 0{a‘l) change in the coupling on jgg(as/N). The renormalization 
scheme independence is also true for ^gj:{as/N), which obeys (1.7), and also for ^^j:{as/N) and 
^^g{as/N) which are both zero. There is also a renormalization-scheme-independent relationship 
between jjj:{as/N) and ^jg{cis/N), which tells us that 





GnboJ 


(4.41) 


where the second term in the brackets is the one-loop contribution to 7 j^(a<j/iV). It is also known 
that the nonsinglet anomalous dimension has no poles at = 0, and neither does the nonsinglet 
coefficient function. Hence, the nonsinglet sector makes very little contribution to the structure 
function at small x, and as such we will ignore it for the remainder of this section. 

A general change in factorization scheme may be expressed by writing an element of the 
transformation matrix U as 


OO OO 

Ua,{N,as{Q^))= Y, E UT<{Q")N- 

n= — k m=l — n 


Y a:{Q^mas{Q^)/N), (4.42) 

n= — k 


with condition on the 17™ such that U obeys Uab = <^ab + Cl(«s)- This flexibility in the factorization 
scheme means that all the above results on the low-order y’s are in principle factorization scheme 
dependent. In particular, the expressions for the y’s in (4.39) may be insufficient in general, and 
should be replaced by 

OO OO 

jiN,Q^)= Y <iQ^) E «nmar(Q")iV—. (4.43) 

n= — k m=l — n 

However, this is only necessary if one makes a change of scheme away from a standard scheme, 
such as the MS scheme, using a matrix U containing terms such as where m > 

n. For simplicity and because, of course, all physical results are ultimately factorization scheme 
independent by definition using any expansion, we will forbid such singular changes of scheme and 
demand that A: = 0 in (4.42). 

However, a change of factorization scheme with A = 0 can still introduce scheme dependence 
into the 7 j(^(as/A^)’s. Using (2.12) we see that 

labiots/N) = E E (4.44) 

c d 
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Insisting that jgg{as/N), ^^^{ag/N), jjj:{as/N) and ^^g{as/N) are unaltered by scheme changes 
leads to the requirement, 


U^fg{N,ag) = 0, t/°/iV,a,) = l, 

oo 

U^ggiag/N) = l + Yl U^g^{agiQ^)/Nr, Ulf{oig/N) = liU^g{ag/N) - 1). 


(4.45) 


n=l 


This requirement also preserves the relationship (4.41) between jjj:{ag/N) and ^jg{ag/N). Again 
for simplicity, and due to factorization scheme invariance of physical quantities, we will only consider 
factorization scheme changes away from MS scheme of the type (4.45). The UJ^f^{aa/N) for n > 0 
will have no restrictions.^^ 

Restricting ourselves to these schemes we may now write 

OO OO OO 

C|(iV,a,) = l +J;<(Q2) ^ Cl^^^aT{Q^)N-^^l + Y,o^:{Q")Cln{as/N), (4.46) 

n=l m=l —n n=l 

and all other coefficient functions as 

OO OO OO 

C“(iV, a,) = ^ <(Q2) ^ ^ Y. <{Q^)ClMslN). (4.47) 


n=l 


m=l —n 


n=l 


All the Cf^{as/N) are both renormalization-scheme- and factorization-scheme-dependent quan¬ 
tities. Indeed, all of the ^ are renormalization scheme and factorization scheme dependent, 
except for the Cf^ ^ which come from the one-loop longitudinal coefficient functions which, as 
we saw in the previous subsection, are totally scheme independent. There are also two renormaliza¬ 
tion and factorization scheme (with our restrictions) independent relationships between coefficient 
functions: 

Cl,{ag/N) = ^[clAo^g/N) - (4.48) 

where the second term in brackets is the one-loop contribution to C 2 i{ois/N), which is itself renor¬ 
malization scheme and factorization scheme dependent (being equal to (Nf/Gn) in MS scheme); 
and 

(cUa./N) - T) = l(c|,.(a,/JV) - ^). (4.49) 

where the second terms in the brackets are the one-loop contributions to C^^{aa/N) and 
Cl ^{ag/N), both of which are renormalization and factorization scheme independent. 

To our knowledge, no-one has yet bothered considering the effects of any transformations of the type 
we have forbidden (of course, there is no point, since physical quantities are factorization scheme invariant). 
A number of the type described above with Ugg{N,as) ^ 0 have been considered (however, once again, 
physical quantities are invariant under such scheme changes). 
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Working in an arbitrary factorization scheme (up to the above restrictions) and using the 
general expressions for the 7 ’s and coefficient functions, we may find expressions for the structure 
functions. The first step towards this is solving the renormalization group equations for the parton 
distributions. The lowest-order part of the equation is, 

2 .^ 2 ^ d ( 0 0 W/d'(iV,Q2)\ 

>da,{Q^) \goiN,Q^) J - ) [go{N,Q^) J ' 

This may easily be solved to give 



go{N,Q'^) = {go{N,Ql) + |/o^(iV, Qg)) exp 




U as{Q^) 


a 2 (g 2 


das{q^ 




(4.51) 

This is analogous to the lowest-order solution within the loop expansion and contains two 
factors, one of which must appear in all the higher terms in the expansion: instead of 
{cuaiQ'^)/ois{Q‘^)y+-^^\ corresponding to the two eigenvalues of in the loop expansion, we 


have exp 




and 1 , corresponding to the two eigenvalues of ^^{ag/N) 


in the leading-ln(l/x) expansion, i.e. ^^g{as/N) and 0. In particular we notice that fQ{N,Q‘^) is 
-independent. 

Let us briefly digress in order to discuss the form of the inputs for the parton distributions. 
In particular, we examine how the terms in our expression change under the change in input scale 
Qg ^ (1 + S)Qq, leading to the change in coupling at the starting scale (4.7). Rather trivially, our 
expression for fQ{N,Q^) = /^(A'", Qg) is unchanged by the change in the coupling at the starting 
scale, and as such can be chosen to be independent of Qq. Hence, as for the lowest order input in 
the loop expansion we may write /if(A^, Q^) = /if(A^, Qg) = /^{N). The expression for the gluon 
involves a little more work. The change in the evolution term under a change in starting scale is 


A exp 


('20)70^(0-^ (g,2)/JV) 


'aRQ 2 ) 

= exp 


das{q^ 


a 2 (g 2 ) 

as((i+<5)Qo) ^^g[as{q^)/N) 


L4qs(Qo) 


a 2 (g 2 


das{q^)) - 1 


exp 


and using (4.7) we find that 


(Qg)^o^(Q,^(g2)/jv) 




a 2 (g 2 


das{q^) 
(4.52) 


A exp 


Wo) 7 °g(a«(g 2 )/iV) 


VJ q;s(Q^) 


a 2 (g 2 


das{q^ 


= {e^p[-6bo^gg{as{Qo)/N)) - l]exp 
+ higher order in ag and/or N. 






Q; 2 (g 2 


dag{q^ 


(4.53) 
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We choose 5'o(-^)Qo) the change in exp dag (g^) can be compensated 

for by a change in goi^^^Qo) up to corrections of higher order. Hence, the gluon input may be 
written as 

9oiN,Ql)=goiN) + (goiN) + lfi{N)) ^ = go{N) + go{N,Ql). (4.54) 

The change in go{]^,Qo) under a change in starting scale necessary to make go{N,Q‘^) invariant 
under changes of Qq up to higher orders may easily be calculated from (4.53), and is equal to 

Ago{N,Ql) = 5bo^ggias{Ql)/N){goiN) + ^fo{N) + 50 (A^,Qo)) + higher order 

(4.55) 

= ^lgg{<^s{Ql)/N)igo{N) + |/o^(lV) + 5 o(A^,( 5 o)) + higher order, 

in the limit of small 5. Thus, an appropriate expression for the input which satisfies this condition 
is given by 

{9o{N) + lfi{N)+~go{N,Ql)) 

= igo{N) + ^f^{N))(l + ^ exp[ln(Qg/Hgg) 7 °g(a,(Qg)/iV)], 

k rn=l ^ ^ ^ 

(4.56) 

where Agg is an unknown scale. The series Ylm=i 9 o,m{o:s{Qo)/N)'^ is at yet undetermined, but is 
potentially renormalization and factorization scheme dependent. It will only be determined when 
we come to construct the structure functions themselves. 

We also consider the form of the W-dependence of our inputs. If we take the point of view 
that any steep behaviour in the parton distributions only comes about due to perturbative effects, 
then we assume that /^{N) and go{N) are both soft, i.e. either flat or even valence-like when 
the transform to x-space is performed (or at most going like a finite, small power of ln(l/x)). 
This requires that they both be analytic for > 0. We also note that the soft parts of the input 
are common for the whole of the gluon input, as shown in (4.54), i.e. go{N,QQ) is just the soft 
{go{N) + |/(f (IV)) multiplied by a series of the form Ylm=i '^miotsiQo)/^)^■ Thus, as in the loop 
expansion, we may think of gQ{N), and /^{N) as fundamentally soft, nonperturbative parts of the 
input. The parts multiplying these are then really determined by perturbation theory. Since we 
are meant to be expanding our solution for the structure functions in powers of both ag and N, it 
might be argued that we should expand /^{N) and go{N) in powers of N. We feel this is not really 
appropriate since it is the perturbative part of the solution for which we are able to solve, and thus 

Of course, we should be trying to find fundamentally nonperturbative inputs for the structure functions 
rather than the partons, as discussed for the loop expansion. However, in this expansion scheme Jq {N) 
and go{N) are trivially related to F 2 ^o{N) and Fl,o{N) as we will see in the next subsection. 
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which we are able to order correctly, and in the expressions for the structure functions the whole 
of the nonperturbative inputs should multiply the well-ordered perturbative parts of the solution. 


Before using our solutions for the parton densities to construct expressions for the structure 
functions we will solve higher-order renormalization group equations in order to determine the 
general form of the solutions for the parton distributions. This is made easier by the fact that 
^^j:{as/N) and Jfg{as/N) are both zero, and hence, 

= -il}f{»s/N)f,^{N,Q^)+^}g{as/N)go{N,Q^)). (4.57) 

All of the above quantities on the right-hand side are in principle already known. Thus, we have 
an expression for {d fi {N, Q^)/das{Q^)) which can be written in the simple form 


2,r,2,dffiN,Q^] 


aiiQ^ 


dasiQ^ 


= - as{Q^n}g{as/N){go{N) + lfi{N) + ~go{N,Ql))x 


exp 




UasiQ^) 




das{q^ 




9 67r6n 


(4.58) 


where (4.41) has been used. Thus, a^{Q^)d ff {N,Q^)/das{Q^) is simply a sum of two power 
series of the form J2m=o'^miois/d^)^ multiplying input densities and our eigenvalue determined 
evolution factors (the second series rather trivially having Um = 0, m > 0). Thus, all terms in the 
expression are of the same order in this expansion scheme. 

Integrating (4.58) we get 




)/o4a') + /f(V<3S) + (9o(iv) + lfi{N) + <MN,Ql))y- 


9 67r6o \as{Q^) 




Ois{q^ 


exp 


(Qg)^o^(Q,^(r2)/jY) 


IJ Ois{q^) 


a2(r2) 


das{r‘^) 


das{q^ 


(4.59) 


We may express the last term differently by integrating by parts. This gives 




j}g{a,{Q^)/N) 


exp 


(Qo)^^g{as{q‘^)/N) 


L4a,(Q2) a^(g2) 

7},(a.(g^)/iV) 


das{q^ 


-as(Qo) 


j0g(as(Qi)/N} 


d ( 2. rfg{as{r)/^) \ 
da,iq^) \ ^jUas(g^)/N)J 


7}gMg^W)\ 




'as(Q^) } \ \gg 

(as{r‘^)lN) 


exp 

■ In 


Uo,s{q^) oi‘l{r‘^) 

4 21V/ 


das{r‘^) 


da,{q^){go{N) + lf^{N) + ~goiN, Q^)) 


9 67r6o V«s(g^) 


/o"(lV) + /f(lV,g2). 


(4.60) 
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^s{l^fg{o-s/N)/^^g{as/N)) is a power series of the form X]m=o inte¬ 

gration by parts may be repeated indefinitely, producing power series which behave like 
Nag{Q'^)J2m=-k^rn{o:s{Q‘^)/N)‘^, with k increasing by one at each integration. Hence, the solu¬ 
tion to the next-to-leading-order renormalization group equation for is not of a single 

consistent order in this expansion scheme. We may write it as 


/f( 7 V,Q 2 ) = (^^(Ar) + lfi{N) + ~go{N,Ql))x 


as{Q 


,yfg{a,{Q^)/N) 


exp 


(Qi^^g{asiq^)/N) 


’as(Q^) 


a2(g2 


das{q^ 


{Qo 


asiQl)) 


^UN,a,{Ql)) 


as(Q^' 


fQ {N) + higher order in Ug and/or N. 


(4.61) 

Let us examine the form of our solution for ff {N,Q‘^). No part of it is of the same 
form as the zeroth-order parton distributions, i.e. either just a constant soft distribution or 
such a constant soft distribution multiplied by powers of as{QQ)/N and the evolution term 


exp 




Thus, the lowest-order corrections to this zeroth-order solu- 


tion are given by the lowest-order parts of (4.61). 

In the same manner as for the zeroth-order gluon distribution we can determine the general 
form of our input, in this case ff {N,Qq), by considering a change of starting scale. Under the 
change of scale leading to the change of agiQo) in (4.7) each of the as((5o)“dependent terms in 
(4.61) will change. The change of go{N,QQ) has already been chosen to cancel (at leading order) 


the change due to exp 


£ 


0‘s(Qo}7gg(°‘s(<l^)/J^) 


da.siq'^ 


and hence the first term in (4.61) is stable to 


Jas(Q^) ofliW 

the change in starting scale. The last explicit term in (4.61) is also stable since the change induced 
in the ln(as((5Q)) term is of higher order in aglQ^)- The second term does however vary; the 
change due to the variation of as(Qo)jjg(as(Qo)/N)/^gg(as(Qo)/N) is of higher order in agiQ^), 
but goi^jQ^) changes as prescribed in (4.55), and hence the second term changes by an amount 
of the same order as itself, i.e. by a series of the type asiQo) Yin '^n{c(s{Qo)/d^)^■ In order to keep 
the whole of (4.61) unchanged at this order, /f(N, Qg) must change under a change in Qg in a 
manner which compensates this change in the second term. Hence, we choose 


ff{N,Ql) = {go{N) + IuHn))N £ j . (4,62) 

The change of ffmiQo) under a change in Qg may easily be calculated using the known change in 
9g{N,QI), i-e. 


Af^iN,Ql) = 5agiQl)j}giagiQl)/N){goiN) + |/g^(iV) + 5 o(iV,Q^)) + higher order, (4.63) 
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in the limit of small 5. Rather obviously, the simplest choice which satisfies the requirement is 


= (5o(iV) + lfi{N))as{Ql 


'70 (a.(Q§)/iV) 


1 + ^ 90,r 


O^siQl) 


ew[HQo/^gghggMQo)/N)], 


= (<?o(iV) + lfoiN)+UN,Ql))as{Ql) 


'y}gias{Ql)/N) 

i%MQl)/N)- 


(4.64) 

However, as when deriving the O^as^Qo)) input for the nonsinglet structure function within the 
loop expansion (4.22), the invariance of the parton distribution under changes in Qq is not the sole 
consideration. We also require that at lowest order the expression for the input for the structure 
function is renormalization scheme independent, and that at all orders it is factorization scheme 
independent. Also, we think of the scale Agg (or A^s etc.) as the value of Qq where the input for 
the structure function becomes just the nonperturbative input, so /f (A", Qo) should be consistent 
with this. If we add a series of the form 


fi{N,Ql) = Nj2R 


OisiQo) 


{go{N) + luHN)) 


(4.65) 


to (4.64) then the resulting expression still satisfies (4.63), but gives us the flexibility to satisfy the 
other requirements above. Thus, similarly to (4.24) we write 


fi{N,Ql) = ( 50 (A) + lf^{N)+go{N,Ql))asiQl 


jfgiasiQD/N) 

'j^ggMQl)/N) 


+ fi{N,Ql) 


(4.66) 


as the appropriate input, where /;f(A, Qq) is potentially renormalization and factorization scheme 
dependent, and is not yet determined. 

Therefore, the leading correction to the zeroth-order input /^(A) is a series of the 
form A^“^q(q;s(Qo)/A)™' multiplying ( 50 (A) + |/|f(A)). The lowest-order parts of the 
Q^-dependent expression for f^{N,Q'^) are a series of the form N /N)^ mul¬ 
tiplying exp *~Q 2 (g 2 das(q^) and a zeroth-order input, and also a logarithm of 

{ois{QQ)/as{Q‘^)) multiplying the other common factor of unity and a zeroth-order input. 
ff{N,Q^) then contains other terms which are subleading, i.e. of higher order in Ug and/or 
A, to these. 

Hence, this solution for the next-to-leading evolution equation for the singlet quark distribu¬ 
tion clearly demonstrates that we will not obtain well-ordered, factorization-scheme-independent 
expressions in the leading-ln(l/x) expansion in as straightforward a manner as when using the loop 
expansion. 
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In order to investigate this further, we now consider the NLO renormalization group equation 
for gri(A^, This is a little more complicated that that for ff{N,Q‘^): 




2^ dgi{N, Q'^] 

das{Q^) 


= - {o:s^lf{as/N)f^{N,Q‘^) + as^lg{as/N)go{N,Q‘^) 

+ 7°^(a,/iV)/f (iV, Q2) + f^^{as/N)g,{N, Q^)). 


(4.67) 


This may be solved, and integrating some of the resulting terms by parts, we obtain 


gi{N,Q^) = 




as(Q^) ois{q^) 

{9o{N) + IfiiN) + ~go{N, QD) + g,{N, Ql 


das{q )x 

as(Qg)^ 0 ^(Q,^(^ 2 )/^) 


exp 




OisiQl] 


L 4 a,(Q 2 ) 02(^2 

+ higher order in and/or N. 


das{q^ 


^9 evrfto 

(4.68) 

The integral in the first term gives a series of the form am(«s/-^)™' plus a term of 

the sort ln(as(Qo)/®s(Q^))- those terms explicitly shown above are the lowest-order terms 
in gi{N,Q‘^). The input 5 'i(A^, Qq) is present in order to compensate for the change in 

das{q^) under a change in Qo- However, this means it is a 
series of the form ^^(Qg) an(as(Qo)/-^)’^’ such is higher order than the other terms 

present, and may be dropped. We see that the solution for gi{N,Q‘^) contains terms of the same 
apparent order as the zeroth-order solutions for the parton distributions. Therefore gQ{N,Q‘^) 
seems to be only part of the lowest-order expression for the gluon distribution, and the terms 
shown above appear to be on the same footing. These are of a similar form to the leading parts of 
ff{N,Q^), with one term being a logarithm of {as{Qo)/as{Q^)) multiplying /(f(iV,Qo); but the 


series multiplying exp 


l»siQ^) a2(g2) 

than the corresponding series in the case of fi{N,Q^). 


ras(Qo) %a(as(g )/^) j ( 2\ 


is a power of N, or equivalently of a^, lower 


So we have the rather undesirable situation that solutions to the NLO renormalization group 
equations for both the quark and gluon distribution functions result in solutions which are of mixed 
order, but fortunately, ones from which we can clearly extract the leading-ln(l/x) behaviour. 
However, the situation gets worse, as may be seen by solving the NNLO renormalization group 
equations. For the quark distribution function this is 


as(Q^ 


dfi{N,Q^) 

das{Q‘^) 


= -{asiQ^)^ffias/N)f^{N,Q^) + asiQ^)^fgias/N)go{N,Q^) 

+ j}f{as/N)f^{N, Q2) + jj^{as/N)g, {N, Q^)) . 

(4.69) 
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As with the equation for {N, ) all terms on the right are already known in principle (it is clear 

that this will be true for the equations for all the f^{N,Q‘^)). Moreover, part of the right-hand 
side of (4.69) is of the same order as the right-hand side of (4.57), as is seen by inserting the leading 
parts of fi{N,Q‘^) and gi{N,Q‘^), i.e. 


2^r.2,dfiiN,Q^) 


OisiQ 


dasiQ'^ 


= -as{Q‘^)j}g{as/N) 


d^°KllgMq‘^)/N) + ^7}g{as{q^)/N)) ^ 


^2 


(goiN) + |/o (^)) + 9o{d^,Qo) ) exp 


as{q^) 

f<^^^Ql)^gg{as{q^)/N) 

Ua,(Q2) 


das{q )x 


a 2 (g, 2 ^ 


das{q^ 


-f^{N)as{Q^ 


9 Gvrfeo 


In 




+ higher order in and/or N. 


(4.70) 

Solving (4.69) we find that ^‘jj:{as/N) and ^‘^^{as/N) play no role as far as the leading part of the 
solution is concerned, and 

l%{as{Q^)lN) VJa4Q2) «s(r) 


{go{N) + |/o^ {N)) + go{N, Ql) ) exp 


as(Qg);yO (q,^(^ 2 )/^) 


LJ ocs{Q'^) 


al{q^^ 


daa{q^ 


+ 


l(^2Nf J asiQl) 
2V9 67r6o 


fQ{N) + higher order in and/or N. 


(4.71) 

Thus, the leading part of /^(A", Q^) has terms of the same order as the leading part of /f (A, Q"^), 
i.e. a sum of the form A J2m=o ^■m{oia/N)'^ multiplying exp d as{q^) and soft 

inputs. We also explicitly include above the parts of /^(A, Q^) which have terms of the form 
ln{as{Qo)/o(s{Q^)) directly multiplying /^{N), even though these are not of the same form as 
in ff{N,Q^). We will take these terms to be part of the lowest-order parton distribution since 
they are not a power of or A higher than any terms already produced. However, we ignore 
terms of the form ln(Q!s((5o)/«s(Q^)) multiplying series of the form Aaj„(as/A)™', such as 
/f(A, Qq), because they are clearly of higher order in A than the term consisting of logarithms 
multiplying /,;/(A). 

Similarly, for the gluon the presence of lowest-order terms is not limited to gi{N,Q‘^] 
also seen in g 2 {N,Q‘^): 


but is 


<72(A,Q2) = 




(Qo)(7lg(a,(g2)/A) + 2 

/ 2 ^ j 

(^s{q^) 


(50(A) + |/o^(A) + 5 o(A, Qo)) X exp 


^^o)^^gg{as{q^)/N) 


U oisiQ^) 


a 2 (g 2 ^ 


daa{q^ 


(4.72) 




+ i/o(^)o n 


2V9 67r6o \as{Q‘^) 


+ higher order in and/or A. 
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This phenomenon persists to all orders, but is independent of the ^""’s for n > 1. It can be proved 
by induction in a straightforward manner that the leading part of the solutions to the nth“Order 
renormalization group equations are 


f^{N,Q^)=as{Q^) 




7 °,(a.(gg/iV) 

1 


exp 




das{q^) 


as(Q") 


(n- 1)! 

{goiN) + |/o^(A^) + go{N, Ql)) exp 


Ois{q^) 

(QoM (a^(g 2 )/^) 


das{q^) X 


n — 1 


iJas(Q^) 


al{q^) 


das{q^) 


+ 


for n > 2, and 




(4.73) 


gn{N,Q^) = 


1 


d^°KllgMq^)/N) + ^j}g{as{q^)/N)) 


n'- 


asiq'^) 


das{q^) X 


{go{N) + |/o*(iV) + ffo(A^,go))exp 


das{q^) 


^s(Qg) ^ 0 ^(o;^(g 2 )/jy) 

UasiQ^) oil{q^) 

(4.74) 

for n > 0. One can therefore find well-ordered parts of the full solutions for f^{N, Q^) and g{N, Q^) 
by summing all such terms, leading to 

pS'^A7Z^2^__ , 4fS/AT\ , ;; I AT 


f^iN,Q^)=asm 


^%iasm/N) 

exp 


igoiN) + lf^{N) + go{N,Qi))x 


^^°^ggias{q^)/N) jlg{as{q^)/N) + ^j}g{as{q^)/N) 


+ foiN) 


Uo^siQ^) 

my 

as{Q^)J 


a 2 (g 2 ) 


+ 


as(g^) 


-das{q) 


a 


4 2 iVj 
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+ fI {N,0^) + higher order in and/or N, 


(4.75) 


and, 


g{N, g2) =( 50 (AT) + + goiN, g^))x 

exp 


f»s{Qo)x^o^(a,{q^)/N) ^ y{asiq‘^)/N) + y}g{a,{q^)/N) ^ ^ 

2 T 2\ ' T 2^ dCXs[Q ) 

UasiQ^) ) «^(r) 


(4.76) 


-IfoiN) 


yoy] 

as(g^)y 


4 2 IV^ 

9 67rbo 


+ higher order in Ug and/or N. 
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We may also write the expression for the rate of change of the quark distribution 


daa{Q‘^ 


=aa{Q^)l)g{aa{Q^)/N){g^{N) + {N) + ,Ql))y< 


exp 


adQo)^o^(a^(g 2 )/jv) ^ ^lg{aa{q^)/N) + ^^)g{as{q^)lN) ^ 2 ^ 


Uas(Q") 


-laamfoiN)- 


Gnbo Vas(Q^ 


4 


2Nffas{Ql)\ 


+ higher order in and/or N. 

(4.77) 

This last expression will be important since at leading order a‘l{Q‘^){dF 2 {N, Q'^)/das{Q'^)) is di¬ 
rectly related to dF 2 {N,Q‘^)/d InQ^ which we will wish to study as well as F 2 {N,Q‘^). 

4.3. The Leading— ln(l/x) Expansion: Strncture Pnnctions. 

It is now possible to examine the form of the solutions for the structure functions. We can 
construct the leading part of the full solutions by combining our leading solutions for the parton 
distributions with the zeroth- and first-order coefficient functions. This gives 

FLiN,Q^) = aaiQ^)Cl,{aa{Q^)/N)igo{N) + I fi{N) + , QD) x 

<'^‘^Q°'^%g{as{q^)lN) ^ llg{ois{q^)/N) + ^x^)g{as{q^)/N) ^ 

Uc^siQ-^) 


exp 


a 2 (g 2 ) 

2\(r^S 4 r^g \ fS 


+ 


+ CamiCUo-ICl,,o)foiN) 


^s{Qo, 

a.s{Q^] 


as{q^ 

4 2 iVj 
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-das{q^ 


-|- higher order in ag and/or N, 

(4.78) 


where (4.49) has been used. Also 


F 2 {N,Q^) = +Cl,{ag{Q^)/N)){g,{N) + {N) + ,Ql)y. 


exp 


dQo)x^o^(ag{q^)/N) jlg{ag{qy/N) + ^jjg{ag{q‘^)/N) 


U ctsiQ^) 

+ foHN) 

We can also write 


a 2 (g 2 ', 


-k 


ag{q^ 


-das{q^ 


(4.79) 


as{Ql] 

as{Qy 


4 2 iVj 
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+ /f (iV, Qo) + higher order in ag and/or N. 


=(a.(Q=)7;,(a.(Q=)/JV)+ 


dag{Q'^ 


ag{Q^)Ciyag{Q‘^)/N)fgg{ag{Q‘^)/N)){g^{N) + {N) + UN,Ql))x 

'■“”(‘^“)7gg(as(92)/iV) ^ ylg{ag{qy/N) + ^x^)g{ag{q^)/N) 

L4a3(Q2) 


exp 


al{q^^ 


+ 


as{q^ 


das{q^ 




^sjQl 

OisiQ^ 


.4 

9 67rbQ 


-k higher order in ag and/or N. 


(4.80) 
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In each case, since we only consider the singlet structure functions in this expansion scheme, we 
drop the superscript S for the structure functions for the rest of this subsection. 

It is clear that each of these expressions will be factorization scheme independent since 
each represents the expression for a physical quantity up to corrections of a form different from 
the terms explicitly appearing and which we have deemed to be higher order in our expansion 
scheme. {al{Q‘^){dF 2 {N,Q‘^)/das{Q‘^)) has as much right to be considered as a physical quan¬ 
tity as F 2 {N,Q‘^), since at leading order it is proportional to {dF 2 {N,Q‘^)/dlnQ‘^). However, we 
would like to find out if each of these expressions may itself be more rigorously ordered and, as 
a by-product, split into factorization-scheme invariant pieces. This would be desirable from a 
theoretical point of view, since the current expressions will soon be seen to be incompatible with 
renormalization scheme consistency. Also, from the practical point of view we do not know all 
of the terms appearing in the explicit parts of these expressions in any factorization scheme, i.e. 
^gg{as{Q‘^)/N) is not at present known, and we would clearly like to have some sort of scheme- 
independent expression for the structure functions involving terms we already know. 

In order to see if we can obtain scheme-independent and thus physically relevant subsets of 
the solutions (4.78)-(4.80), we shall examine in detail the form of these solutions. Since it is the 
least complicated, and because we will find it useful, we begin with the expression (4.78) for the 
longitudinal structure function. Since this expression contains two terms of rather different form, 
i.e. the first depending on the factor exp ^ ) and the second having an 

evolution unenhanced by leading-ln(l/x) terms, it would be surprising if these terms were not 
separately factorization scheme invariant. It is simple to check that this is indeed the case. 

Beginning with the second term, it is clear that this is factorization scheme independent. 
The power of (as(Qo)/<as(Q^)) comes from the one-loop contribution to jjg{as/N) which is both 
factorization scheme and renormalization scheme invariant and which we have already written in 
its actual numerical form. as(Q^) and /^{N) are both clearly renormalization and factorization 
scheme dependent, and the one-loop contributions to the longitudinal coefficient functions are 
also, i.e. g = ( 2 / 37 r) and C'£ g = {2Nf/QTr) in all schemes. Hence the whole term is both 
factorization scheme and renormalization scheme independent, as it must be, depending only on 
leading-order quantities. 

The first term is not as simple to deal with. It must be factorization scheme independent; 
however, many of the pieces appearing in the term are clearly not, e.g. jjg{as/N) and ^{as/N). 
So, to a certain extent it must be the interplay between the 7 ’s, the coefficient function and 
the parton inputs which leads to a factorization-scheme-independent result. Remembering our 
discussion of the solutions for the structure function in the loop expansion, we see that there is at 
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least one way of writing this term as the product of two factorization-scheme-independent pieces, 


as{Q^)Cl^{as{Q‘^)/N){g^{N) + |/o^(iV) + Q^)) x 


^^^°hlg{as{q^)/N) ^ %g{a8{q^)/N) + ^^)g{as{q^)/N) ^ ^ ^^ 2 - 

t„‘2\ do:s[q , 


'as{Q^) 




,(Qg )^0 (Q,^(^ 2 )/jy) 




daa{q^) exp —in 


'«.(Q 2 ) 


(Qg) 7^g(«^(g^)/A^) + ^l)g{as{q^)/N) ^ d 
02 ) as{q^) das{q^) 


Cl,{aa{q^)/N) 


O^siQl) 

as(Q^) 


daa{q^ 


(4.81) 

The left-hand side is factored into its value at Qq a term which gives the evolution from this 
value. Clearly these are physically distinguishable and as such must be individually factorization 
scheme invariant. Also, we already know that exp "(q°) j^ ^ ^factorization 
scheme invariant on its own, and exp[— hi{aa{QQ)/oigiQ'^))] clearly is. Hence, we make the dehnition 


(4.82) 

where ^t{Q'^,Qo) must now be factorization-scheme independent (as can be checked using the 
rules (2.10)-(2.12)) and {Qq,Qq) = 1. 

Having isolated the factorization-scheme-independent parts, we may factorize (4.78) com¬ 
pletely into factorization-scheme-invariant input and evolution parts, i.e. 


dQo)7i^(a^(g2)/iv) -h ^7}giaaiq^)/N) d 


/qs(Q^) 


das{q^ 


CiAaa{q^)/N) 


9 J^T 

Fl{N,Q^) = as{Ql)-^{go{N) + |/o^(iV))exp 


dQol^O (q,^(^ 2 )/jy) 


'a4QT 


das{q^) X 


Cl,{aa{Ql)/N) 

{^9 


5o(a^, Ql) 


9o{N) + IfiiN) 


exp k>+(Q2,Qg) - In 


o^siQl) 

as{Q‘^) 


-f a^(Q^)|^ "" 27n ’^ + higher order in «« and/or N. 

° (4.83) 

It is now possible to attach direct physical significance to each of the factorization-scheme- 
independent pieces appearing in this expression. 

We hrst consider the inputs. Going to (4.79) for the moment we see that /^{N) is the only 
term in F 2 {N,Qq) which is zeroth-order in our expansion scheme. As such it is the zeroth order 
input for F 2 {N,Q‘^), and we may write 




f^{N) = F2 ,o{N), 


(4.84) 
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and this is one of our two fundamentally nonperturbative inputs. Also, from (4.83) we see that 
the total as((5o)"iiidependent input for Fl{N,Qq) (once we have divided out a single power of 
«s(<5o)/(2vr)) is 

+ ifiiN)) + (;hzli^)/S(iV), (4,85) 

This whole expression may therefore be written as Fl^q{N), and we have 


OAT 9 

FlAN) = -^i9o{N) + 


(4.86) 


Fl, oiFl) is our other fundamentally nonperturbative physical input. We note that the expression for 
F 2 fi{N) in terms of the parton inputs is the same in this case, i.e. (4.84), as in the loop expansion 
(4.36). However, the definition of Fl^q{N) in terms of the parton inputs, (4.86), is not the same as 
in the loop expansion (4.37). Thus, the definition of the nonperturbative gluon, gQ{N), is different 
in the two expansion schemes. 

We also make a similar definition for the part of the input for Fl{N, Q^) which is of the form 
(^niQl)iasiQo)/Ff)"- multiplying {go{N) + ^fo (N)), or more correctly multiplying {Fl,o{N)- 
((36 — 8Nf)/27)F2^o{N)), to complete our definition of the input in (4.83) and write 


2Nf 


{go{N) + lf^{N) + go{N,Qi)) 


CiJas{Ql)/N) 


/^g 


= Fl,o{N)- 


36 - 8Nf 
27 


(4.87) 


F2,o{N)+Fl,^{N,QI) 


This whole expression must be both factorization scheme and renormalization scheme independent, 
facts which reveal information about the form of the gluon input. C'£ ^(A^,as(Qo)) both renor¬ 
malization and factorization-scheme dependent, so the scheme dependence of go{^,Qo) must be 
precisely so as to cancel this out, i.e. 


{9o{N) + lfi{N) + ~goiN,Ql)) 


^(5o(iV) + |/o^(iV))f 1+E 

^ m=l 





eMHQo/^LLhgg{as{Ql)/N)] 


= (MN) + |/o^(iV))go(iV,g2)(c£^^ o/C£^i(a,(Q2)/iV))exp[ln(g2/yi^^)^o^(«^(g2)/^)]^ 

(4.88) 

where go(A',g^) is scheme independent. Hence, there is no reason for Fi^o(«s(go)/-^) to depend 
on the leading-order longitudinal gluon coefficient function at all, and indeed, a natural choice 
seems to be that (1 -|- Ylm=i 9 o,m{c(siQo)/Ff)'^) is chosen equal to (C'£ ^ i(-^) CKs(f5o))) (it is 

difficult to see what else it could be chosen equal to), i.e. go(A^, g^) = 1) and therefore that 


FlAN) - 


36 - 8Nf 


27 

= (fl,o{N) - 


F2An) + FlA^,qI 

36 - 8Nf 


27 


(4.89) 


F2AN) ) eMHQAALL)Agg{as{Qt)/N)]. 
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Hence, we have a prediction for the input for the longitudinal structure function at small x in terms 
of the nonperturbative inputs and some scale All (where All = Agg from the previous subsection). 
As with Ajls earlier, All is the scale at which the input is equal to the nonperturbative input alone, 
and hence we would expect it to be typical of the scale where perturbation theory starts to break 
down. As we have already stressed, Qq is a completely free parameter. We have constructed the 
solution to be insensitive to Qq at leading order, but there is clearly some residual Qg^dependence. 
Hence, there will also be some optimum Qq to choose as the starting scale. 

We may now examine the terms governing the evolution, continuing the convention started 


in (4.82). exp 




is the factorization-scheme-independent factor coming 
from the eigenvalue of the zeroth-order anomalous dimension which will govern the small-x growth 
with Q^. As such we make the definition 




(Ql) 


>as{Q^) 


a 2 (g 2 


das{q^ 


(4.90) 


Of course, the other evolution factor resulting from the eigenvalues of the zeroth-order anomalous 
dimension was simply unity and as such ^q{Q^,Qq) does exist, but is implicitly zero. However, 
there is a correction to this factor of unity in (4.78)-(4.80), and we make the definition 




9 Ovrftn 


OisiQ^ 


(4.91) 


Having made these factorization-scheme-invariant definitions for the inputs and evolution we 
may write the solution for Fl{N,Q‘^) as 


Fl{N,Q^) = 


“s(Qo)/p_ /An , AT /^36 SNf 


27r 


-{FL,o{N)FFL,o{N,Qt)- 


27 


exp 


Ql) + Ql) - In' 


as{Q^ 


i^ 2 , 0 (iV))x 


+ ^4^f 36 )F 2 ,o(iV)exp 


27r V 27 
+ higher order in ag and/or N. 




as{Q‘^ 


(4.92) 


It is now relatively obvious how we may separate out the “leading part” from this expression. 
[<!>/"(Q^, Qq) — ln(as((5g)/as((5^))] contains the same type of terms as ^q{Q‘^,Qq), but each is 
a power of N higher. Thus Qg) “ Ai{as{Q‘^)/as{Q‘^))] is subleading to <I>(/((5^, Qg), and 

indeed 4>/"((5^, Qg) is a renormalization-scheme-dependent quantity, as it must be in order to 
absorb the change in resulting from a change in the definition of as{Q‘^) under a 

change in renormalization scheme when working beyond leading order. Hence, we should factor 
,Q‘l)—hi{as{Q‘l) / aaiQ"^))] out of the first term. Since [4>^((5^, (5g)-ln(Q;s((5g)/as(Q^))] 
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is of the same form as the zeroth-order-in-iV part of then it should also be factored 

out of the leading-order expression for Fl{N,Q‘^). Thus, we are left with 


Fl{N,Q^) =^f^[FL,o{N)+~FL,o{N,Ql) - 


+ 


27r 

as{Qo) f 36- 8Nf 


36- 


2^)F2,o(iV))exp[$+(Q2,g2)] 


27r 


277r 


(4.93) 


F2,o(iV). 


Using the one-loop running coupling constant, as is appropriate for a leading-order expression, the 
whole of (4.93) is not only manifestly factorization scheme independent, but also renormalization 
scheme independent, as one would hope. Also, FQo) is constructed precisely so as to make 
the expression unchanged, at this order, under a change in starting scale. Hence, within this 
expansion scheme (4.93) is genuinely the leading-order expression for Fl{N,Q‘^). 


We now turn our attention to the more phenomenologically important case of the struc¬ 
ture function F 2 {N,Q‘^). Here we immediately have an ambiguity: should we consider the ex¬ 
pression for the structure function itself in this expansion scheme, or that for its CKs-derivative, 
a‘^{Q‘^){dF 2 {N, Q'^)/das{Q'^)) (which, using the definition of the running coupling (2.5), is directly 
related to {dF 2 {N, Q‘^)/dlnQ‘^) at leading order)? It may be argued that in certain senses the latter 
is more natural because it is a “real” perturbative quantity, beginning at first order in as does 
Fl{N,Q^). 

There is in fact a distinction between F 2 {N, Q'^) and {dF 2 /d InQ^) in the usual loop expansion. 
Differentiating a fixed-order expression for F 2 {N,Q‘^) and using the /3-function evaluated to the 
appropriate order in Ug results in the fixed-order expression for {dF 2 /d InQ^)^® plus terms of 
higher order in as{Q‘^), which depend on the /3-function beyond lowest order (and thus are absent 
when working to leading order only). Hence, the size of these extra terms is of the same order 
as the renormalization scheme uncertainty, and therefore the distinction between the fixed-order 
expressions for F 2 {N,Q‘^) and {dF 2 /d InQ^) is of similar magnitude to the distinction between 
renormalization schemes. 

There is also a distinction when using the small-x expansion, and it appears more graph¬ 
ically, and is not only dependent on terms in the /3-function beyond lowest order. We distin¬ 
guished between F 2 {N,Q‘^) and al{Q^){dF 2 {N,Q'^)/das{Q^)) in (4.79) and (4.80) because the 
expression to a given order in Ug and N even for al{Q‘^){dF 2 {N,Q‘^)/dag{Q‘^)) is no longer sim¬ 
ply obtained by differentiating F 2 {N,Q'^) at given order, and conversely given order in F 2 {N,Q‘^) 
is not obtained just by integrating al{Q‘^){dF 2 {N,Q‘^)/das{Q‘^)) with respect to Ug at fixed or¬ 
der. This was clearly illustrated when finding the leading part of fi{N,Q‘^) by solving the equa¬ 
tion for a‘^{Q'^){d ff {N,Q'^)/dag{Q‘^)) earlier, and it leads to a distinction between even the LO 


(dF 2 {N,Q^)/d InQ^) is analogous to the longitudinal structure function, i.e. it has an overall power of 
cts{Q^) and the lowest-order expression consists of first order-in-Q;s((5Q) inputs multiplying the eigenvalue- 
determined evolution terms )^+.- 
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{dF 2 /d InQ^) and the derivative of the LO F 2 {N,Q‘^) in this expansion scheme. Thus, even at 
leading order we have to decide which of the two expressions to use, although we would hope that 
there is not too much difference between the choices in practice. We will postpone the decision for 
the present and will derive the form of the leading-order expressions for both F 2 {N,Q‘^) and its 
derivative. We can then examine these expressions and use them in order to help us decide. 

For technical simplicity we begin with the expression for a‘l{Q‘^){d F 2 {N, Q‘^)/das{Q‘^))- Using 
the definitions introduced in our discussion for the longitudinal structure function we may write 
(4.80) as 




da.s{Q^ 


CUo 


iFL,o{N) + FLAN,Ql)- 


2Nf\Cl,{a,m/N) 

+ ^ti 

+ higher order in and/or N. 


36 - 8Nf 
27 


i^2,o(iV))x 


exp[4>+ (Q2, Ql) + $+ {Ql, Q2)] - asiQ‘^)^^^F 2 ,o{N) exp[4>i (Q^ Ql)] 


(4.94) 

The factorization scheme independence of this complete expression guarantees that the term 
il}g{as{Q‘^)/N) + C'2,i(«s(Q^)/^)7gg(«s(Q^)/^))(C'£,i,o/C'£,i(«s(Q^)/^)) u factonzatiou- 
scheme-invariant quantity, of the form X]m=ostraightforward to 
verify this. Indeed, it was shown by Catani and Hautmann that {'yjg{as{Q'^)/N) + 
Cfq {asiQ"^)/N)^gg{as{Q‘^)/N)) and C'£ i{as{Q‘^)/N) could each always be expressed in terms of 
the product of a factorization-scheme- and renormalization-scheme-independent factor (which 
they calculated) and a scheme-dependent factor, where the scheme-dependent part was the same 
for both [18]. Using their results for the scheme-independent parts it is a trivial matter to find 


that 


o / (J9 

= + E(4<“-(‘3")/"))"). 

n=0 ^ 

(4.95) 

which is clearly both factorization scheme and renormalization scheme independent. Thus, making 
the definition 


asiQ^H^iasm/N) = 


27Tbn 


oo 

- + 

n =0 


99 


>sm/N)r 


(4.96) 
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our expression for the explicit part of (4.80) is entirely in terms of factorization-scheme-invariant, 
and hence physically meaningful quantities, i.e. 




dF2iN,Q^ 

da.siQ^) 


=as{Q^)lldc^s{Q^)/N){FL,o{N)+FLfl{N,Ql) - 


2, /36-8iV/ 


27 


F2,o{N))x 


ex' 


:p[4>+(g2, Ql) + 4>+(g2, g2)] _ 4 exp[4>^(g^ Qo)] 

y OTTOo 

+ higher order in and/or N. 

(4.97) 

We no'w have an expression ■which is analogous to that for Fl{N,Q‘^) in (4.92), except that it 
has still not been explicitly separated into inputs and evolution terms. In order to do this we must 
rewrite (4.97) as 

=«*(Qo)72'l(««(Qo)/^^)(^l,o(A^)+ ^L,o(A^,g^) - 


da.siQ'^ 


4 2iV/ .^2 


27 


ex' 


p[ch+(g2, Ql) + Q^) + ^i{Ql Q^)\ 


9 Gvrfer 


aa{QQ)F 2 ,o{N) exp 


(g^go) - In 


as{Ql) 

as(g^) 


+ higher order in as and/or N, 

(4.98) 


where <I>|/(g^, g^) is a series of the same form as <I>|/(g^, g^), defined as 

fO‘s{Qo) fj 

4>^(g^,go) = - / — T^Mo‘s{q^)7lLiFl,as{q^)))das{q‘^). (4.99) 

JcsiQ^) das{q^) 


It is now clear how we obtain the “leading part” of this expression. All of the inputs are 
of leading order, but we must factor out the subleading parts of the evolution. [4>/“(g^,gg) + 
^^(g^,gg)] is subleading to (g^,gQ) (we note that ^t{Q^,QQ) is a renormalization-scheme- 
independent contribution to this subleading evolution), and must be factored out. Since it is of 
the same order, so must [<I>/(g^,gg) — ln(as(gg)/as(g^))]. This leaves us with the leading-order 
expression 

=a,(g^)72^i(as(g^)/IV)x 

(Fz.,o(A^) + ~FL,oiN,Ql) - (^^-^^)F2,o(iV))exp(4>+(g2,g2)) 

(4.100) 

Again, this expression is insensitive to changes in starting scale, up to higher order, and using the 
one-loop coupling constant, is renormalization scheme independent as well as factorization scheme 
independent. 
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Finally we consider the expression for F 2 {N,Q‘^) itself. It is now a relatively simple matter 
to write this in terms of factorization-scheme-independent quantities and in terms of inputs and 
evolution terms. Using the definitions we have already made and defining {Q^, Qo) by 


^^(Q^Qo) = - 

we may write (4.79) as 


/ 

J O' 


(Qo) 


d f as{q^hl^{as{q^)/N)^ 

^ - (TTFt -- ]das{q ). (4.101) 

„^(Q 2 ) da,{q^) V 7 gg(^,as(r)) / 


F2iN,Q^)=asiQt,) 




lUasiQD/N) 


{Fl,o{N) + Fl,o{N,QI)- 


36 - 8Nf 
27 


i^2,o(iV))x 


eXp[4>+(Q2, ^ g2) ^ |,+ (g2^ g2)] ^ (g2^ g2)] (4.102) 

+ fi {N, Qq) + higher order in a* and/or A^. 

So the factorization-scheme-independent input F 2 {N,Qq) is 


i"2,0(iV) +«.(g^)?^|^4§|^(FL,0(iV) + ^L,0(iV,g2) _ 
Iggy^siQoj/N) 


36 - 8 iV/ 
27 


i"2,o(iv)) + /f(iv,g 2 ). 

(4.103) 

By construction it is guaranteed that the change of (4.103) under a change in starting scale 
will cancel the change in the evolution in the first term in (4.102) under a change in start¬ 
ing scale up to higher orders. It is also clear that (4.103) is both renormalization scheme 
and factorization scheme independent, as we require, as long as /f^(A^, Qo) is scheme indepen¬ 
dent. In fact, the requirement that if Qg ~ ^ll the input reduces to the nonperturbative 
input F 2 fi{N) determines /f(A^, Qq) uniquely. It must be the scheme-independent quantity 

3 6-8Aff \ 


-0's{Ql){'ylLio's{Ql)/N)/^^g{as{Ql)/N)){FL,o{N) - jF 2 ,o{N)). Hence, the input for 

F 2 {N,Q^) is 


F2{N,Qi) = F2 ,o{N) + aM 


2\12l(^s{Q‘q)/N) ^ t^T r^2 


lgg{c^s{Ql)/N) 


Fl,o{FI,Qq) + higher order in and/or N, 


(4.104) 


and the expression for F 2 {N,Q‘^) is 


F2{N,Q^) =aM) 


2\ 72l(-^) Ots{Q‘l)) 


^UN,as{Ql)) 


{FLfl{N) + FL^^{N,Ql)- 


36 - 




27 7 


F2,o{N))x 


exp[4>+ (Q^ Ql) + 4>+ {Ql, + |,+ (g^^ F 2 ,o{N) exp[4>i (g^ gg)] 

'36 - 8Nf' 


- as{Qo 


2a72l(^> as(Qo)) 


7° (iV,a«(g^)) 


{Fl,o{N) - 


27 


(4.105) 


^2,0 (iV)) 


-|- higher order in and/or N. 


By comparing (4.104) with (4.93) we now see that there is a very direct relationship between 
the inputs for our two structure functions at small x, i.e. we have a definite prediction, up to 
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additive nonperturbative parts which are flat at small x, for one in terms of the other, as well as 
approximate predictions for the form of each. Also comparing with the form of (4.100), these two 
inputs for the structure functions are directly related to the slope of {dF 2 /dh\Q‘^) for small x at 
Qq. As already mentioned, we do not yet know at what Qq it is most appropriate to choose the 
inputs, but it is a nontrivial requirement that the inputs for the three expressions are of the correct 
form and related in the above manner at any Qq. 

Our expression for F 2 {N, cannot be split into leading and next-to-leading pieces in quite 
such a clear and symmetric way as our previous two examples, essentially because it begins at 
zeroth order, not at first order in as- Looking at (4.105) it is clear that there is only one LO (in 
our expansion scheme) input multiplied by a LO evolution, and that is F 2 ^o{N) multiplying unity. 
Hence we take the leading-order expression to be given simply by 

hAN,Q^) = F2 ,o{N). (4.106) 


This is obviously completely independent of ag, and is rather trivial. At next-to-leading order, 
or equivalently, at leading-cTs-dependent order, we include the whole of (4.105) except that we 
factor exp[<l>^(g^, Qq) + Qq)], out of the first term, i.e. we have NLO inputs multiplying 

LO evolution, and vice versa, as well as F 2 ,o{F',Q‘^)- Hence, at leading-ag-dependent order we 
have 


F2,o{N,Q^) = F2,o(iV)exp[ci>r(Q",Qo)] 


+ «s(Qo) 
-as(<3o) 


AlMQD/N) 

Ag{as{Ql)/N) 

lAMQD/N) 

AaMQl)/N) 


{FlAN)+~Fl,o{N,QI) - 
(FL,o(iV)-(^^^^)F2,o(A^)) 


F2,o(iV))exp[cl>+(Q2,Q2)] 


+ higher order in and/or N. 

(4.107) 

In this expression there are clearly no terms which mix if we were to make a change in definition 
of the coupling + ea^, and hence we can consider it as a leading-order expression. If the 

one-loop coupling is used, it is both factorization scheme and renormalization scheme independent. 


We now have the full set of LO expressions in the leading-ln(l/x) expansion scheme. We 
could obtain the correct scheme-independent expressions for the structure functions at higher 
orders within this expansion scheme. We choose to finish at leading order, however. The labour 
required to obtain higher-order expressions becomes progressively greater and we would obtain 
expressions requiring unknown anomalous dimensions and coefficient functions. Working to NLO 
we would need to calculate all the NLO evolutions and all the NLO inputs. This would require all 
the NLO expressions for the anomalous dimensions and coefficient functions, i.e. the T^a’s, 7/a'® 
and C'/ 2 ’s. There is optimism that these NLO terms will soon be known [35], and once this is so 
the full NLO scheme-independent expressions should be calculated. 
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We should make some comments about our LO scheme-independent expressions for Fl{N, Q'^), 
F 2 {N,Q‘^) and {dF 2 {N,Q‘^)/dln{Q‘^)). First we note that all depend on factorization-scheme- 
independent combinations of the 7 ga’s, 7 )a’s and C^ys (along with the input parton distributions, 
where g'o(Ai, Qq) is also factorization scheme dependent). There is, however, no terribly simple 
prescription for how one uses the anomalous dimensions and coefficient functions in order to ar¬ 
rive at the expressions. One does certainly not take the known anomalous dimensions, solve the 
full renormalization group equations, and combine with the known coefficient functions. This in¬ 
troduces terms we do not even include in (4.78)-(4.80), let alone in (4.93), (4.100) and (4.107). 
The additional terms will depend on the factorization scheme used, and can be very large. The 
method of determining LO scheme-independent expressions is also more complicated than solving 
the renormalization group equations to a given order and combining with coefficient functions to 
a given combined order, as in the loop expansion. The only way to obtain the correct expressions 
using the parton model directly seems to be to calculate carefully and keep all terms of a given 
type, as explained. 

We should also make some mention of why factorization scheme dependence can be very large 
in this expansion scheme. In order to do this let us consider (4.81) as an example. A representative 
example of the way in which factorization-scheme-dependent calculations are done is to consider 
this expression evaluated with yg{as{Q'^)/N), /N) and ^{as{Q'^)/N) known in some 

particular factorization scheme, but yg{as{Q‘^)/N), which is unknown, either set equal to zero, or 
guessed by imposing some ansatz such as momentum conservation (which may well give completely 
the wrong answer). 

As a first comment we consider the input. The input in terms of parton distributions 
is multiplied by ^{asiQ^/N), a series of the form X]m=o which, in gen¬ 

eral, becomes singular at N = A((5o)- Under a transformation of the type (4.45) this series 
will be multiplied by {Uy{as{QQ)/N))~^ in order to compensate for the change in the input 
parton distributions. A number of scheme transformations that are considered have the series 
Uy{as{Ql)/N) = I + also becoming nonsingular at A^ = A((5o) (e.g. 

[24][25]). This nonsingular behaviour leads to powerlike growth of the form as x ^ 0, 

but the magnitude of the powerlike behaviour and the manner in which in which it is approached 
depends on the strength of the singularity and/or on the precise behaviour of the coefficients in the 
series. Hence, changes in definition of the form in (4.45) with this type of singularity can lead to 
very marked differences in the form of the gluon at small x, or if the gluon is kept roughly constant 
(e.g. it is attempted to predict the form of the input structure function by assuming a form for 
the input gluon), to significant changes in the form of Fl{x,Qq). Very similar considerations also 
hold for F 2 {x,Q^) because the input depends strongly on Cy ^{asiQy/N), which transforms in 
the same way as C'£ ^ ^{asiQy/N) under changes in factorization scheme; i.e. the parton input 
(or prediction for the structure function if the parton input is assumed, e.g. to be flat) depends 
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very strongly on scheme. We also note that some of the scheme changes away from the standard 
MS scheme (e.g. the SDIS scheme) involve Ugg{as{Qo)/N) with larger coefficients in the power 
series than in ^ (oisiQo)/^) defined in the MS scheme. This leads to C'£ ^{as{QQ)/N) defined 
in the new scheme to have negative coefficients and therefore to Fl{x,Qq) developing a powerlike 
behaviour with negative magnitude unless either the input gluon or quark (or both) has a powerlike 
growth of the form itself, with a large enough multiplicative factor, in order to counter 

this effect. This must be borne in mind when using such schemes. 

Examining (4.81) we can also see how the evolution may be strongly factorization scheme 
dependent. As we have already mentioned, the whole of ^i{Q^,Qo) must be used in order to 
obtain a factorization-schemeMndependent expression. If 'ygg{as{Q‘^)/N) is omitted, or the wrong 
'ygg{as{Q‘^)/N) used, then this can give a completely misleading result for the evolution. The 
integrand in ^^{Q^,Qo) is a series which is a power of as{Q‘^) down on the series 'ygg{as{Q‘^)/N). 
However, in many popular factorization schemes the coefficients in the incomplete, or incorrect 
series for this integrand are much larger than those in 'jgg{as{Q‘^)/N) (helped by the fact that 
many of the early coefficients in 'ygg{as{Q‘^)/N) are zero), e.g. they commonly behave roughly like 
(12In2/7r)"'n“^/‘^, whereas the coefficients in 'ygg{as{Q^)/N) behave roughly like (12In2/7r)”n“^/^. 
Hence, the incorrect ^t{Q^,Qo) can have a dominant effect on the evolution. Under changes of 
factorization scheme the coefficients in the factorization-scheme-dependent series can change by 
amounts similar to their own magnitude, i.e. 'ygg{as{Q‘^)/N) is very unstable under factorization 
scheme changes as can be seen from (2.12). Therefore, the evolution of the structure functions in 
terms of a given input can appear to have a very strong factorization scheme dependence. Once 
again, this is true for the evolution of F 2 {x, Q'^) as well as for Fl{x, Q'^): the influence of the incorrect 
(Q^, Qq) can be more important than that of 4 >q {Q^, Qq) and 72 l(A^, as(Qo)) combined, where 
for the latter the coefficients in the series are again relatively small. 

Simply using an incorrect calculational procedure, such as solving for the parton distribution 
using the renormalization group equations up to some order, and then combining with the coefficient 
functions to a given (combined) order leads to expressions which are not only similar to (4.78)- 
(4.80) with incorrect or missing 'ygg{as{Q‘^)/N), but which have additional factorization-scheme- 
dependent terms. (Solving by using the full known anomalous dimensions and combining the 
resulting parton distributions with all known coefficient functions is even worse.) These will be 
formally of higher order than the terms in (4.78)-(4.80), but again can have very large coefficients 
in the series expansions. This can lead to even more dramatic effects than those outlined above. 
One clear example of such incorrect effects is found in [22]. 

Once 'ygg{N, asiQ"^)) is known in a given scheme and (4.78)-(4.80) can be calculated correctly 
there is no guarantee that the correct <I>£ (Q^, Qq) is not larger than <!>[]" (Q^, Qq)- If this is the case, 
^liQ^iQo) will then have a large, but at least definite, effect. However, because it is a formally 
NLO correction to the structure functions, if introduced then the full set of NLO 
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expressions, both evolution and input factors, must be calculated at the same time. The correct 
calculational method respects renormalization scheme independence as well as factorization scheme 
dependence. As mentioned earlier, this requires many more terms than just Qq). Hopefully, 

the complete NLO expression, as well as being factorization scheme independent, will also cause 
only fairly small changes to the LO expressions. 


Having obtained our full set of leading-order expressions, we can also now examine the dif¬ 
ference between {dF 2 {N,Q‘^)/d ln((5^))o and as{Q‘^)dF 2 ^o{N, Q'^)/d ln((5^), i.e. the LO derivative 
of F 2 {N,Q‘^) and the derivative of the LO F 2 {N,Q‘^). If we differentiate (4.107) with respect to 
cks(Q^) we obtain 




dF2,o{N,Q^) 

das{Q‘^) 


=ois{QI)iIl{ois{QI)/N)(^ 
{FLfl{N) + pLfl{N,Ql) 


7°,(a«(Qg)/iV)r 

- ^^^^ ^F2,o(AI))exp[^+(Q^gg)] 


- iQ^)^^F 2 ,o(N) exp[chr {Q^, Ql ]• 

(4.108) 

This is clearly not exactly the same as (4.100), the difference being due to three additional terms in 
the above expression as compared to (4.100). These are the factor (7°g(Q;s(Q^)/A')/7gg(as(Qo)/-^)) 
in the first term, and the factors (as(Q^)/«s(Qo)) exp[4>j“(Q^, Qq)] in the second term. All 

of these factors are unity at the boundary of the evolution, and the two expressions are therefore 
identical in this limit, i.e the inputs are the same. Therefore, it is the evolution terms which are 
different when comparing (4.100) and (4.108). Writing 


\i%MQl)/N)) 



d 

das{q^) 


Mllg{o^s{q^)/N))da^{q^) 


(4.109) 


and 


a«(Q') 


as(<5o)exp 


— In 


Q«(Qo) Y 


(4.110) 


we see that the terms present in (4.108) but absent in (4.100) are NLO evolution terms. Thus, as 
in the loop expansion, the difference between the fixed-order expression for {dF 2 /dlnQ"^) and the 
In derivative of the fixed-order F 2 consists of terms of higher order. However, in the leading- 
ln(l/x) expansion this difference exists between even LO expressions. Close to the boundary the 
effect of the additional evolution terms in (4.108) is very small, but they clearly become more 
significant for » Qq or << Qq. So, as we might hope, despite the formal difference 
between the derivative of the LO structure function and the LO expression for the derivative of the 
structure function, within the region of expected applicability, i.e. x small and relatively near 
Qq, the two are very similar in practice. As one goes away from the boundary, especially at large x, 
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the terms begin to differ markedlyd^ This highlights the fact that the main region of applicability 
of the leading“ln(l/x) expansion is small x and ~ Qo- 

Finally, we should also make some comment about longitudinal momentum conservation. The 
first moment of the parton distributions is interpreted as the fraction of momentum carried by 
that type of parton. As such it is usually required that + g{l,Q^) = 1, i.e. the total 

momentum carried by the partons is that of the proton. For this to be true for all then, of 
course, + g{l,Q‘^))/dlnQ‘^) = 0, and from the renormalization group equations this 

is true if 

7//(l,a«(Q2)) +7g/(l,a«(Q2)) = q, jfg{l,a,{Q‘^)) + jgg{l,asiQ^)) = 0- (4-111) 

This is assumed to be true for the all-orders anomalous dimensions, and thus momentum will 
be conserved for the all-orders parton distributions. When expanding the anomalous dimensions 
order by order in as{Q‘^), i.e. as in (1.3), it is easy to specify that (4.111) be true for the anomalous 
dimensions at each order, and to define a wide variety of factorization schemes which maintain this. 
It is not difficult to see from §2.2 that this guarantees that the fraction of momentum carried by 
the nth~order parton distributions is conserved at each n. However, it does not necessarily tell us 
anything about the amount of momentum carried by the nth~order inputs, for any particular n. 
Sometimes, all the momentum is designated to be carried by the zeroth order part of the solution, 
but this need not be the case, and indeed, we see no good reason why it should be. Most often the 
input is implicitly assumed to be the all-orders input. In this case it is true that it must carry all 
the momentum, but this method destroys the strict ordering of the solution. 

As always seems to be the case, the situation is not as simple for the leading-ln(l/a;) expansion. 
As can be seen from the matrix in (4.50), the leading-order 7 contains entirely positive entries for 
N = 1, and is clearly not consistent with momentum conservation; /(f((5^) carries a constant 
amount of momentum while that carried by gaiQ"^) is constantly increasing with Q^. In a general 
factorization scheme we assume that there is no reason that working to a finite higher order will 
restore the relationship (4.111). 

Two general methods have been proposed to restore momentum conservation. The first multi¬ 
plies the known y'^’s and 7 ^’s in some factorization scheme by some finite power series in N which 
vanishes at N = 1 [21], e.g. the simplest example is (1 — N). The evolution equations are then 
solved using the whole of these anomalous dimensions and momentum conservation is guaranteed. 
However, this prescription destroys any sense of ordering the solution correctly and is extremely 

Similarly, in the loop expansion significant differences appear between the derivative of fixed order 
F 2 {N,Q^) and {dF 2 {N,Q^)/d\nQ^) at the same fixed order when higher powers of as{Q^) become more 
important, i.e. at low Q^. Hence the differences are small in the expected region of validity of the loop 
expansion, i.e. large Q^. 
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scheme dependent, causing the sort of large factorization scheme variations described above. Also, 
since the full anomalous dimensions have singularities at = — 1 the power series expansion about 
= 0 for given order in does not even converge at A^ = 1. Thus, it seems inappropriate to 
demand that the first few terms in the expansion about A^ = 0 should approximate the correct 
value at A^ = 1. An examination of the first few terms in the expansion of the two-loop anomalous 
dimensions about A" = 0 shows that a truncation of each after a few terms is not at all similar to 
the above prescription. 

An alternative method [20] is to assume that the relationship 


as{Q^)7}f{l,as{Q^)) +-fgf{l,as{Q^)) + as{Q‘^hlf{l,as{Q‘^)) = 0 

asiQ‘^hfg{l,asiQ'^)) +7gg{l,as{Q‘^)) + asiQ^hlgi'^,OisiQ^)) = 0 , 


is satisfied, and to determine the unknown 'y^^{as/N) this way. In fact it has been proved that 
it is always possible to choose a factorization scheme where this is true [36], and suggested that 
this might be some sort of preferred scheme. This will guarantee momentum conservation if one 
truncates the series for the y’s at the 7 ^’s and solves the whole renormalization group equation 
using this truncated 7 . However, this does not lead to a well-ordered solution for the parton 
distributions in any sense. If one were to believe that solving the renormalization group equations 
order by order and combining with the coefficient functions up to a given combined order led to an 
ordered solution for the structure functions in this expansion scheme, as in the loop expansion (and 
as is used in some calculations), then the zeroth-order solutions for the partons do not conserve 
momentum, and even if (4.112) is true then adding (4.50), (4.57) and (4.67) still leaves us with 


E 


-as{Q‘^) 


dasiQ"^) 


7g/(l,«s(Q^))/f (1,(3^) +7gg(l,as(Q^))5'l(l,<3^ 


(4.113) 


as well as terms coming from the difference between 7 ^^ and 7 ^^. Thus, momentum violation will 
not be zero in general, and may be quite large. Thus, it seems that in order to enforce momentum 
conservation strictly within this expansion scheme one must make a guess at the full anomalous 
dimensions (and coefficient functions) in terms of some truncated form of them, and hope that this 
is a good approximation to the full solution. This sacrifices any possibility of making a well-ordered 
expansion for the structure functions. 

Of course, in this paper we advocate that the most sensible approach is to obtain a well- 
ordered solution for the structure functions, and thus definitely eliminate any questions of scheme 
dependence; i.e. regard a correct treatment of the physical quantities as of paramount importance. 
Doing this, there does not seem to be any way to ensure that what we choose to define as the 
parton distributions within our final expressions are such as to conserve momentum. Thus, we 
simply take the hint offered us by the zeroth-order anomalous dimension, and accept the fact that 
momentum is not conserved order by order in this method of expansion. How badly it appears to 
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be violated, however, will depend very much on which factorization scheme we claim to work in. It 
will always be possible to choose one where momentum violation is very small and, if one wishes, 
one may choose to think of this as a “physical scheme” for the partons. However, the structure 
functions themselves will be completely unaffected by this choice, and the real physical relevance 
is therefore rather questionable. Because momentum is not in general conserved order by order 
in this expansion scheme the amount carried by a certain order will vary with Q^. Hence, using 
an arbitrary factorization scheme there seems little reason to demand that the momentum carried 
by the zeroth-order inputs for the partons should sum to one; it seems more sensible to share the 
momentum amongst the different orders in the inputs. It is clear that the momentum carried by 
the zeroth-order parton distributions, for example, will increase quickly with Q^, and thus it makes 
sense when choosing their inputs to choose distributions which sum to less than unity. As higher 
order corrections to the evolution come in, acting to curb this growth in momentum of the parton 
distributions (hopefully countered by increased growth with of the first moment of the structure 
functions coming from the effects of the coefficient functions, leading to the overall behaviour of the 
structure functions being largely unchanged) they can bring inputs carrying positive momentum 
with them. If we could work to all orders then the momentum carried by the inputs would eventually 
sum to unity, and stay at this value for all Q^, but this is of course not possible. 


We can also discuss the relationship between our scheme-independent solutions and the ones 
which would be obtained using Catani’s physical anomalous dimensions. One can solve for the 
LO structure functions using these effective anomalous dimensions in exactly the same way as 
we solved for the parton distributions in the previous subsection. In the same way that we have 
the relationships between the LO anomalous dimensions, (1.7) and (4.41), we have relationships 
between the effective anomalous dimensions (3.13), i.e 

TUa,m/N) = (4.114) 


and 


rl,{asm/N) = - lasiQ^ 


2Nf 


(4.115) 


where the second term in (4.115) is the one-loop contribution to T 2 L{as{Q‘^)/N). Using these rela¬ 
tionships it is straightforward to follow through the steps in (4.50)-(4.77) to obtain the analogous 
expressions to (4.75), (4.76) and (4.77): 


F2{N,Q^) =a,iQ^) 


2,rl^{as{Q^)/N) 




{Fl,o{N) + Fl,o{N,QI) - 


36 - 8Nf 


27 


F2,o{N))x 


exp 


U as{Q'^) 


9/ o\ 4" / o\ (IOIs[Q 

2 ^ af(r) 


asiq‘^ 


- F2,i{N,Qo) + F2,o{N) 


CtsjQl 

as{Q^ 


4 2Nf 
9 67rbQ 


-|- higher order in and/or N, 


(4.116) 
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Fl{N,Q^) ={Fl,o{N)+~Fl,^{N,QI) - )i"2,o(A^))x 
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UasiQ^) Oil{<P^ 
36 - SNf 
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as{q^ 


-das{q^ 
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+ higher order in and/or N, 


(4.117) 


and 


-a 2 (g 2 


dF2{N,Q^^ 

da.siQ'^) 


=a,{Q^)TlAc^,iQ^)/N){FLAN) + FL,oiN,Qt) - 


2 , f36-8Nf 


27 


i^2,o(iV))x 


exp 


L4a4Q2) 


a 2 (g 2 ) 


25^')fiJa,(g2)/iV) 


as(g'^ 


-das A 


-|a,(g 2 )F 2 ,o(iV); 


67r6o Vas(Q^) 

where Fl^o{N,Qq) is a function of rather than 7 °^. 
We see that, once we make the identifications 


.4 


2Nf f asiQo)\ 


+ higher order in and/or N, 

(4.118) 


A,,{asm/N) = rAAsm/N), 


(4.119) 


iAAsm/N) = rAMQ^)m, 


(4.120) 


and 


Agg{as{Q^)/N) + IAfg{as{Q^)/N) 


d ( fci,Asm/N)\\ 
das{Q^)\\ Cl,^o )) 

= fAAs{Q^)/N) - (^^^^)fija,(g2)/iv), 


(4.117) is identical to (4.92), (4.118) is identical to (4.97) and (4.116) is identical to (4.105). Of 
course, the identifications (4.119)-(4.120) are exactly what we obtain from the definitions of the 
physical anomalous dimensions in §3. 

Thus, we are able to reach these expressions for the structure functions somewhat more directly 
by using the physical anomalous dimensions, and do not have to worry about problems with factor¬ 
ization scheme dependence (though we do have to calculate the physical anomalous dimensions in 
terms of known coefficient functions and anomalous dimensions of course). Once we have obtained 
these expressions using the physical anomalous dimensions we may then separate each of the terms 
into input parts and evolution parts (where more of this has already been done automatically when 
using the physical anomalous dimensions) and keep the most leading parts, obtaining once again 
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the LO expressions (4.93), (4.100) and (4.107). So, using the physical anomalous dimensions leads 
us in a rather more direct manner to the correct leading-order expressions. If we were to work to 
higher orders, the amount of simplification obtained by using the physical anomalous dimensions 
rather than working in terms of parton densities would increase significantly. However, we stress 
that one will always automatically obtain factorization-scheme-independent answers by working 
to well-dehned orders in physical quantities even when working in terms of partons. Also, even if 
one uses the physical anomalous dimensions, care is still needed to obtain expressions which are 
consistent with renormalization scheme dependence. 

In this section, we have derived well-ordered, factorization-scheme-independent expressions 
for structure functions in the leading-ln(l/a;) expansion (which should be useful at small x and 
moderate Q^) up to the order which is useful at present. This expansion does, however, sacrifice 
any attempt to describe the structure functions at large x at any reasonable distance from Qq (we 
will discuss later what large x and small x turn out to be), in the same way that the loop expansion 
should show signs of failing at very low values of x. We would hope there is some expansion scheme 
which will be useful at all values of x. In the next section we will show that there is indeed an 
expansion scheme which satishes this criterion, and argue that it is the only really correct expansion 
scheme. 

4.4. The Renormalization—Scheme—Consistent Expansion. 

In order to devise an expansion scheme which is useful at both large and small x we would a 
priori expect that we would need to use the known anomalous dimensions and coefficient functions 
at low orders in both and in the leading-ln(l/a;) expansion. There have already been various 
methods along these lines, and the phrase “double leading expansion” was coined in [20]. However, 
these methods have all suffered from scheme dependence. As when deriving our expressions for the 
structure functions when using the leading-ln(I/x) expansion we will ensure that we obtain results 
which are invariant under changes of factorization scheme and, as a stronger constraint, demand 
complete consistency of our expressions for physical quantities with renormalization scheme invari¬ 
ance (which in itself automatically guarantees factorization scheme independence). Consequently, 
our approach will be rather different from those used previously, and the results and conclusions 
will also be somewhat different. 

To begin, let us consider what we have meant by “consistency with renormalization scheme 
dependence” so far in this paper. In both the loop expansion and the leading-ln(I/x) expansion we 
demanded that once we had chosen a particular renormalization scheme and chosen to work to a 
particular order in this renormalization scheme then we would include all terms in our expressions 
for the structure function which were of greater magnitude than the uncertainty due to the freedom 
of choice of renormalization scheme (i.e. the uncertainty in the definition of the coupling constant). 
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and no others. In both cases the leading-order term consisted of the lowest-order inputs multiplying 
the lowest-order evolution terms. If working with the n-loop coupling constant in a particular 
renormalization scheme the uncertainty in its definition is of order Thus, the uncertainty 

of the input or evolution when working to nth"Order is the change in the leading-order input or 
evolution if the coupling changes by Sag = eag~^^, i.e > q;s( 1 + ). Hence, the uncertainty in 

the whole structure function is of the order of the change of the leading-order part under such a 
change in the coupling. Therefore the nth“Order renormalization-scheme-independent expression 
includes all complete terms smaller than this change. 

This definition does give us a well-defined way of building up an ordered solution to the struc¬ 
ture functions, but relies upon the definition of a given expansion scheme. It leaves an ambiguity 
about how we define the leading-order expressions and in how we define the order of terms com¬ 
pared to this leading-order term. Our two examples, i.e. the loop expansion, where the size of a 
term is determined simply by its order in and the leading-ln(l/x) expansion, where ln(l/x) is 
put on an equal footing to a^, are just the two most commonly used examples of expansion schemes 
(even though the ln(l/x) expansion has not previously been presented in the quite same way as in 
this paper). Both have potential problems; in the former one does not worry about the fact that 
the large-ln(l/x) terms can cause enhancement at small x of terms which are higher order in 
and in the latter one does not worry about the fact that at large x, especially as increases, it is 
the terms which are of lowest order in ag that are most important. Hence, one would think that 
both have limited regions of validity. 

The shortcomings of these two expansion schemes come about because, even though any given 
order contains no terms which are inconsistent with working to the same given order in a particular 
renormalization scheme, in neither case does it include every one of the terms which are consistent 
to working to a given order in the renormalization scheme. In each expansion scheme some of the 
terms appearing at what we call higher orders are not actually subleading in ag to any terms which 
have already appeared. Thus, despite the fact that for a given expansion method these terms are 
formally of the same order as uncertainties due to the choice of renormalization scheme, they are 
not terms which can actually be generated by a change in renormalization scheme.^® 

In order to demonstrate this point more clearly we consider a simple toy model. Let us imagine 
some hypothetical physical quantity which can be expressed in the form 

OO OO OO OO 

,ag{Q^)) = Y,<^s{Q^) Yl = Y<iQ^) E 

m=l n= — 7n i=0 j=l — i 

where the expansion in powers of N about = 0 is convergent for all N. The first way of writing 
H{N,as{Q‘^)) as a power series corresponds to the loop expansion, where we work order by order 

Similarly, they cannot be generated by a change in renormalization scale. 



as{Q 

N 


(4.122) 
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in m, out to m = k, and use the /c-loop coupling. The second corresponds to the leading~ln(l/x) 
expansion where we work order by order in i, out to i = I, and use the (Z + l)-loop coupling. 
Let us, for a moment, consider the LO expression in the loop expansion, 

The coupling is uncertain by and hence the uncertainty of the leading-order expression 

(i.e. the change due to a change of the coupling) is ~ 0 (g(Q‘^) ^in-ZV*^. There is no change 

with powers of N less than —1, and hence any such term is not really subleading. Similarly, the 
uncertainty of the leading-order expression in the leading ln(l/x) expansion contains no terms at 
first order in (or with positive powers of N), and such terms are not really subleading either. The 
full set of terms contained within the combination of both leading-order expressions is genuinely 
leading order, and is therefore renormalization scheme independent by definition. 

Perhaps the best way in which to write our expression for H{N, as{Q‘^)) in order to appreciate 
these points is 

OO OO OO OO 

H{N,a,{Q^)) = Y, + E E (4.123) 

m= —1 n=l m=2 n=m 


i.e. as an infinite number of power series in as{Q‘^), one for each power on N. Each of these 
series in as(Q^) is independent of the others, and the lowest order in as(Q'^) of each is therefore 
renormalization scheme independent and part of the complete LO expression for H(N,as(Q^)}. 
The full LO expression for H{N,as{Q‘^)) is therefore 


Ho{N,as{Q^)) = E + E 


m= — l 


m=2 

OO 


= a 


:(Q") E «0niV” + E^O 

1=2 


n= —1 


Ois (^Q 
N 


2\\3 


(4.124) 


Hence, the combined set of terms considered LO in both our expansion schemes comprise the 
full set of renormalization scheme invariant, and thus truly leading-order, terms. By considering 
H[N,as{Q‘^)) written in the form (4.123), and considering how the coefficients in the expression 
must change in order to make the whole expression invariant under a redefinition of the coupling 
constant, as(Q^) + C>(a^((5^)), we see that the nth-order expression for H{N, o-siQ'^)), 

which should be used with the n-loop coupling constant, consists of the sum of the first n terms 
in each of the power series in as((5^). Thus, the full nth-order expression always consists of the 
nth“Order expression in the loop expansion plus additional terms with inverse powers of N greater 
than n. 

Similar arguments have already been applied to the anomalous dimensions and coefficient 
functions, for example [18] and particularly [20]. The latter claims that one may, but need not, use 
expansions of the above form for the anomalous dimensions and coefficient functions, and moreover, 
in practice expresses the terms beyond fixed order in as functions of x/xq for x < xq and sets them 
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to zero otherwise (xq < 1, and is in general simO.l), thus reducing their effect (see also [28]). Here we 
take a strong, inflexible viewpoint and insist that the complete renormalization-scheme-consistent 
expressions, with no artificial suppression of leading“ln(l/x) terms, must be used. Furthermore, 
and very importantly, the expressions used must be those for the physical structure functions, 
not for the factorization-scheme- and renormalization-scheme-dependent coefficient functions and 
anomalous dimensions. 

When considering the real structure functions the situation is technically a great deal more 
complicated than our toy model, but the principle is exactly the same. This can be seen by 
examining the the LO expressions for the structure functions in the two expansion schemes already 
considered. There is some overlap between the LO expressions for the structure functions when 
using the loop expansion and when using the leading-ln(l/x) expansion, but each contains an 
infinite number of terms not present in the other. However, we were previously happy to use the 
one loop coupling for both expressions. The uncertainty in the definition of this coupling is O(a^). 
Considering the change of each of our leading-order expressions under a change of coupling of 
0(ag), the changes in the expressions are rather complicated. However, it is not too difficult to 
see that, as with our toy model, the change in the LO structure functions in the loop expansion 
contains no terms in the LO expressions in the leading-ln(l/a;) expansion, and vice versa. Thus, 
none of the terms contained within each of the LO expressions are generated by uncertainties at 
higher order in the opposing expansion scheme. Therefore, they should really all be regarded as 
genuinely LO, and be included in the full expressions for the structure functions which use the 
one-loop coupling constant. 

Hence, as with the toy model, there should be some combined expansion-scheme-independent 
expressions for the structure functions which we can genuinely call the “leading order” expressions. 
Since these expressions will contain all the parts of the one-loop expressions, and also contain 
leading-ln(l/x) terms as well, they should be able to describe the data over the full range of 
parameter space (except very low Q^, of course), as we would like from our correct LO expressions. 
We shall now demonstrate how we obtain these expressions. 

There are two main complications when considering structure functions in comparison to our 
simple toy model. One is that the structure functions are combinations of perturbative evolution 
parts and input parts (which are viewed as partly perturbative with nonperturbative factors), 
rather than one simple power series in as{Q‘^). The other is that in general the physical anomalous 
dimensions, out of which the perturbative parts are constructed, are nonanalytic functions which 
cannot be expressed as power series about a particular value A^o for all N. The physical anomalous 
dimensions have singularities at = 0 (in the case of the singlet structure function only), and also 
at negative powers of N (as well as possible Q!s((5^)-dependent nonanalyticities due to resummation 
effects, e.g. the branch point in T® ^(Ai, as(Q^)) at A^ = X{Q‘^)). We will deal with this second 
complication first. 
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Let us consider the perturbative parts of the expressions for the structure functions. The 
singularities at negative integer values of N mean that we cannot write any physically meaningful 
quantity as just a power series about iV = 0 (or about {N + 1) for the nonsinglet case). Any such 
power series expansion will have a radius of convergence of unity, whereas the physical moments 
of any x-space quantity will exist for all real N above some minimum Nmin, which depends on 
the asymptotic form of the structure functions as x ^ 0. A series expansion which applies over 
this whole range of N does not exist: the valid expression must include the nonanalytic functions 
explicitly. This seems to make it impossible to order the moment-space solution in such a simple 
way as in (4.122). 

In order to overcome this problem let us consider making the inverse transformation of some 
physically relevant perturbative function A{N, as{Q‘^)) to x-space. The inverse of the Mellin trans¬ 
formation (2.1) is 

1 pc-\-ioc 

A{x,as{Q^)) = — x-^A{N,as{Q^))dN, (4.125) 

J c—ioc 

where the line of integration is to the right of all nonanalyticities. Making the substitution ^ = 
ln(l/x) this becomes 

1 pc-\-ioc 

A(x,as((5^)) = — / ex.p[^N]A{N,aa{Q‘^))dN. (4.126) 

J c—ioo 

Since A(N,as{Q‘^)) has, in general, singularities for all nonpositive integers, this whole integral 
may be evaluated by performing an infinite series of integrals, each with a contour centred on a 
given singularity, and not extending as far as unity from this singularity, i.e. not reaching any of 
the other singularities. Within each of these contours the function A{N, as{Q‘^)) may be expanded 
as a power series about the singularity, i.e. we may write 

1 °° f 

A{x,as{Q‘^)) = — '^x‘^ e-xp[^{N + n)]An{{N + n),as{Q'^))d{N + n), (4.127) 

n=0 

where An{{N + n), as{Q‘^)) denotes A{N,as{Q‘^)) expanded as a power series about N = —n. The 
integrals will produce functions of which do not sum to integer powers of x, and hence each of 
the integrals in (4.127) will be independent and physically relevant in its own right. 

So this is the solution to our problem of how to order the moment space expressions for physical 
expressions which are related to the real structure functions as powers series in and N. We must 
consider the complete moment-space expression as an infinite number of expressions of the form 
(4.123), each one having power series expansions in terms of {N + n), where n = 0 ^ oo. The 
expression for each n is then related to the part of the x-space expression behaving ~ x”. Of 
course, in practice, unless we want to examine the details of the perturbative calculation of the 
structure function for x very close to 1, we can ignore all n greater than a finite, relatively small 
constant. 
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Thus, when we calculate the expressions for the perturbative part of the singlet structure 
functions, we will only be concerned about any LO terms beyond lowest order in for the specific 
case of n = 0. For n > 0 we take the whole LO expression to be the one-loop expression. The 
terms we ignore by making this necessary decision are those which are LO in ln(l/x) at first order 
in X. Although these terms grow like a™ ln^™'~^(l/x), as opposed to a^ln”*(l/x) at zeroth order 
in X, there is no evidence that their coefficients are any larger than those for the zeroth-order- 
in-x logarithms. Since the resummed terms at zeroth order in x only begin to make an impact 
as X falls to ~ 0.1 (as we will see), and only become dominant for x much smaller than this, the 
effect of terms like xa^ ln^"*~^(l/x) should be extremely small in comparison. Indeed, the effect 
of those terms of the form xa^ ln^™'~^(l/x) which are actually known, i.e. for m = 2, can indeed 
be shown to be negligible. In a similar manner, we will only consider the one-loop expressions for 
the nonsinglet structure functions in practice: the other LO parts of the expressions again lead to 
small-x enhancement of the form xa™ ln^™'~^(l/x), which is very small compared to the leading 
singlet small-x enhancement, and there is only detailed data at very small x for the total structure 
function.^® 

Hence, we only really need to consider calculating a full LO renormalization-scheme-consistent 
(RSC) expression for the perturbative contributions to the structure functions expressed as power 
series in ag and in N about = 0, as for our toy model. However, we now have to return to 
our first problem, i.e. the fact that the structure functions are expressed in terms of both inputs 
and evolution parts. Using the results we have already obtained in the earlier subsections it is not 
too difficult to construct the full LORSC expressions for the inputs and for the evolution parts 
of structure functions. For the case of the nonsinglet structure functions the construction of the 
LORSC expression is then just the product of these two terms. Let us discuss this as an example 
first. 


4.5. The Renormalization—Scheme Consistent Nonsinglet Structure Functions. 


We consider a nonsinglet longitudinal structure function. For the nonsinglet structure func¬ 
tions the physical anomalous dimensions contain no singularities at iV = 0, so the leading-ln(I/x) 
behaviour comes from singularities at N = —I. Expanding about N = —1, the full LO physical 
anomalous dimension can be written in the form, 
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amiN + iy 


+ bmiN + I) ^ 
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\{N + iy 


m 



(iV + I)2 


(4.128) 


Also, the full leading-order-in-ln(I/x) physical anomalous dimension is not yet known for the nons¬ 
inglet structure functions, as will be discussed below. 
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The first sum is just +1)> one-loop anomalous dimension expanded in powers of {N +1). 

The second sum contains the leading singularities in (N + 1) for all other orders in as{Q‘^)- The 
final sum is included because, despite the obvious fact that it is a power of {N + 1)“^ down on 
the second sum, a series of this form cannot be created from the second sum by a change in the 
definition of the coupling of 0(ag(Q‘^)). Therefore, the third sum is not subleading in aa{Q‘^) to 
the second sum, and must be renormalization scheme independent. Integrating (4.128) between 
Qq and Q^, and including the overall power of asiQ"^) for the longitudinal structure function we 
obtain the leading-order evolution 


Us{Q^)J 


exp 



\N + 1 



/ \ daa{q^) ' 

\{N + l)y JboaUq^)/ 


(4.129) 


In the loop expansion the lowest-order input for F^^{N,Q^) was as{Qo)C^f 
where F^q{N) is a nonperturbative factor and C^fi{N) has an expansion in powers of {N + 1) 
beginning at zeroth order. Factoring out the nonperturbative part, our lowest-order input is a{Qo) 
multiplying a power series in {N + 1) which starts at zeroth order, i.e C^fi{N + 1). In order 
to construct the full LO input we must consider how the evolution term (4.129) changes under a 
change in starting scale Qq ^ (1 + 6)Qq, and therefore how the input must change in order to 
compensate for this. The change in the first term is oc Q!s(go)(^L o + 1) ~ !)• This can be 
absorbed into a change of the input at order a‘^{QQ), and hence the second-order-in-Q!s(go) input 
has a part cx (r^Q;(A^ + 1) — l)C'^f;(A^ + 1). The power expansion of this expression in terms 
of {N + 1) begins at order —1, and a term of this type cannot be generated by a change of the 
order-as(gQ) input under renormalization scheme changes. Hence, this part of the a‘^{Qo) input, 
i.e (X a^{QQ)a-iC^f i{N = —l)(iV + 1)“^, belongs to the LORSC input. The rest of the Ci(a^(gg)) 
input is genuinely subleading to the 0{as{QQ)) input and is renormalization scheme dependent. 
Extending this argument, and considering the form of the input required to compensate for the 
change of the whole of (4.129) under a change in starting scale (as with goiNjQ^) in subsection 
4.2), we can see that the full LORSC input is 


»siQl)Fi^oHN){C^LiN + 1 )+C^Un = -l){HQl/ADaaiQl)x 


jjj. 

(4.130) 

where bo = a_i and cq = 0. The first term is just the lowest-order input in the loop expansion, 
while the second includes all the leading-ln(l/x) terms in a simple form which is compatible with 
making the full expression invariant under changes in starting scale Qq. This second term only 
depends on the part of the one-loop coefficient function at zeroth order in (iV + 1); the higher- 
order-in-(A^ + 1) parts multiply the part consisting of the leading-ln(l/x) terms to give the type 
of terms which can be generated from (4.130) by changes in the renormalization scheme. 


+ ^ ^ (Cm + ^ ^ ^ bnbm-n ^^{Qo/ 


«s((3o) 
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Hence, the full LORSC expression for the nonsinglet longitudinal structure function, expanded 
about N = —1, is 



(4.131) 


As we have already argued for the singlet case, as far as the expansion about the singularities at N to 
the left of the rightmost singularity are concerned, we may as well just take the one-loop expression. 
In practice, we will only use the one-loop expression for the nonsinglet structure functions for all 
the singularities in the anomalous dimensions. This is because of the phenomenological reasons 
given at the end of the last subsection, and also because of lack of knowledge of the full physical 
anomalous dimensions. In the MS scheme the terms ~ a'^{N + parton anomalous 

dimension are all known [37]. It does not appear as though the coefficient functions contribute 
to these sort of terms in the physical anomalous dimension (though there is no formal proof of 
this), and hence, it is believed all of these terms are known. However, there is little knowledge 
yet of the terms of the sort ~ a'^{N + 2 )^™'“^. We have argued that these are an intrinsic part 
of the LORSC expression for the nonsinglet structure function, and they should be calculated and 
included in order to give a true indication of the effect of leading-ln(I/x) terms. Hence, we believe 
that calculations of the nonsinglet [23] [38] (and polarized [39]) structure functions which claim to 
include leading-ln(I/x) corrections are incomplete, even at leading order, until the terms of the 
form a'P{N + l) 2>^-2 known. 


We now consider the nonsinglet structure function F^^{N,Q^). This leads us back to our 
previous question of whether we should use the full RSC expression for F 2 {N,Q‘^) or that for 
{dF 2 {N,Q‘^)/d InQ^). In order to illustrate the difference between the two, and help us make our 
choice, we consider the simpler nonsinglet case before the singlet case. We also pretend for the 
moment that the inverse powers of {N + I) in the expressions do not increase as the power of 
increases, i.e. there is no small-x enhancement at higher orders in a^. This being the case, the LO 
term in the evolution is just 


as(Qc 


(4.132) 


Of course, the input may be written as a power series in UsiQ^), as we saw in §4.1, and is of the 
form 


F27iN^ Ql) = + ^siQl) HQl/A77j^f{N + I) + higher order in a«(Q^)]. (4.133) 
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Therefore, we have the problem that as well as the lowest order a*-, and hence Qo“dependent part 
of the input there is also the “sub-lowest-order”, Qo'i^dependent part. These two terms are clearly 
of different order, but under a change in renormalization scheme, and hence in the definition of 
the coupling, both remain unchanged and both should therefore appear in the LO definition of the 
structure function. This mixing of orders seems rather unsatisfactory, and comes about because for 
the structure function still exists in the formal limit of a* ^ 0, i.e. the parton model 
limit, being equal to the simple Qg-independent function F^q{N). Hence, it is not a perturbative 
quantity in quite the sense way as F^^ or (dF^^{N,Q‘^)/dlnQ‘^), both of which vanish in 
this limit. 

So in our simplified model the LORSC expression for F^^{N, Q^), obtained by combining the 
LO input and evolution, is 




as{Ql)HQl/ADr%^\N + 1)] 


a 


m 


a 






(4.134) 


This consists of two parts which are clearly of different magnitude, i.e. one is a power of as down on 
the other with no small~x enhancement. This seems against the spirit of a well-ordered calculation. 
Indeed, the second part of this LO expression is of the same order of magnitude as part of the LO 
input multiplying the NLO evolution, i.e. 




as{Q^))t^A^\N + l) 


( as{Ql) 


fyNS 
^ 2,0,i 


(7V+1) 


(4.135) 


Even when we take into account the higher inverse powers of (A^ + 1) at higher powers of ag, (4.134) 
will be part of the LO expression for F^^{N, Q^), and this unsatisfactory behaviour remains. It is 
also clear that the same effect will be seen for the singlet structure function. 

If we instead consider {dF^^^{N,Q^)/d InQ^) (again ignoring small-rc enhancement for the 
moment) the situation improves. In this case the full expression is 


dFi^^{N,Q^) 

dlnQ 2 


=r 2 ,ns{N + 1, as{Ql))Fi^^{N, Ql) x 




ns(,N + 1 , aa{q‘^)) — 


d In q‘ 


■ln{T2,Ns{N + l,as{q‘^))) 


dlnq'^ . 
(4.136) 


Hence, the input may be written as 


/dFi^,^N,Q^ 
V dlnQ^ 


Ql 


- pNS 
~ ^2,0 


(iV)[a,(g2)r2.04(^r+ 1) + 
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(4.137) 



Of course, the 0{asiQQ)) piece is renormalization scheme independent by definition. The 
0 (q;^((5q)) piece is renormalization scheme dependent is renormalization scheme de¬ 

pendent) in order to absorb changes in the 0{as{QQ)) piece under a change in the coupling of 
0 (q;^((5o))- So this time we have a LORSC input which is of a given order in UsiQ^). The full LO 
expression for {dF 2 {N,Q^)/dlnQ'^) is then 


fdFi^,^{N,Q^)\ 

V dlnQ^ J, 




NS 


asm) 


(4.138) 


This is rather more satisfactory than the renormalization-schemeOndependent expression for 
F^^{N,Q^) itself, and hence we choose {dF 2 {N,Q^)/dlnQ^) to be the perturbative quantity we 
calculate, in both the nonsinglet and singlet case. 

If we wish to calculate the structure function F 2 {N,Q‘^) itself to a given order we will do this 
by integrating the given order expression for {dF 2 {N, Q‘^)/dlnQ'^) between and Q^, and adding 
it to F 2 {N,Q‘q) evaluated to the same order. For example, in our simplified nonsinglet model we 
would integrate (4.138) and add this to the explicitly written part of (4.133). This results in the 
effective LO expression 


F2%cA^,Q'")=Fi:o{N) 


NS/ 


a 


■{QV> 


a 




- 2 , 0,1 


(N+l) 


+ asiQ't,) HQo/A^A^NViN + 1)]. (4.139) 


This is perhaps more sensible than (4.134), since the two terms are now of a more comparable size. 
Moreover, this expression is more stable under changes in Qq. 

Of course, this whole discussion of FA^{N,Q‘^) has been rather simplified by the assumption 
that the higher-order-in-Os terms in the physical anomalous dimension do not contain higher sin¬ 
gularities in (N -|- 1). Recognizing that they do, we obtain expressions which are more complicated, 
as in the case of FA^{N,Q^). With a little work it is possible to see that the full leading-order 
expressions are 




1 


d lng2 

(pm + {N + l)Cm) 


AT + l' 


m=l 


as{Ql) 
{N + lf 


"" «-i«s(go) 1 /p 2 //.ivs\ 1 ( as{Ql) 


( as{Ql) 

wm 


exp 



1 

N + l 



f asiq^) 
({N + iy 


\ das{q^) ' 

)boam_ 


( 1 + Em=i mm m/jN + 1 )^-) \ \ 

V 1 +Er=i bmiam m/{N+im)) 

(4.140) 
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and 


j^NS 

^2,RSC,0 


{N, Ql) = ( 1 + \n{Ql/A^^)aM 


-pNS 

^ 2 , 0,1 


{N + l)+ 


1 SF I f ^s{Qo) 

iV + 1 ^ WiV + 1)2 

m=l ^ ^ ^ 


m oo m 

+ ^ (cm + 5 ^ Kbm-n ln((5o/yl^^)) 
m =0 n .=0 


( as{Ql) 
\{N + iy 



(4.141) 

where for F^^{N,Q^) the coefficients in the series in {N + l)“^are not necessarily the same as 
for the longitudinal structure function, hence the slightly different notation. However, a_i = d_i 
and, if the hypothesis that in the MS renormalization and factorization scheme all the contribution 
to the a'^/{N + i)2>^-i singularities in the physical anomalous dimensions is due to the parton 
anomalous dimension, bm = bm as well. The Cm are not equal to the Cm in general though, so 
there is no guarantee that q{x, Q"^) and F^^{x, Q^) will behave in the same way in the 

small-x limit. 


When all the bm and Cm are known they can be used to present an argument for the form of the 
small-x behaviour of nonsinglet structure functions, e.g. to explain any discrepancy between the 
input power at small x for these structure functions and predictions from Regge physics. This seems 
to us to be an interesting project, and we look forward to the calculation of these coefficients. Until 
this happens our discussion of the LORSC calculation of the nonsinglet structure functions is rather 
academic. However, it has enabled us to discuss many of the issues in a simpler framework than 
if we had gone directly to the singlet structure functions. We will discuss these singlet structure 
functions next. 


4.6. The Renormalization—Scheme Consistent Singlet Structnre Functions. 

When calculating the singlet structure functions we cannot just construct the complete LO 
evolution and input and combine these to obtain the LO expression because the evolution mixes the 
two different structure functions. Each of the component parts of the LO expressions for F 2 {N, Q^) 
and Fl{N,Q‘^) (we omit the superscript S in this section) must consist of LO input parts and 
evolution parts, but it is not obvious what these are. In order to find the full LORSC expressions 
for the singlet structure functions we will have to work in steps. We will consider only the full LO 
expression with the perturbative factors expanded about the particular value of iV = 0 (the non- 
perturbative inputs are the full nonanalytic expressions for Fl{N) and F 2 {N)), and the simplest 
way to proceed is to work directly with the physical quantities, solving the evolution equations in 
terms of physical anomalous dimensions and structure functions. We have already proved that in 
the loop expansion the LO expressions only depend on the one-loop physical anomalous dimen¬ 
sions, and in the leading-ln(l/x) expansion the LO expressions depend only on r°^(Q;s((5^)/A'"), 
r® 2 (as(Q^)/.^), ^ 2 Li^siQ‘^)/N) and T 22 iois{Q‘^)/F^)- Hence, it is only the combination of these 
anomalous dimensions which is considered in our solution. 
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We cannot simply write the physical anomalous dimension matrix 


f as{Q^)T^£[{N) +Tl^{as{Q^)/N) 

V a. mrli (N) + a, (Q2)ri^ (a, (Q^)/N) 


as{Q^)T%{N) + Tl,{a,{Q^)/N) \ 

a^mrliiN) + asmTUa,{Q^)/N)) 


(4.142) 


(where T\]^{as{Q‘^) /N) = T\i^{as{Q'^) /N) with the one-loop component subtracted out, etc.), 
and solve the renormalization group equations. With this anomalous dimension matrix there is no 
simple closed form for the solution of these equations, and the full solution contains terms which 
are not properly of leading order. We must choose some way of solving for the structure functions 
systematically which enables us to extract the true LO behaviour in as simple a manner as possible. 

In order to do this we take account of the fact that the one-loop solutions for Fl{N,Q‘^) 
and {dF 2 {N,Q‘^)/dlnQ'^) must be part of the complete LORSC solutions. Hence, we split our 
anomalous dimension matrix up into the form 

„,o2d’’°4W+ RiWl + l r»jo.(02)/jv) rij(o,(Q2)/jv) 

T°i{N))*\a.{Q^)TUa,{Q-^)/N) a.(Q^)vUa.{Q^)/N) ] ' '' ’ 


and solve by treating the second matrix as a perturbation to the first. Doing this we obtain the 
one-loop solutions as the lowest-order solutions and can systematically calculate corrections to 
this, extracting the parts of these “corrections” which are leading order. 

So, first let us consider the solution to the renormalization group equation 


d (f^,:\n,q‘^)\ 2 ^r°’i(iv) r%{N)\ (fI’\n,q^)\ 

dlnQ2 VF°’'(iV,Q2); ’\rliiN) r°^'(iv) j [f^’\n,q^) ) 


(4.144) 


with boundary conditions F^’\n,Qq) = Fl{N) and F 2 ’\n,QI) = F 2 {N). We may write the 
solution for the longitudinal structure function as 


F°’'(iV,Q^) = F^ 


01 , 


-(N) 


a 


■{Q 


a 




^0,1,+ 


(N) 


+ F°’'’-(iV) 


a 


m 


a 


:(Q^ 


-{N) 


(4.145) 


where f*^’^’"*“’“(A^) are the two eigenvalues of the zeroth-order physical anomalous dimension matrix 
(which are the same as the eigenvalues of the zeroth-order parton anomalous dimension matrix), 
and f2’^'~^{N) + f2’'''~(N) = Fl{N). Having chosen to write the lowest-order solution for the lon¬ 
gitudinal structure function in this way we may then write the lowest-order solution for F 2 {N, ) 
as 


F°’'(1V,Q2) =e+(lV)F°’'’+(iV)(^^^j +e-(iV)F°’'’-(iV) 




where 


e+(iV) = 


roi,+(Ar)_rO.[(Ar) 


r°’2(iv) 


e-(iV) = 


asiQ'^ 

ro.^.-(Ar)-r»’[(Ar) 

T%iN) 


, (4.146) 


(4.147) 
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and e+(iV)F°’'’+(iV) + e-(iV)F°’'’“(7V) = F 2 {N). In practice 

e+(iV) = N/6 + 0{N% e-{N) = + 0{N), (4.148) 

and 


po,i,+^N) = FLiN)-^^^^F2{N) + 0{N), F°’'’-(iV) = ^^—^F2{N) + 0{N). (4.149) 

It is then simple to see that 

F°’'’+(iV) = e+(iV)F°’'’+(iV) = ^(^Fl{N) - ^^(iV) + ©(iV^), 

F^’^’-{N) = e-{N)Fl’^’~{N) = F2{N) + 0{N). 


The first correction to the one-loop solution may be obtained by solving the equation, 

^0,1 ( -\T\ T'O,/ , 

d ln( 


d ( FE\N,Q^)\ ^ 2. TAiN)\ (Fp{N,Q^)\ 

[nQ 2 (iv,g2) ) ArAiN) r°^'(Ar) ) \Ff{N,Q^) ) 


+ 


fA{N,Q‘^ 


TAia^m/N) -^J2z^YA{aPQ^)/N) 

as{Q^)TA{aPQ^)/N) -^^^F^as{Q^)TA{as{Q^)/N) ) \fA{N,Q^)) ’ 

(4.151) 

where we have used the relationships in (4.114) and (4.115) in order to simplify the second matrix. 
We proceed as follows. First we define the vectors 


= («+(%)• ^^('^>=(e-(‘iV))- 


F 


■ ^ “ VF|i(iv,g 2 ) j ’ 

(4.152) 


and write 


F^^ {N, g^) = e+ (iV)F^i’+ {N, Q^) + eT (iV)F^h- (jv, g2). 
We also define projection operators p'^{N) and p~{N) by 


p+(iV) •e+(iV) = 1, 
p~{N)-e+{N) = 0, 
which in practice gives 

P+(iV) = 8n}-36 ) +0(N), 


p+(fV) •e“(A^) = 0, 

p-(N)-e-(N) = l, 


p-(N) = ( 27 ]+0(N). 

8Nf-36 


(4.153) 


(4.154) 


(4.155) 


Multiplying (4.151) byp"'“(A^) now leads to the straightforward first-order differential equation 


dF^^’+(N, g 2 ) 
d lng2 

p+(N) ■ 


a,(g 2 )rO’'’+(iv)F"i-+(iv,g 2 ) + 

/ r8^(a,(g2)/iv) _(25gfMr6^(«^(g2)/iv) \ 

[aPQ^)rA(aPQ^)/N) -(^^^aPQ^)rA(aPQ^)/N) J 


fFA(N,Q^)] 

[fA(n,q^)J' 

(4.156) 


68 



We also write the zeroth-order solution in the form 

= e+{N)F°’^’+{N,Q^) + e-{N)F^’^’-{N,Q^), (4.157) 

where F^’’'’~^{N,Q^) = {N){as{Qo)/as{Q^))^°''^'^^^^ and similarly for {N,Q‘^). Doing 

this (4.151) becomes 

dF^^’+{N, Q 2 ) 


d lnQ 2 


= as{Q^)r^’^’+{N)F‘^^’+{N,Q^)+ 


p+{N) 
e+(iV)F°’'’+(iV)( 


rL(a«(Q")/iV) -^^^^rl^{as{Q^)/N) ] 


®s(Qo^ 


+ e-(iV)F°’'’-(iV)( 


a,(Q2)y ■ - V ; L V \a,{Q2) 

(4.158) 

This equation can now be solved by using the power series expansions of (N) and (N) in 
terms of N. Only a small part of the overall solution contributes at LO. Considering the contraction 
of the matrices between the two vectors in (4.158) we obtain 


p+(A^) • 


■^55^™rL(a,(Q2)/iv) 


as{Q^)Tyasm/N) -(2^^a,(Q^)Tl^(asm/N) 


e+(iV) 


(4.159) 


= T}^i^{as{Q )/N) + higher order, 


and 


p+(Af) 


^25z^rL(a,(Q2)/^) 


rl^iasm/N) 

as{Q^)Tl^{as{Q^)/N) -(2^^as{Q^)Tl^{as{Q^)/N) 


\e-{N) 


= 0 a,iQ^) 


CisiQ 

N 


2\\^ 


(4.160) 


Therefore (4.158) becomes 

dF^^’+{N,Q^) 


d lng 2 


= a,{Q^)T^’^’+{N)F^^’+ {N,Q^) + 

2 s f°'‘’ + (iV) 


Tl^ia,m/N)Fl’^’+{N) 


as(Qo ) 

OisiQ^ 


(4.161) 


+ higher order. 


Thus 


r-a.(Qg) rl^(asiq^)/N) 

>C^s(Q^) 


F^h+(^r,g2) ^ r - ^ ) da,{q^)Fl'^’+{N) 

J oc 


al{q 


dsiQ 

Ois{Q‘^ 


2 . p.'. + (iV) 


(4.162) 


+ higher order. 


This method of solution demonstrates the particularly nice feature that the leading-order “cor¬ 
rections” to the one-loop solution which are proportional to e^{N) are completely independent of 
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^o,i, Since this is due to the form of the “correction” matrix contracted between p'^{N) 

and e~{N), we can write the equation for the nth”Order correction to Q"^) as 


dF^^’+{N,Q^) 

dlnQ 2 


a,(Q2)rO-«.+(iV)F"".+ (iV,g2) 


+ p+(iV) 


/ rL(a.(Q 2 )/iV) 

\as{Q^)Tlj^{as{Q^)/N) 




and precisely the same argument goes through each time. Therefore, 


(4.163) 


dF^^’+{N, Q2 


= a,{Q^)T^’^’+{N)F^'^’+{N, Q2)+r0 ^(a,(g2)/^r)F"(—i)-+(iV, Q^)+highei order. 

(4.164) 


d lnQ2 

Thus, it is easy to prove by induction that 


F^-’+(iV,Q 2 ) = 


Oil{q 




+ higher order. 


(4.165) 


This leads to the straightforward expression for the whole of the leading-order part of F^{N, Q^) 

^s{Ql) fojiV,a,(g2)) 


p5C,oPP) = P’^(P 


asiQ‘^ 


exp 




Uas{Q^) 


ois{q^ 


(4.166) 


Of course, we have not yet considered the corrections to the one-loop input. We could have built 
up these corrections to the input at the same time as we built up the evolution. However, it is 
easier to simply examine the change of the whole of (4.166) under a change in starting scale, and 
choose the input necessary to make the expression insensitive to such changes. In order to do this 
it is simplest to rewrite (4.166) in the form. 


77 + 

^ RSC,0 


{N,Q^) = 


a 


■{Ql 


a 


■{Q‘^ 


^ 01^1 + 


(AT) 


Fl{N) - 


36 - SAT/ 
27 


F^{N) 


exp 



"-(Qg) fi^(a,(g2)/iV) 

^(Q 2 ) a 2 (g 2 ) 


das{q^) 


+ higher order in N, 


(4.167) 

where the second term contains those parts of the one-loop input which are higher order in N. The 
factor {Fl{N) — F 2 {N)) must be multiplied by exp[ln(gQ/ylii)r°^(as(gQ)/A^)] in order 

to absorb changes in the evolution term under a change in starting scale Qq ^ (1 + 6)Qq. This 
does not absorb the whole of the change of the evolution term, but all other changes required in the 
input are subleading to this. Similarly the change in the input of the “higher-order-in-A^” terms 
needed to absorb the change in the evolution is entirely subleading in as(go) to the one-loop input. 
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or to {Fl{N) — (36 exp[ln{QQ/ALL)T^j^j^{as{Qo)/N)]. Hence, the full LO expression 

for F~^{N, Q^) is 


^RSC,oi^^Q ) ~ 


Fl{N) - 


36 - 8Nf 


27 


+F^’^'^{N) 


as{Q‘^ 


2 r“-‘.+(jv) 


exp 


F2(iV)J iexpMQl/ALL)TlLMQl)/N)] - 1) 

rl^(as{q^)/N) 


UasiQ^) 




-dasiq 


(4.168) 


We can solve for the corrections which are proportional to e (N) in exactly the same manner 
as above. Multiplying (4.151) by p~{N), instead of p^{N) leads to 


dF^^’-{N,Q^) 


= as (Q2)r0’*’- (iV)F"^-- (iV, Q") + 


dlnQ 2 


p-(iV). 


cl, -I 


^^z^\rl^{as{Q^)/N) 


\ 


Tl^{as{Q^)/N) ^ 2 , 

yas{Q^)rl^{as{Q^)/N) -(^^^^jas{Q^)Tl^{as{Q^)/N) j 


fI\n,q^) 

F^\N,Q^) 


(4.169) 

Again considering the contraction of the matrix between the two vectors in the last term we obtain 


p-{N) 


Tl^{as{Q^)lN) 

a,(g2)rija,(Q2)/jv) _ 


2 ^)rija,(g2)/Ar) 




^-^^\as{Q^)T\^{as{Q^)/N) 
_ f36-8Nf\, . __ N 


27 


27 


e+(lV) 


(as(g^)r^^(Q; 5 (g^)/iV) - —Tl^{as{Q‘^)/N)) + higher order, 

6 

(4.170) 


and 


p-(iV). 


TlMQ")/N) 

a,(g2)rL(as(g")/iV) 


36-87Vf 


^]Tl^{asm/N) 


27 


2^^'a,(g2)ri^(a,(g2)/iv) 


27 


e-(iV) 


(4.171) 


= 0 ' 




N 


Therefore we can write (4.169) in the form 

dF^h-(iv, g 2 ) 


d lng 2 


= a,(g2)rO’'’-(Ar)F"i’-(iV,g2) + 


36-8A^/ \ / «,(g^ ) 
27 J\as{Q^ 


2 . P3. + (JV) 




(asmTUasm/N) - -Tl^{as{Q^)/N))( F^iN) - 
+ higher order. 


36 - 8Nf 
27 


F2iN) 

/ 

(4.172) 
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Rather than solving this equation, it is easier to solve for the whole of the correction to F^’^’ {N, Q^) 
in one go. Using (4.170) and (4.171) we easily obtain 


MQ^)rlL{asm/N) - jrl^{as{Q^)/N))(^FL{N) - x 

//■_ /Arv r f|^(as(g2)/iV) 2\] 


2 U'^’ + CiV) 


exp[ln(Qo/^LL)r'^i(a*(Qo)/iV)] exp 
+ higher order. 


'as{Q^) 


-dusiq 


(4.173) 


The solution to this is relatively simple if we wish to keep only the LO parts. Letting 
^ we obtain, 

^c/3«,-g 2 ^ ^ g 2 ) 


exp[ln(Q2/^^^)rO Ja,(Q2)/^)] 


a 2 n\ f°’''+(iV) 




'o^siQ^) 


dois{q^ 




exp[ln(Q^/^LL)r^i(a,(Q^)/iV)] 


a 2 n\ f°’''-(iV) 


+ higher order. 


(4.174) 


Adding to the solution at one-loop, the whole of the LO part of F {N, Q'^) is 


_ f36-8Nf 


F'ysc,oi^^Qn= Fy^’-iN)- 


Fl{N) - 


F2{N) X 


(ayQl)Tl^iayQl)/N) - f L^ Ja,(Q2)/^)x x /«^(q2)x f«''.-(iv) 

V Tl^{ayQl)/N) ))\asm) 

36-8Nf\ fayQl)Tlj^{ayQl)/N) - fTl^{as{Ql)/N)\ f ^ ^ f36-8Nf\ ^ 

I - Tt.ia.mm - ) " 


exp[ln(Q2/A^^)rO^(a,(Q2)/^)](^_^ 

+ higher order. 


a 2X\ f“-'' + (iV) 


p(«o) fL(a,(g2)/jV) 
/aRQ 2 ) as(g^) 


da.siq'^ 


(4.175) 

where, in order to make the expression invariant under changes in starting scale, and also ensure 
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that for Qq =All we have F {N,Qq) = (N), we make the choice 


/ 36-8iVA f asiQDrl^iasiQD/N) - f hj \ 

V 27 A Tl^MQl)/N) J 

(^Fl(N) - (^ 36 -8^/ y^(jv)^(exp[lri(QgMz.L)r»L(«.(Qo)/A^)] - !)• 

(4.176) 


We have been able to obtain, without too much difficulty, the leading parts of F~^{N,Q‘^) 
and F~{N,Q‘^). We must now use these in order to obtain LORSC expressions for the structure 
functions. The way in which we have set up the calculation makes this very straightforward for the 
longitudinal structure function: we multiply F^g^ q{N,Q‘^) and F^g^ q{N,Q‘^) by (as((5^)/27r). 
However, we notice that all the parts of F^^^ ^{N, Q'^), except the one-loop expression, are sublead¬ 
ing to F^g^ ^{N, Q^): all the terms in the inputs in the former are a power of OsiQ^) higher for the 
same power of N. Thus, the only part of F^g^j which contributes to the LORSC expres¬ 

sion for Fl{N,Q^) is the one-loop part. Adding this to (4.168), and multiplying by (as((5^)/27r), 
we obtain 




27r 


as{Q‘^ 


exp 




L4a,(Q2) 


as{q^ 


-das{q 


F 




f’+(iV)+(Fz.(iV)- 
2 ' (N) 


36 - 8Nf 
27 


OsiQ'^, 


F,iN)]{exp[HQ^,/ALL)nLMQt)m - 1 ) 


(4.177) 


It is slightly more involved to find the LORSC expressions for F 2 {N,Q‘^). In this case we 
wish to find the LORSC expressions for the derivative {dF 2 {N,Q‘^)/d InQ^) and for the input 
F 2 {N,Qq). We consider the former first. Using the form of e+(A^) and e~{N) in (4.148) it is 
clear that, besides for the one-loop contributions, the LORSC expression will come from (27/(36 — 
8Nf)) X {dF 2 {N,Q‘^)/d InQ^) and from N/6 x {dF^{N,Q‘^)/d InQ^). Explicitly we obtain 


dF2{N,Q^ 

d lnQ2 


= OisiQl) 


RSC,0 


e-(lV)rO’^’-(iV)F°’'’-(iV)(^^ 


2 . 


+ e+(iV)r°’'’+(iV)F^’‘’^(iV) -r2’'/(iV)( FL(Ar) - 


^ 0 ,/, 


1,0 


36 - 8Nf 
27 


F2iN) 




exp 








OisiQ'^ 


(4.178) 
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where we have used the relationship 


^{Tl^{a,{Q^)/N) - Tl^{as{Q^)/N)) = /N) - T\^{as{Q^)/N)), (4.179) 


and T 2 ’^j^{N) ^Fl{N) — y^^-^^jF 2 {N)j is the part of the input which is common to both the 
one-loop and the leading-ln(l/x) input, and is subtracted from the former in order to avoid double 
counting. Similarly the input for F 2 {N, Q'^) is given by the relatively simple form 


F2,rscAN,QI) = F2{N) 

+ " (36 {eMHQl/^LL)rlLMQl)/N)] - 1) 

+ HQl/ALL)as{Ql)(^e+{N)T^’^’+{N)FlF+^N,Ql) + e-{N)T^’^’-iN)Fl’^’-{N,Ql^ 

-r^:°(iV)(Fz.(iV)- ~ 27 ^^ )^ 2 (JV)^ 

(4.180) 

The third term is the renormalization-scheme-invariant order-as((5Q) input, which must compen¬ 
sate for changes in the one-loop evolution under changes in Qq . Again we explicitly extract a term 
oc ln((5o/^)r2 ’L(A") ^-^l(A^) — -^ 2 (A^)^ in order to avoid double counting. The fact that 

this term can be thought of as appearing from two different sources leads us to choose both our 
unknown scale constants in the input equal to the same value We note that our choice of 

inputs not only ensures QQ-invariance up to higher orders, but are also such that our expressions for 
Fl,rsc,o{F',Qq) and F2^rsc,o{F',Qq) reduce to the nonperturbative inputs if Qq = All. Having 
obtained our LORSC expressions for {dF 2 {N,Q‘^)/d InQ^) and F 2 {N,Qq), then as already argued, 
in order to obtain our expression for F 2 {N,Q^) we integrate {dF 2 {N,Q‘^)/d lyiQ‘^)Rsc,o from Qq 
to and add to the input T2 ,rsc,o(A^) Qo)- 


Thus, we have our complete leading-order, including leading-ln(l/x) terms, renormalization- 
scheme-consistent expressions for the structure functions. These are significantly different from 
both the one-loop expressions and the leading-ln(l/a:) expressions, although they do reduce to them 
in the appropriate limits; i.e. to the one-loop expressions for very large x, and to the leading-ln(l/x) 
expressions in the limit of very small x and near the boundary of evolution, ~ Qq. Indeed, 
once we include the 0 (q!s((5o)) inputs for F 2 {N, Q^) in the definition of the leading-order inputs in 
the loop-expansion (as we should), each of the terms in the inputs and evolution terms in (4.177), 
(4.178) and (4.180) contains a part which appears in both the LO expression in the loop expansion 
and the LO expression in the leading ln(l/x) expansion. Our full LORSC expressions are obtained 
rather more easily than above by simply letting each of the input and evolution terms become the 
combination of the terms in the two expansion schemes. In a sense this result is obvious, but it is 
necessary to verify this by deriving the expressions as above. 
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Let us comment on the form of our final LORSC expressions. We note that as for the leading- 
ln(l/x) expansion, the LORSC expansion still leads to predictions for all the small-x inputs: 
predictions of each in terms of the nonperturbative inputs (which we imagine should be quite fiat) 
and the nonperturbative scale All, and also stronger predictions for the relationships between the 
inputs (although the scale Qq at which they should be chosen is not determined). 

Finally, we notice that each of the terms appearing in our expressions is manifestly renormaliza¬ 
tion scheme invariant, and it is clear that no terms are subleading in as to any other terms, either in 
the input or in the evolution. If we had simply solved the renormalization group equations using the 
whole of the anomalous dimension matrix (4.142) then we would have obtained many terms which 
do not appear in our full leading-order expressions (4.177), (4.178) and (4.180). These would still 
be renormalization scheme independent (and trivially factorization scheme independent of course), 
since our input anomalous dimensions are renormalization scheme independent. However, all of 
these extra terms would be of the same form as terms which must be renormalization scheme depen¬ 
dent in order to absorb the changes of the leading-order expressions under a change of the coupling 
coming from a renormalization scheme change, i.e. as{Q‘^) as{Q‘^) + 0(ag(Q^)) (e.g. we saw in 

4.2 that the subleading-in-ln(l/x) evolution <1>^ (Q^, Qq) has a manifestly renormalization-scheme- 
independent part depending on T^Lio^siQ"^)/]^) as well as a renormalization-scheme-dependent 
part depending on T\]^{as{Q‘^)/N)). Thus, these terms should be dropped, and (4.177), (4.178) 
and (4.180) are the correct expressions for the structure functions to be used with the one-loop 
coupling constant. Finally, as already stated, when considering the expressions for the structure 
functions expanded about negative integer values of N we will simply use the one-loop expres¬ 
sions, due both to lack of knowledge of the explicit resummed anomalous dimensions which would 
occur, and also because it almost certainly makes practically no difference to do this as far as 
phenomenology is concerned. 

Now that we have our desired expressions, (4.177), (4.178) and (4.180), we can see how they 
compare to the data. The first step towards this is to consider the solution in x space. After briefly 
doing this we will consider detailed fits to the data. 


5. X—Space Solutions. 

We shall now discuss how we use the expressions (4.177), (4.178) and (4.180) in order to obtain 
our expressions for the x-space structure functions and ultimately compare with data. The data 
on J- 2 {x,Q^) exist over a range of from 0.2GeV^ to SOOOGeV^, so it is clear that we will 
need expressions for the structure functions which cross quark thresholds. Thus, before presenting 
details of the x-space solutions we must consider how we treat these quark thresholds. 

In practice we will impose a lower cut on the data of = 2GeV^, except for the HERA data 
where we choose = 1.5GeV^ simply in order not to lose some of the very low x data. The 
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threshold of heavy quark production is = 4m^, where W is the invariant mass of the hadronic 

system created from the struck proton (or neutron) and in the limit of zero proton (and/or neutron) 
mass is given by = Q^(x~^ — 1). mjj is the mass of the heavy quark. Thus, we clearly work 
in the limit where the up, down and strange quarks are effectively massless. However, we cross the 
h-quark threshold of ~ 20GeV^, and at the lower end of our range are in the region of the c- 
quark threshold of ~ 2GeV^. The correct treatment of these heavy quark thresholds is less well 
established than the treatment of effectively massless quarks, and is certainly more complicated 
(especially when considering leading-ln(l/a;) terms). Hence, we will delay a more sophisticated 
treatment to a future article, and in this paper use a relatively simple treatment. 

We use the prescription for treating heavy quarks outlined in [40]. This involves treating all the 
quarks as massless, but only allowing the heavy quarks to become active above the simple threshold 
= m^. Hence, the value of Nf appearing explicitly in any expressions changes discontinuously 
at this threshold. The running coupling constant, which itself depends on Nf, is dehned to be 
continuous at the thresholds. It is determined by the relationship 

as,n(Q^) = as,n+i(Q^) ^ , (5.1) 

where the central as{Q‘^) with Nf = 4 is dehned by 

/ 2 \ ^ 27r . , 

^ ^ (33-2.(fV, = 4))ln(QVA|cD,4)’ 

This leads to a kink in as{Q‘^) at the thresholds, with the coupling used below a threshold being 
larger than the continuation downwards of the coupling used above the threshold. 

This complete prescription for treatment of the heavy quarks is consistent with the decoupling 
theorem, as it is guaranteed to provide the correct expressions far above or below any threshold, 
as discussed in [40]: the increase in the coupling below a threshold compensating for the absence 
of virtual heavy quarks in calculations below this threshold. It is clearly going to be rather un¬ 
satisfactory in the region of the threshold, with heavy quark structure functions having an abrupt 
threshold in rather than the physically correct smooth threshold in IT^. Work to rectify this 
is in progress, and will certainly involve the use of the heavy quark coefficient functions at leading 
order in ln(l/x) already calculated [13] [30]. We note, however, that the treatment of quark thresh¬ 
olds in this paper is no less rigorous than in most of the other calculations of structure functions 
currently performed, and is perhaps better than some. 

The number of active quark flavours, Nf, appears in (4.177), (4.178) and (4.180), in a number of 
places other than the implicit dependence in It appears in the one-loop physical anomalous 

dimension eigenvalues (^N), in the eigenvector factors e^^~\N), and in the definition of 
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Fl{N)', i.e. C£ ^ i{Fl) oc Nf and hence, from the definitions (4.36) and (4.37) which define the 
nonperturbative inputs for the structure functions, we have 

h{N) - Cl,^i{N)F 2 {N) = Fl{N) - 3(/^^) ^2(A^) oc Nf. (5.3) 

The first two forms of dependence have very little impact in practice, the -dependence being 
relatively weak. The last dependence, the proportionality of what is in practice almost all of Fl{N) 
to Nf, has a large impact. It means that the expressions for the longitudinal structure function 
and for the InQ^ derivative of F 2 {N, Q^) are discontinuous at and at = ml- Also the 

definition of the input F2{N, Qq) is sensitive to where Qq is chosen with respect to the thresholds. 

This dependence of Fl{N) on Nf is obviously rather unsatisfactory. When treating the heavy 
quark thresholds more rigorously we should subtract out of the full structure functions those con¬ 
tributions where there are particles with heavy quarks in the final state. Doing this, we could 
define intrinsic light quark structure functions where the final state particles do not contain heavy 
quarks, and which have a fixed number of active flavours. Our definition of the nonperturbative 
inputs could come from these light quark structure functions alone. There would then also be heavy 
quark structure functions which could be calculated separately from and entirely in terms of the 
light quark structure functions. However, using our somewhat simplistic prescription for turning on 
the heavy quark contributions, even if we were to define heavy quark structure functions completely 
separately from the light ones, the expressions for Fl ^h{N,Q‘^) and {dF 2 ^h{N,Q‘^)/ dhi{Q‘^)) would 
be discontinuous at the thresholds. Thus, in this paper we feel that the treatment of quark thresh¬ 
olds makes it more appropriate to define the whole structure function in terms of a discontinuous 
number of active quark flavours. The contribution to this complete structure function from the 
production of heavy quarks can then be extracted straightforwardly. Hence, we treat the problem 
of quark thresholds as described above. We will discuss the effect of the discontinuities at the 
thresholds on fits to data in detail in the next section. We note, however, that the expressions for 
the whole structure functions in terms of the singlet and nonsinglet components in (3.1) depend on 
the average value of the squared charge of the active quarks, and this will also change discontinu- 
ously at the thresholds. This means that overall the effects of heavy quark thresholds on a best fit 
to F2{x,Q‘^) data are in practice rather small. The effects of the 6-quark contribute only like the 
charge squared, i.e. 1/9, and hence have little impact. The effects due to the c-quark threshold 
are proportionally much larger, but because they are at ~ 2GeV^ they only affect a very small 
proportion of the complete set of data. 

Having defined our treatment of quark thresholds we can now discuss the form of the x- 
space solutions for the structure functions. As already discussed it is the perturbative part of the 
moment-space structure functions for which we can produce well-ordered, RSC expressions and 
the nonperturbative parts of the expressions for these structure functions, F2{N) and Fl{N), will 
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be nonanalytic, complicated functions of N. Our complete moment-space expression for a general 
structure function (or derivative of a structure function) will be 

Fi{N,Q^) = Pi,2{N,as)F2{N) + Pi,L{N,as)FL{N), (5.4) 

where the q;s)’s are the calculable perturbative components of the complete expressions 

which can be expanded as power series in both and N. Hence, the Pij{N,asys are an example 
of the physically relevant perturbative functions we discussed in §4.4. Taking the inverse Mellin 
transformation of (5.4) back to x-space we obtain 

xFi{x,Q^) = Vi,2{x,as) ®F2 {x) + T’i,L(x, a^,) ig) J^l(x), (5.5) 


where g) denotes the convolution of two quantities, i.e. 


/■^ d z 

A{x) B{x) = X / — A{xjz)B{z). 

/ /y. Z 


(5.6) 


Making an ordered calculation of the structure function is then equivalent to making an ordered 
calculation of the Vij{x, as)’s. As proved in §4.4, this is equivalent to calculating these perturbative 
parts of the expressions for structure functions as RSC expansions about each of their nonanalyt- 
icities in N, and then taking the inverse Mellin transformation of each of these expansions by 
integrating around curves enclosing each nonanalyticity. 

Our moment-space expressions (4.177), (4.178) and (4.180), are part of the full expressions for 
the structure functions in the form in (5.4). They contain the complete nonanalytic expressions for 
F 2 {N) and Fl{N), but the factors multiplying these are LORSC expansions of the perturbative 
parts of the expressions about N = 0. Thus, in a sense the expressions (4.177), (4.178) and (4.180) 
are misleading. They need to be interpreted correctly. We cannot simply take an inverse Mellin 
transformation of these expressions because the result would depend on all the singularities in 
the nonperturbative functions, whereas the perturbative factors in these expressions are defined as 
power series expansions in N about = 0 with radii of convergence of only unity. The expressions 
do not have a direct meaning for general values of N. Instead, we must interpret them as providing 
us with the perturbative factors multiplying F 2 {N) and Fl{N) which are LORSC when expanded 
about N = 0. 

Hence, we use (4.177), (4.178) and (4.180) to obtain the x-space structure functions as follows. 
We can factor out the perturbative part of each of these expressions which is proportional to either 
Fl{N) or F 2 {N), and take the inverse Mellin transformation of this perturbative function by 
integrating around a contour encircling = 0. This provides us with the leading-in-ln(l/x), at 
lowest-order-in-x parts of the perturbative V{x, as)’s as well as all the lowest-order-in-cTs parts at 
lowest-order-in-x. In order to obtain the leading-order-in-ln(l/x)-and-Q!s, at lowest-order-in-x, 
part of the perturbatively calculated structure functions we then convolute with the whole of the 
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nonperturbative Fl{x) and F 2 {x) which are obtained by a complete inverse Mellin transformation of 
Fl{N) and F 2 {N), including contributions from singularities at many values of N. In principle, the 
full LORSC x-space structure functions are calculated by repeating this procedure for the LORSC 
moment-space structure functions where the perturbative parts are expanded about N equal to all 
negative integers. This would systematically include all leading-order“in-ln(l/x)-and-a!s terms at 
progressively higher orders in x. As already explained, in practice one only really needs include the 
one-loop, or leading-in-a^ part of the full LORSC perturbative parts for all negative integers, and 
also only need to work to a relatively small negative integer unless one is concerned with values of 
x very close to one. Thus, it is this type of calculation we aim to make in order to compare our 
expressions with data. 

In practice the calculation of the structure functions is performed in a rather different manner. 
The structure functions are calculated using a modification of the computer program that is used 
by Martin, Roberts and Stirling in their global fits to structure function data. This works in 
terms of parton densities and also directly in x-space. Input parton densities are specified at 
some scale Qq and the evolution is calculated on a grid in x and Q^. This evolution is 
obtained by integrating up the complete renormalization group equations involving the complete 
parton distributions and full specified splitting functions. The structure functions are calculated 
by numerically performing the convolutions of the resulting Q^-dependent parton distributions 
with coefficient functions. Thus, in fact we obtain the factorization-scheme-independent LORSC 
structure functions by working in terms of parton densities obtained from the full solution of the 
renormalization group equations and choosing coefficient functions and splitting functions (i.e. an 
effective factorization scheme) which will reproduce the correct expressions as closely as possible. 

In order to do this we choose to work in a DIS-type scheme for simplicity, and take the 
longitudinal coefficient functions to be just the one-loop expressions. By breaking the expressions 
for the LORSC structure functions in terms of the structure functions themselves into expressions 
involving partons with the above choice of coefficient functions, and also using the power series 
expansions of the one-loop anomalous dimensions and eigenvalues in terms of N it is possible 
to choose effective splitting functions which when used in the computer program will reproduce 
the correct expressions for {dF 2 {x,Q'^)/dln{Q‘^)), .F 2 (x,( 5 o) and Fl{x,Q'^) to very good accuracy. 
We note that this choice has to be different when calculating F 2 {x,Q‘^) from when calculating 
Fl{x, Q^). Due to the method of calculation used by the MRS program, it is extremely difficult to 
obtain the exact behaviour one desires of the structure functions, i.e. a given evolution away from 
a particular input at some fixed Qq, with any choice of splitting functions. However, we aim to get 
as close as possible to this exact behaviour with relatively simple choices. 

We will briefly describe our choice of effective splitting functions as follows. At first order in 
as we choose the normal parton splitting functions. We then add corrections to these in order to 
reproduce our desired results. Denoting these corrections by A'yab{N,as{Q‘^)), for a given parton 
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to parton splitting function, and expressing these in moment space and in terms of previously 
discussed quantities for simplicity we get for Q^) 


A7/3(iV,a.(Q')) = 0 
A7//(iV,a,(Q2)) = o 

A7,/(iV,a,(Q2)) = lTl^{as{Q^)/N). 

This form is not too difficult to understand by looking at e.g. (4.177); the only leading~ln(l/x) 
enhancement of the one-loop expression comes from the evolution in the longitudinal sector which 
is directly related to the enhanced gluon evolution and hence corrections to the gluon anomalous 
dimensions. For the corrections are a little more involved: 


A-ifg{N, as{Q‘^)) 


= asmCl,iiN)TUas{Ql)/N) 




8A/ 

27 


“ Q ) + 


rk(«s(g^)/A) 


A^gg{N,a,iQA) =TIlMQA/N) - lAjfg{N,asiQA) 
A^gf{N,a,{QA) = ItIMQA/N) - ^A^ff{N,as{QA)- 


(5.8) 


The corrections to the quark anomalous dimensions are functions of rather than , ex¬ 
cept for the single power of as{QAi because T\^{as/N) appears only in the input terms in 
(4.178). In the small~x, or small-A limit we have the simplifications that A'-)fg{N,as{QA) 
as{Q‘^){2Nf/‘i)T\^{as{Ql)/N) and A'yff{N,as{QA) ^ |A 7 /g(A,^^(Q^)). Thus, in this limit 
these corrections to the one-loop quark anomalous dimensions take on the standard form of 
as{Q‘^)'yjgAs{Qo)/^) as{Q^)AfAs{Qo)/^) minus their one-loop components. The pref¬ 
actors multiplying these terms, which depend on the one-loop longitudinal coefficient functions, 
ensure that the leading-ln(l/x) enhancement of the rate of growth of T 2 {x, Q'^) is directly coupled 
to the longitudinal structure, not to the gluon. This delays the enhancement to slightly smaller 
values of x, and reduces it a little. The contributions to the gluon anomalous dimensions beyond 
LO in ln(l/x) are present to counter the exponentiation of the corrections to the quark anomalous 
dimensions, i.e. this is similar to choosing '^^^{as/N) equal to —^AfgAs/Af) in (4.78)-(4.80) in 
order to ensure that the exponentiation of 7 j^(q;s/A) has no effect in these expressions. The effect 
of these terms is in practice rather small. 

We may now perform checks to see if the particular choice above does indeed lead to an 
accurate representation of the correct expressions. In order to do this we obtain solutions for 
the structure functions which are the same as those obtained by the MRS program for our given 
choice of splitting functions, by analytically performing the same calculation as performed by this 
program, but in moment space, where the calculation is possible. The resulting expressions can 
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be compared with the correct LORSC moment-space expressions in the region of relatively small 
X {x ^ 0.1) by expanding both expressions in powers of N, and an accurate approximation of 
the differences between the two found.These differences between the expressions can then be 
estimated by transforming back to x-space. With the choice of the splitting functions above the 
differences between the expressions are very small. 

One may also take a more direct approach. For a simple, but physically reasonable, choice of 
input parton densities, and hence input structure functions, one can use the types of techniques 
outlined in [22] to begin with the exact expressions (4.177), (4.178) and (4.180) in moment space 
and then transform to x-space obtaining analytic solutions for the structure functions. These 
expressions are very complicated, but are almost exact for very low values of x. They begin to 
become approximate for larger values of x, particularly at high due to the effects of truncating 
power series in N. The range of x for which the expressions are reliable depends on the severity 
of the truncation but it is not too difficult to have accurate expressions for x 0.1. These enable 
one to compare the results obtained from the inverse Mellin transformation of our correct moment- 
space expressions with the calculations performed using the computer program in x-space. For 
X ^ 0.1 our choice of the anomalous dimensions leads to agreement with the analytic expressions to 
an accuracy much better than the errors on the data in any appropriate range of parameter space. 
There is no reason to believe that input distributions which are very similar in form, but which need 
to be described by more complicated expressions, should lead to deviations much different from our 

The ease of obtaining an accurate expression at small x is aided immensely by the fact that the 
expansions about = 0 of the one-loop parton anomalous dimensions, and subsequently their eigenvalues, 
eigenvectors and projection operators have coefficients which stay roughly constant as the power of N 
increases, or even where the lowest powers have the largest coefficients. For example, the expansion of 
^0’*"+(7V) about = 0 goes like {&/N — 5.64 + 0.24Af — 2.S7N^ + 3.38A^^ + •••)• This enables one to 
truncate these series after a very small number of terms when working at small x, or conversely work 
to quite large x without keeping too many terms in the series. It is also the reason why the “double 
asymptotic scaling” formula, which effectively truncates the above series after the second term, is such a 
good approximation to the full one-loop solution (with flat parton inputs) up to quite high values of x. In 
the MS renormalization and factorization scheme, or the DIS factorization scheme, the two-loop anomalous 
dimensions, and in the former case coefficient functions, do not exhibit nearly such a nice behaviour, e.g. 
the MS scheme jlg{N) expanded about N = 0 goes like (26.7/N — 57.7 + 104A'' — 176A^^ + • • ■), and 
the expansion of the coefficient function is even worse, ^ i(N) ex (1.33 — 6.16Af + 9.2N^ — 12N^ + •••). 
Hence, truncations of the two-loop anomalous dimensions and coefficient functions expanded as series about 
N = 0 give a far worse approximation to the effects of the complete anomalous dimensions and coefficient 
functions at small x than in the one-loop case. This has been explicitly checked, and indeed, evolution and 
convolution of a given set of parton distributions using the series expansions for the two-loop anomalous 
dimensions and coefficient functions truncated at lowest order in N lead to structure functions up to a 
factor of two larger than those obtained from calculations performed using the full anomalous dimensions 
and coefficient functions. 
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particular choice. Hence, using this direct test and the more general, if slightly more approximate 
check above, we are confident that our choices do lead to the correct x-space expressions for x ^ 0.1 
up to very small corrections. 

At higher x the check is simpler, in fact almost trivial. Simply setting to zero all terms in 
the MRS program other than those coming from one loop, we observe that the expressions for the 
structure functions for x > 0.1, particularly at high are very close to the one-loop expressions 
alone (becoming essentially identical for x > 0.3). This result agrees with calculations taking the 
inverse Mellin transformation of the LORSC moment-space expressions. Hence, the direct x-space 
calculation must be the same as our desired expression in this relatively high x range. Therefore, the 
correct choice of effective splitting functions leads to the MRS program producing a very accurate 
approximation to our correct LORSC expressions for the structure functions over the complete 
range of parameter space. 

We note that with our choice of splitting functions and definitions of parton densities then 
momentum is not conserved by the evolution: in the best fit to J- 2 {x,Q‘^) discussed in the next 
section the total momentum carried by the partons at = 2GeV^ is 87% and at = SOOOGeV^ 
it is about 94%. Hence, the amount of momentum violation is at the level of a few percent. We 
have already defended this violation of momentum conservation in subsection 4.3. We now also 
point out that starting with our definition of the partons we could define new parton densities 
and splitting functions by defining non-zero C(^f\as/N) and C(}f\as/N) beyond one-loop, and 
use the transformation rules in §2 to keep the structure functions unchanged. If these coefficient 
functions had negative coefficients then, compared to our prescription above, the low parton 
distributions would need to be larger, and would hence carry more momentum. As increased 
the effect of these new coefficient functions would decrease, and the extra amount of momentum 
carried by the new parton distributions compared to the original ones would decrease. Thus, the 
effect of such a redefinition of parton distributions would be to increase the amount of momentum 
carried by the partons at low and also to slow the growth of momentum with (or even to 
turn it into a fall). With a judicious redefinition of coefficient functions of this sort (with perhaps 
some dependence on UsiQ^) as well as as{Q‘^)) it should clearly be possible to find an effective 
factorization scheme where the structure functions are identical to those given by our prescription, 
but where the momentum violation will be extremely small. As already mentioned, one may think of 
this as a “physical scheme”. We have not seriously investigated this redefinition of parton densities 
in any quantitative manner since it will not affect any physical quantities. 

Now that it has been determined that our choice of splitting functions is correct, one can 
vary the input parameters in the standard way to obtain the best fit. In practice the starting 

In fact the agreement between the result obtained from the computer program and the “correct” result 
is certainly better than the agreement between NLO calculations performed using different programs and 
different prescriptions for truncating the solution at NLO [34]. 
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scale for the numerical evolution is chosen to be Qq = so that the charm contribution to the 
structure functions is guaranteed to turn on at this scale. Evolving upwards, the bottom quark 
contribution turns on at Thus, the parton distributions and, via convolutions with 

coefficient functions, the structure functions are input at Qo- However, this does not mean that 
the true starting scale, Qg, has to be the same as Qo- This true starting scale is the scale at which 
some of the couplings in the anomalous dimensions in our full LORSC expressions have to be fixed. 
When performing the numerical integration upwards form a lower scale we simply interpret the 
values of {dT 2 {x^Q‘^)/ dlnQ"^) and T 2 {x,Q‘^) at the fixed scale Qq as our input values 

for these structure functions. We then demand that at Qq the relationship between these quantities 
is the same (up to an accuracy much better than the errors on the data) as that demanded by the 
expressions (4.177), (4.178) and (4.180). We also demand that the value of each of these quantities 
at Qq is compatible with them being written as the convolution of our completely determined, up 
to the scale All, perturbative parts and our nonperturbative soft inputs. In practice we demand 
that these nonperturbative inputs must be flat at small x, i.e. they must be described well by a 
function of the form 

Ti{x) = Fi{l - x)’^*(l + eiX^-^ F^ix). (5.9) 

It is by satisfying these requirements on the inputs, while simultaneously choosing Qq as the fixed 
scale appearing in the anomalous dimension or splitting function dependent input terms in our 
expressions, which enables us to identify Qq as the true starting scale. The evolution which takes 
place upwards from Qq{= rn^) to Qq would be identical if it were to take place backwards from Qq 
to Qq. The former manner of performing the evolution is simply more convenient in practice. 

We should perhaps comment more on this starting scale Qq, and on its relationship to All- 
All has been defined as the particular choice of Qq for which the inputs for the structure functions 
reduce to the purely nonperturbative inputs. Hence, we imagine it is a scale typical of nonperturba¬ 
tive physics, i.e. All ^ IGeV^. If we choose Qq different from (in practice higher than) All, then 
All is not precisely the scale at which the full expressions for the structure functions reduce to 
the nonperturbative values, although it will remain fairly close to this. Hence, All is simply some 
phenomenological parameter which we will fine tune in order to get the best fit, though it does 
maintain its interpretation as a scale typical of the transition between perturbative and nonper¬ 
turbative physics, and it would be surprizing if it were much greater than say IGeV^. As already 
mentioned, the choice of Qq is undetermined. However, it should be such that it does make our full 

As already mentioned, a more sophisticated treatment of the heavy quarks would involve calculating 
heavy quark structure functions separately from the light quark structure functions. More sophisticated 
treatments within the usual loop expansion, e.g [41], have the thresholds for heavy quark production built 
into the coefficient functions, and it is not important to begin numerical evolution at, or below any particular 
value of in this case. 
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expressions relatively insensitive to its value. (We also reduce this sensitivity by letting All be a 
free parameter for each choice of Qq, though it should not vary a great deal.) This insensitivity to 
Qq imposes two restrictions on its value. Firstly, Qq should not be too small, otherwise the value 
of asiQo) be very sensitive to our choice of Qq, making it difficult to make our expressions Qq 
insensitive. However, we cannot simply choose very high Qq. This is because terms beyond LO in 
the full RSC expressions will contain terms of higher order in \n.{QQ/A ll)-, albeit along with higher 
orders in asiQ"^)- Hence, we do not wish to choose Qq to be a great deal larger than the scale All- 
Thus, Qq is probably some intermediate value. It will, of course be determined in practice by the 
quality of the fits to data using the expressions (4.177), (4.180) and (4.178) evaluated at particular 
values of Qq. 


Now that we know precisely how we will perform our calculations of structure functions we 
can make some comment on the general form the structure functions have to take. The expressions 
for Tl{x,Q'^) and J^ 2 ix,Q'^) depend on r^^(as/N') and T^Liois/^)- Both of these series have 
coefficients which behave like n“^/^(121n2)/7r for very large n. This leads to a cut in the N- 
plane at N = Q;s41n2 in both cases, and to the structure functions having asymptotic behaviour 
T{x,Q‘^) ~ (lnx)“^/^x“^''^^“'' as X ^ 0, where as = as(Qo) Q^ > Qq. 

However, this is only for the strict asymptotic limit x ^ 0. It was convincingly demonstrated 
in [20] (and discussed from a different point of view in [22]) that one only need keep a finite number 
of terms in the leading-ln(l/x) series if working at finite ln(l/x). In practice, if one works down 
to X ~ 10“® and the nonperturbative inputs behave roughly like (1 — x)^, then keeping terms up 
to 10th order in the series for r°^(Q;<j/A^) and T^Liois/N) is more than sufficient (truncating at 8th 
order leads to only tiny errors). The series up to this order have the explicit form 
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It is clear that the low-order terms in both these series have coefficients which generally grow 
far less quickly than the asymptotic relationship a„_|_i = 41n2ajj; some fall, and in the case of 
are even zero. Thus, in the range of parameter space we are considering we will have 
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rather less steep behaviour than the asymptotic limit oi ~ (Inx) 4in2as^ Using the 

techniques in [22] to derive analytic expressions for the structure functions in the small-x limit, it 
becomes clear that the extreme smallness of the first few coefficients in the series for r^^(as/A^) 
leads to the small-x form of J-l{x,Q‘^) being driven largely by the first term, with the 4th“Order 
term having a reasonable effect at the smallest x values. Similarly, the small-x form of J- 2 {x,Q‘^) 
is determined largely by the series T\^{as/N) and the first term in again with the 

4th"Order term in r°^(Q;<j/iV) playing some role at the smallest x values. 

One can be a little more quantitative. Choosing = 25GeV^ and making a guess that 
All = O.SGeV^ we obtain All) ~ 3.4. We also choose nonperturbative inputs which behave 

like Fi{l — x)^ for i = L, 2. As argued in [22], for values of x ^ 0.01 it is a good approximation to 
replace the (1 — x)® behaviour with the behaviour 0(0.1 — x)?^ Doing this we can derive simple 
expressions for the form of the inputs for x ^ 0.01 which will be accurate up to a few percent error. 
Using the series (5.10) and (5.11) and the expressions (4.177) and (4.180) we can take the inverse 
Mellin transformations (using the techniques of [22]), obtaining 


J^l{x,QI) ~ as{Ql)/2'K 
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where ^ = ln(0.1/x), and as{QQ) ~ 0.2 if Aqcd,4 ~ lOOMeV. Putting in values of Fl = 2.5 and 
F 2 = 1, choices which, as with the (1 —x)® behaviour, are roughly compatible with the high x data, 
we have a rough estimate of the form of the input structure functions at Qq = 25GeV^. We see 
that the coefficients in the series in agiQ^)^ multiplying (Fl + 4/27+2) grow a little more quickly 
for + 2 (x,gg) than for Fl{x,Qq). However, in the former case this contribution which rises with 
falling X is accompanied by the flat + 2 , whereas in the latter it is accompanied only by the nearly 
insignificant —4Q;s(gg)/(547r)F2. Thus, we find that for x w 10“"^ Fl{x,25) behaves approximately 
like 0.1x“°'^ but that +2(x,25) is slightly flatter, behaving approximately like 0.6x“*^'^®. These 
powers of x are clearly somewhat less steep than the asymptotic x~^'^. Gomparison of the estimate 
of +' 2 (x, 25) with the data in [1] and [2] shows a very reasonable qualitative agreement. Of course. 


The argument of the 0-function depends on the power of (1 — x). Higher powers would require the 
step to occur at lower x and vice versa. 
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there is no data for ^i(x,25) for values of x anything like this low. Being rather more general, 
we hnd that for any Qq between lOGeV^ and lOOGeV^, and with any sensible choice of All (e.g 
0.2GeV^ > All > 2GeV^), then both F2{x,Qq) and Fl{x,Qq) behave roughly like for 

0.01 > X > 0.00001. Hence, this type of behaviour can be taken as a prediction of the theory. 

One can also use the same techniques to make an estimate of the structure functions at values 
of away from 25GeV^. The expressions involved are relatively straightforward to derive but 
become rather less concise than (5.12) and (5.13), and we will not write them explicitly here. 
With the particular inputs above they do lead to good qualitative agreement with the data below 
X = 0.01. Therefore we have every reason to feel encouraged by our results and believe that we 
really are proceeding correctly. However, the real test of our approach will be a complete global fit 
to the available data for F 2 {x,Q‘^) using the rather more accurate calculations, especially at large 
X, of the MRS program. We will therefore discuss these detailed fits next. 


6. Fits to the Data and Predictions. 

As explained in detail in the previous section we indirectly use the expressions (4.177)-(4.180) 
to calculate the lowest-order-in-x, RSC part of the perturbative part of the solution for the x-space 
singlet structure functions, and convolute this with the nonperturbative parts. The-higher-order- 
in-x parts of the perturbative solution for the singlet structure function, along with the nonsinglet 
structure functions are calculated using the normal one-loop prescription. Phenomenologically this 
procedure is practically identical to using the full LORSC solution for the structure functions. Let 
us now explain how we use the procedure to obtain fits to structure function data. 

Once we have the general LORSC expression then by combining the singlet and nonsinglet 
components and varying all the free parameters {All, the soft inputs for Fl{x, Q'^) and F^ix, Q'^) 
and the soft nonsinglet inputs), we obtain the best fit for the available T 2 structure function data 
using a particular starting scale Qo- We note that at this starting scale the input .F^(x,(5o) and 
the evolution {dT 2 {x,Q'^)j are forced by (4.178) and (4.180) to be trivially related at 
small X. This is not the case when working at fixed order in and indeed, is a new feature of 
the approach to the small~x structure functions in this paper, and in particular of the manner of 
determining the inputs that is used. Choosing the renormalization scale to be the one-loop 
value for A 7 v_f =4 is fixed at lOOMeV, thus giving as(M|) = 0.115.^^ This precise value is not 
determined by a best fit, but a value near to this is certainly favoured. 

Of course, since we are using a genuinely leading-order expression, any change in renormalization scale 
Mh ^ L exactly countered by a change in AiVj =4 of the form ANf =4 —> k~^ANf= 4 - However, it is 
encouraging that making the simple choice jFn = leads to a value of as{Mz) which is nicely compatible 
with the usual value. 
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There are some further details we should mention. Firstly, when obtaining a fit using the 
above approach, the values of used are not precisely those published in [1] and [2], This 

is because it is not T 2 ix, Q'^) that is measured directly at HERA, but the differential cross-section. 
This is related to the structure functions as follows 


(Pa 


27ra^ 


dxclQ'^ Q'^x 


2 — 2y + 


y 


1 + R 


Mx,Q^ 


( 6 . 1 ) 


where y = Q'^fxs ^ Q^/QOOOOx for the HERA experiment and R = P^l{x,Q'^) j{T 2 {x,Q‘^) — 
Rl{x, Q^)). Over most of the range of parameter space y is very small and Rl{x, Q^) is likely to be 
small. Hence, the measurement is essentially directly that of J- 2 {x,Q‘^). However, at the lowest x 
values, especially as increases, the value of J^ 2 {x^ Q'^) must be extracted using some prescription 
for the value of J^l{x,Q‘^)- Both the HI and ZEUS collaborations, roughly speaking, obtain their 
values of Rl{x, Q^) from predictions coming from NLO calculations of J-' 2 {x, Q'^); HI use the GRV 
structure function parameterization while ZEUS perform their own fit to the data on J- 2 {x, Q‘^)P^ 
Using these values of Pl{x, Q^) can lead to an increase in !F 2 {x, Q^) of about 12% for the highest 
values of y compared to the values of p 2 {x, Q^) obtained using the assumption that Tl{x, Q'^) = 0. 
Since the approach in this paper leads in practice to a somewhat lower prediction of Tl{x, Q^) than 
the standard NLO-in-as approach, the values of T 2 {x, Q'^) used in the fit must be altered to take 
account of the fact that our predictions for Rl{x,Q'^) are not the same as those used by HI and 
ZEUS. Thus, the R 2 {x,Q‘^) values are a little (at most about 6%) lower for the largest values of 
Q'^/x than in those in [1] and [2]. In practice, the best fit is first obtained using the published values 
of R 2 ix,Q‘^), a prediction obtained for J^l{x,Q^), the values of R 2 {x,Q^) altered accordingly and 
the best fit obtained once again. The values of R 2 {x-,Q‘^) are not altered a second time using the 
corrected prediction for Rl{x,Q^) since this changes by only a small amount, leading to further 
changes in T 2 {x, of much less than 1%. 

We should also comment on our choice of for the best fit. This is chosen to be equal to 
4GeV^ in order to obtain a reasonable description of the available data on the charm structure 
function coming from EMC [42] and also from preliminary measurements at HERA [43]. The 
quality of the fit is shown in fig. 1. As one can see the fit is of a fair quality, with the predicted 
J- 2 {x, Q"^) perhaps being a little large in general at large x, and a little small at small x. This result 
is qualitatively consistent with the fact that we have used a threshold at = m^, rather than at 

We note that both NLO calculations of T 2 {x,Q^) use a smooth treatment of the charm threshold 
while the predictions for Fl{x,Q^) use a LO formula which uses NLO parton distributions and assumes 
4 massless quarks in the expression for the longitudinal coefficient functions. Thus, there seems to be an 
internal inconsistency in this method of determining R 2 {x,Q^). Fortunately the errors due to using the 
LO formula in general cancel with those from using 4 massless quarks. The net effect is an overestimate 
of at most about 2 — 3% on the extracted value of !F 2 {x,Q^) at the largest values of y. For most of the 
parameter space the effect is negligible. 
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= Q^(x~^ — 1) = 4m^, since in the latter case at very low x the threshold will effectively occur 
at lower values of Q^, and hence .F|(x, Q^) should increase relatively, while at large x the threshold 
will be at higher values of Q^, and Q^) will decrease relatively. Hopefully, a correct treatment 

of the charm quark threshold will support these qualitative conclusions. Since our treatment of 
the charm threshold is so primitive we do not regard our value of as a determination of the 
charm mass. Rather, rric is used simply as a phenomenological parameter, ensuring that the charm 
contribution to the full structure function is qualitatively correct. Of course, = 4GeV^ is 
somewhat high compared with the values obtained from reliable determinations. We will comment 
on this value later. 

The strange quark is treated as being massless in our calculations, which is presumably a good 
approximation for the values of considered. However, we insist that the strange contribution 
to the structure function is 0.2 of the singlet structure function minus the valence contributions 
(i.e. the “sea structure function”) at = m^. This ensures compatibility with the data on 
neutrino-induced deep inelastic di-muon production obtained by the CCFR collaboration [44]. 

Finally we consider the form of the gluon at large x. Within our effective factorization scheme 
we would expect the gluon distribution to be quite similar to that in the MS scheme at NLO-in-Og 
for very large x. Thus, we demand that our gluon is qualitatively similar to that obtained from 
the WA70 prompt photon data [45] at x > 0.3. This means that the gluon distribution must be 
roughly of the form 2.5(1 — x)® at = 20GeV^ for values of x this large. Encouragingly, this is 
the type of large x behaviour that the best fit chooses for the gluon, and no strong constraint is 
needed. 

The fit is performed for a wide variety of data: the HI [1] and Zeus [2] data on 
with 0.000032 > x > 0.32 and 1.5GeV^ > > 5000GeV^; The BGDMS data [46] on J^ 2 ^{x,Q‘^) 

with 0.07 > X > 0.75 and 7.5GeV^ > > 230GeV^; the new NMG data [47] on jr^^(x,(5^) 

and T'^'^(x,(5^) with 0.008 > x > 0.5 and 2.5GeV^ > > 65GeV^; NMG data on the ratio of 

T'^"'(x,(5^) to 7^2^{x,Q‘^) [48] with 0.015 > x > 0.7 and 5.5GeV^ > > 160GeV^, GCFR data 

[49] on ^ 2 ^{x,Q‘^) and T^^{x,Q‘^) with 0.125 > x > 0.65 and 5GeV^ > > 501.2GeV^; and 

the E665 [50] data on jr^^(x,(5^) with 0.0037 > x > 0.387 and 2.05GeV^ > 64.3GeV^.^® For 
each data set the lowest x bins cover a range of from the minimum up to somewhat less than 
the maximum, and similarly the highest x bins cover a range of from somewhat higher than 

It has recently been realized that the MRS fit program has been treating the errors for the E665 data 
incorrectly; the errors as a fraction of the value of A 2 (x,Q^) have been interpreted as the absolute errors, 
and the errors used have therefore been on average about 3 times too big. This has resulted in the very 
low in, for example, [26] and [31]. In this paper the correct errors are used, though the MRSR fits are 
not performed again. As will be seen, including the correct errors raises the for the fit to the E665 data 
by a factor of about 8, as we would expect. However, it does not really affect any comparison of the best 
fits since all the fits we perform give descriptions of this data which are almost identical in quality. 
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the minimum up to the maximum. We use the NMC data on the ratio J- 2 ^{x, Q‘^)lJ- 2 ^{x, Q^) as 
well as on T 2 ^{x, Q"^) and Q'^) since the parameter space covered by the two types of data 

sets is by no means identical. One can see that the full range of data used in the fit covers an 
extremely wide range of both x and and thus provides a very stringent test of any approach 
used to describe it.^^ 

The result of the best fit to this data using the leading-order, including leading-ln(l/x) terms, 
RSC expressions (henceforth referred to as the LO(x) fit) with Qq chosen to be equal to 40GeV^, 
= 4GeV^ and mf = 20GeV^ is shown in table 1. Also the result of the fit to the small-x data 
is shown in fig. 2. As one can see there is a very good quality fit to the whole selection of data, 
and thus over the whole range of x and Q^. Only the NMG Q^) data give a 

of much more than one per point, but this is due to the scatter of the data points and is true for 
any global fit. Overall, the fit gives a of 1105 for 1099 data points.^® 

The fit shown is for the particular starting scale Qq = 40GeV^, but the quality of the fit is 
extremely insensitive to changes in this scale, as we expect from the method of construction of the 
solutions. The fit is essentially unchanged over the range 20—80GeV^, and we choose 40GeV^ as the 
geometric mean. When Qq drops below 20GeV^ the fit immediately gets markedly worse. This effect 
is due to the bottom quark threshold, and we can explain it briefly. Both dJ- 2 {x,Q'^)/dlnQ‘^ and 
Tif {x,Qq) are nearly proportional to T£(x, Q^q)) hence nearly proportional to Nf at small x, as 

can be seen from (4.180) and (4.178) and from (4.36) and (4.37). Thus, when starting at Qg above 
= 20GeV^ they are both oc. Nf = 5. As we go below the bottom quark threshold, = 20GeV^, 
dT 2 {x,Q'^)/dlnQ‘^ becomes smaller, and is oc 4/5.F|'(x, Qg). However, if Qq < 20GeV^, then 
dT 2 {x,Q‘^)/dlnQ‘^ is proportional to {x,Qq), with the same constant of proportionality as 
above, when below = 20GeV^, and is oc 5/4.T^(x, Qg) above = 20GeV^. In the latter 
case, choosing the correct value of fo fff fh® data leads to dN^{x,Q‘^)/dlnQ'^ being too 

large at very small x. The quality of the fit then gets continuously worse as Qq lowers further. 

We have not included ZEUS data for > 2000GeV^ since this was not included in the MRSR fits 
[26] against which we compare our results. The fit to this high x, high data is practically identical for 
all the fits discussed in this paper, giving a x^ of about 17 for the 7 data points. This fairly large value of 
X^ is due to the large scatter of these points. 

The LO(x) fit is a little different from that in [31] since rric = 4GeV^ is chosen here in order to give 
a rather better description of the charm structure function data. In fact this value of rric leads to the best 
overall fit for the data also: lower leads to a slightly worse fit to the HERA data with the fit to the 
other data remaining more or less unchanged, while higher leads to a slightly worse ht to the large x 
data. We also point out that there was a problem in the computer program when at < mi in [31]; the 
slope of dd^ 2 {x, Q^)/dlnQ^ was a little too small in this region. Thus, in fig. 1 in [31] this slope should be 
a little larger below = 3GeV^, and hence the kink at this value of becomes less obvious. As far as 
the ht is concerned, the effect of this error is minute since the error is small, and because there are so few 
data points below = 3GeV^. In practice the overall should be ~ 2 better than quoted in [31]. 
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becoming completely uncompetitive long before reaching We expect that a correct treatment 
of quark thresholds would lead to similar results to above for Qq 40GeV^ {(HF 2 {x,Q‘^)/dlnQ^ oc 
at the input here, but {x,Q'^)/dlnQ‘^ falls off towards 3/4 of this as falls to 
~ 4GeV^ and the charm quark effects die away) and then a smooth falling off of the quality of the 
fit as Qq was lowered, with it again beginning to deteriorate somewhere in the region of 20GeV^, 
due to asiQo) becoming too large below this value. 

There are 18 free parameters used in the fit; ^qcd, 4 : (which we choose to describe as a 
parameter since, although it is fixed at lOOMeV, and small variations would no doubt improve 
the fit slightly, it can certainly not vary much); four parameters for each of the valence quark 
contributions, which are described by functions of the form (4.38) where the normalization is set by 
the number of valence quarks; four parameters for the nonperturbative inputs for the two singlet 
structure functions, which are of the form (5.9); and the unknown scale All, where we allow All 
to be a free parameter for each Qg. We do not consider Qo a free parameter, since it can take a 
very wide range of values. The parameter All-, which we have argued should be a scale typical of 
soft physics, turns out to be 0.55GeV^ for the fit starting at Qo = 40GeV^. This decreases a little 
as Qg increases and vice versa. For Qg = 40GeV^ the soft inputs for the fit are roughly 

Ff (x) w 3(1 - x)^(l + 0.1x°-^ - 0.2x), F^{x) (1 - - 0.65x°-® + 4.5x), (6.2) 

where they have been forced to be flat as x —> 0. These nonperturbative inputs lead to complete 
inputs of F£(x, Qq) k. 3x“‘^'^^ and F^ix, Qq) ^ 0.65x“°'^ for 0.01 > x > 0.0001, with the effective 
A increasing in both cases for even smaller x. Instead of forcing the nonperturbative inputs to be 
flat as X —> 0 we could allow an asymptotic behaviour x~^, where A ^ 0.08. This leads to an 
equally good fit as for the flat nonperturbative inputs. 

Thus, the LO(x) fit seems to be a success. However, in order to gauge its true quality it 
is helpful to have some points of comparison. Hence, we will discuss some alternative fits. As 
in [31], the most recent MRS fits Ri and R 2 are shown. These are obtained using the standard 
two-loop method, where Ri allows to be free (giving A^~"^ = 241MeV) and R 2 fixes 

= 344MeV to force a better fit to the HERA data. The new NMG data for F|‘^(x, Q"^) and 
F^^(x, Q"^) are used rather than the previous sets which were used in the MRS fits. This new data 
is not used directly in the MRS fits, i.e they are not performed again using these new data, the 
values for the new NMG data sets are simply calculated using the same input parameters as in [26]. 
However, there is little indication that these would be changed much by the new data. The MRS 
fits are useful for purposes of comparison for a number of reasons. Firstly the treatment of the 
errors in this paper is identical to that in the MRS fits, so consistency in this respect is guaranteed. 
Also, any systematic differences due to differences in computer programs is guaranteed to be absent. 
Finally, the cuts in Q^ for each data set are chosen to be the same in this paper as for the MRS 
fits. 
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The number of free parameters in the NLO-in-a^ fit is the same as in the LO(x) fit. There is 
one more parameter in the NLO-in-a^ fit due to the powerlike forms of the input gluon and singlet 
quark at small x being independent^® whereas the small-x shape of the inputs for both structure 
functions is completely determined by the one parameter All in the LO(x) fit. However, there 
is one less parameter in the NLO fit due to the normalization of the gluon being determined by 
momentum conservation. Essentially, the LO(x) ht is a little less constrained at large x than the 
NLO fit, but rather more predetermined at small x than the NLO fit. Indeed, the LO(x) fit has 
a certain amount of predictive power in the low x region, as seen at the end of the last section, 
whereas the NLO has none. 

Comparing the LO(x) fit and the MRS fits, it is clear that the LO(x) scheme-independent fit 
is much better for the HERA data (even when compared to R 2 ), much better for the BCDMS data 
(even when compared to Ri) and similar in standard for the rest of the data. The overall ht is 
~ 200 better for the whole data set; clearly a lot better. However, this direct comparison with the 
MRS hts is rather unfair, since there are a number of ways in which it differs from the LO(x) ht. 
Eirst, and not terribly importantly, the normalizations of the data sets are allowed to vary more 
in the LO(x) hts than in the MRS hts, as described in the caption to table 1. (There is a small 
systematic difference in the normalization required by the HI data and ZEUS data.) Rather more 
importantly the MRS hts use the SLAC data [51] on J^ 2 ^(x, Q^) in the hts as well as the data sets 
mentioned in table 1. This covers the range 0.07 > x > 0.65 and 2GeV^ > > 22.2GeV^, i.e. the 

same sort of high x range as the BCDMS data, but rather lower in general. Hence the SLAC 
data is not included in the LO(x) ht due to much greater sensitivity of this data to potential higher 
twist effects than any of the other data sets. If it is included in a ht, then in certain regions of 
parameter space it is incompatible with the BCDMS data, as can be seen in hg. 4 of [26]. Hence a 
best ht to both the BCDMS data and the SLAC data must, to a certain extent, be a compromise 
between the two. Therefore the ht to the BCDMS data in the MRS hts is worsened somewhat by 
the need to ht the SLAC data. If the SLAC data were ignored the hts to the BCDMS data would 
improve by about 40. The LO(x) ht also differs from the MRS hts in the starting scale of the 
evolution and in the treatment of the quark thresholds. The charm contribution to the structure 
function in the MRS hts turns on at Qq = IGeV^ but is suppressed by a phenomenological factor 
of a function of (3.5GeV^/(5^) in order to give a reasonable description of the charm data. The 

We claim that the allowed powerlike form of the inputs at small x is against the spirit of a well- 
ordered perturbative expansion, and that the small-x form of the inputs should be determined in terms of 
the nonperturbative inputs and the requirement that the expressions are insensitive to the choice of Qq, as 
already discussed in §4.2. However, enforcing this rule would mean that the quality of the fit to the small-x 
data was very poor. Thus, we allow the more usual, unjustihed choice for the small-x form of the inputs. 
In fact, at Qq = IGeV^ the gluon is strongly valence-like, while the sea-quark distribution ~ a:-i-° i"at 
small X. 
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bottom contribution to the structure function, which is very small, is not included at all. Finally, 
the fact that the quark and gluon distributions are input at IGeV^ in the MRS fits is not helpful 
to the quality of the fit. If the best fit is obtained for the HERA data using the NLO calculation 
with massless quarks with inputs at Qq ~ 4GeV^, as in [27], and evolution performed downwards 
in then the gluon becomes negative at very small x by the time = IGeV^ is reached. Hence, 
the valence-like input gluon distributions in MRSRi and MRSR 2 are not the ideal distributions 
for fitting the HERA data. A parameterization which would allow the gluon distribution to turn 
over and become negative at very small x would be better, as would simply starting at larger 
where the gluon distribution is happy to be positive everywhere.^*^ 

Thus, in order to make a more meaningful comparison of the LO(x) fit with a NLO-in-ojs fit 
we have ourselves performed a NLO-in-as fit, called NLOi, allowing the normalizations to vary 
in the same way as in the LO(x) fit, with exactly the same treatment of quark thresholds as the 
LO(x) fit, with the input parton distributions chosen at Qq = and with the fit to exactly the 
same data, i.e. to those sets shown in table 1, and with the constraint on the gluon from the WA70 
prompt photon data.^^ As in the MRS fits the values of J^ 2 {x,Q‘^) at small x are those quoted 
in [2] and [1]. Strictly speaking they should be altered to take account of fact that the predicted 
values of Tl{x, Q^) from these hts are not precisely the same as those used in [1] and [2] to extract 
Q^) from the measurement of the differential cross-section, i.e. an iterative procedure should 
be used as with the LO(x) fit. However, the differences between the values of J^l{x, Q^) in [1] and 
[2] and from these NLO-in-ojs fits are far smaller than the differences between those in [1] and [2] 
and the LO(x) fit. Hence, in practice the errors introduced by not altering the values of Tl{x, Q^) 
are very small (the for the fit to the HERA data might improve by a couple of points if the 
formally correct procedure was used). The value of is chosen to provide a good description of 
the data on the charm structure function. The value needed is = 2.75GeV^, somewhat lower 
than in the LO(x) fit, and the fit to the charm structure function data is shown in fig. 3. The 
value of A^E~^ determined by the fit is intermediate to those chosen by the MRS fits, being equal 
to 299MeV. The quality of the fit is shown in table 1. 

Starting the evolution at Qq < IGeV^ quickly leads to much worse fits to the data. Also starting at 
fairly high Qg, i.e. Qq ^ 20GeV^, leads to a worsening of the fit due to the given form of parameterization 
not being adequate to describe the parton distributions at these scales, see e.g. [28]. Thus, Qq is no less a 
parameter in a NLO fit as in the LO(x) fit. 

Once again we allow the small-x form of the inputs to be unjustified powerlike behaviours. At 
Qq = 2.75GeV^ the gluon is quite flat while the sea-quark distribution ~ 

As in the LO(x) fit the choice of rric which gives a good description of the charm data is also the 

o f ^—4 

choice which gives the best global fit. Raising the value of leads to the value of A^^ chosen going up, 
which improves the fit to the HERA data, but the fit to the rest of the data deteriorates by more than this 

r\ r - ^ 

improvement. Lowering the value of leads to the value of Al^ chosen going down, improving the fit 
to the large x data, but again leading to overall deterioration. 
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As one can see, the different treatment of the quark thresholds and the higher starting scale 
for evolution has produced a better fit to the HERA data than the MRS fits, even though the value 
of as((5^) is lower than in the R 2 fit. However, the best fit comes from allowing the normalization 
of the HI data to be at the lower limit allowed by the error on the normalization^^, and is still not 
as good as the LO(x) fit to the HERA data. Not including the SLAC data, and in the case of the 
R 2 fit the lowering of leads to a much better fit to the BCDMS data than the MRS fits, 

but again this is clearly not as good as the LO(x) fit. The fit to the NMC data is much the same 
for the LO(x) fit as for the NLOi fit, and the NLOi fit is a little better for the CCER data. The 
overall quality of the NLOi ht is 1169 for the 1099 points, and thus is 64 worse than the LO(x) fit. 
Therefore, the LO(x) fit is still clearly better than the NLOi ht; but not nearly as convincingly as 
appeared to be the case when compared to the MRS fits. Nevertheless, it is encouraging that, while 
the overall fits to the relatively high x data, i.e. the BCDMS, NMC, E665 and CCER data are 
similar in the NLOi ht and the LO(x) fit, it is the fit to the small-x HERA data that is definitely 
better for the LO(x) fit, as we would expect. This can be seen even more clearly if we simply 
examine the quality of the fit by separating the data into two sets: one where x < 0.1 and one 
where x > 0.1. This is shown in table 2, and demonstrates that the LO(x) fit is superior at small 
X, while not quite as good as the NLOi ht at large x. 

This qualitative result above is exactly what we would expect. The importance of the leading- 
ln(l/x) terms in the LORSC calculation can be quantitatively judged by how they affect the fit. 
If, after obtaining the best LO(x) fit, all terms other than those in the one-loop expressions are set 
to zero, the quality of the fit is unchanged above x = 0.3, begins to alter slightly below this, and 
is clearly worse by the time we reach x = 0.1. Thus, the leading-ln(l/x) terms are important by 
this value of x. However, much of this effect is due to the terms at so the NLO expression 

at fixed order in ag should be insensitive to higher-order-in-Og leading--ln(l/x) terms down to 
X somewhat lower than 0.1. Therefore, above x = 0.1 the NLO fit should in principle be better 
than the LO(x) fit since it contains terms at NLO in ag which are important at large x. However, 
the NLO fit should be considerably worse at small x since it does not contain many important 
leading-ln(l/x) terms. This is qualitatively in agreement with the comparison of the NLOi ht and 
the LO(x) fit. 

However, the above comparison is somewhat incorrect because in the process of obtaining 
the best fit for all the data the NLOi ht may choose some parameters, particularly A^L~^, such 
that they mimic the effects of the leading“ln(l/x) terms, and a decent fit for the small x data is 
obtained to the detriment of the fit to the large x data. In order to check this hypothesis we have 
also performed a NLO-in~as fit with A^L^^ fixed at 250MeV, which we will denote by NLO 2 . The 

The fit to the HI data continues to improve slightly for a normalization going down to 0.96 if this is 
allowed. 
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results of this fit in terms of the different data sets is shown in table 1. The fit clearly improves 
compared to the NLOi fit for the BCDMS and CCFR data, and gives the best overall fit for the 
high X data sets. It worsens somewhat for the NMC data, and also for the HERA data, and overall 
is slightly worse than the NLOi ht, having a of 1184 for the 1099 data points. Nevertheless, it 
is perhaps a truer representation of a real NLO-in-cis ht than NLOi since, as can be seen in table 
2, it gives a better ht to high x data, but not such a good ht to the small-x data which presumably 
require the leading”ln(l/x) terms. Whether one believes this argument or not, it is certainly clear 
that the LO(x) ht does provide a better description of the data than any standard NLO-in-a<j ht. 

We can be even more general in our NLO-in-o;* type ht. We have already argued that, when 
working at hxed order in ag, despite the fact that the amount of momentum carried by the partons 
can be constant, it does not have to sum to unity, i.e. higher-order inputs can carry some fraction 
of the momentum. Thus, we have also performed a ht where the total momentum of the partons 
is not constrained. When doing this we have again allowed the small~x form of the parton inputs 
to be oc which, as we have already stressed, we do not believe is correct procedure for a 

NLO-in-ag calculation. However, we will allow as much freedom in this ht as possible. Doing so 
we obtain a ht which is nearly, but still not as good as the LO(x) ht. The ht to the data for x > 0.1 
is about 25 better than the LO(x) ht, i.e. similar to the NLO 2 ht, but still about 45 worse for 
X < 0.1. The value of chosen by the ht is 230MeV. The amounts of momentum carried 

by the quarks are similar to the other NLO-in-cTs hts; these are constrained by htting the data 
near to the input scale. The momentum carried by the gluon is 57%, far higher than in the other 
NLO-in-CKs hts. This large amount of gluon momentum is necessary in order make the gluon large 
at small x and hence provide a good ht to the small-x data for even with the small 

value of The large amount of gluon momentum leads to the total momentum carried by 

the NLO partons being 114% of the proton’s momentum. We do not believe the validity of this 
ht, having already argued that the LORSC expressions are the theoretically correct expressions, 
and that leading~-ln(l/x) terms are necessary for the correct description of small-x data. Indeed, 
this ht leads to very different predictions for other quantities from those obtained using the LO(x) 
ht, as we will soon see. However, we do believe that the lack of the momentum constraint on 
the inputs allows this NLO-in-CKg ht to be perhaps the most realistic at high x: the fact that the 
NLO-in-Q;<j ht is not designed to work at small x due to the lack of leading-ln(l/x) terms can be 
largely compensated by an overly large input gluon distribution at small x, and the high x part of 
the ht can be as good as possible without having to compromise itself to match the low x data. In 
fact, at X > 0.1 the details of this ht are very similar to the NLO 2 ht, but the gluon is much larger 
at smaller x in this ht. 

Let us comment briehy on other possible NLO-in-a* hts. As already mentioned it is difhcult 
to compare directly with global hts from groups other than MRS since there may be differences 
in treatment of errors or systematic differences in the outputs of the computer program. However, 
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we can make some comments. The latest CTEQ fits [52], e.g. CTEQ4M, are performed in a very 
similar manner to the NLOi fit, though with different cuts on data, and the results also appear to 
be very similar.Hence, there seems little more to say than that the NLOi fit is similar to the 
CTEQ results. 

We also feel we must comment on the ht performed by HI in [Ij. This takes account of the 
charm quark threshold rather more correctly than in the present paper, i.e. there is a smooth 
threshold at = 4m^. However, it is a fit to only the HI, NMC and BCDMS data. It is difficult 
to make a direct comparison with this fit due to the completely different way of performing the 
calculation. It is however possible to perform a NLO-in-Os fit to the same data as the H! fit using 
the same prescription as for NLOi and NLO 2 . This has been done, and many similarities to features 
of the HI fit noticed. It is possible to obtain a NLO-in-cts fit for the HI, NMC and BCDMS data 
of a comparable quality to that in the LO(x) fit. However, if one simply inserts the values of the 
ZEUS and CCER data after finding the best fit one finds that the for the ZEUS fit is over 300, 
i.e. there is a certain degree of incompatibility between finding the best fit to the HI data and the 
ZEUS data, and the fit to the CCER data is very poor indeed. Perhaps more importantly, since 
the fit does not put any direct constraint on the gluon, then as in [1] one obtains a gluon which at 
= 5GeV^ behaves like 2x~^'^ at small x, which helps the fit to the HI data, but behaves like 
2(1 — x)® at large x. This leads to a gluon like 2(1 — at large x at = 20GeV^ rather than 
the 2.5(1 — x)® required by the WA70 prompt photon data. Eor values of x 0.4, 2(1 — x)®'® is 
clearly a great deal smaller than 2.5(1 — x)®, and so the gluon obtained from this type of approach 
is hopelessly incompatible with the prompt photon data. 

We feel that similar, if less dramatic, considerations are probably true of the type of approach 
adopted in [20], [27] and [28], where fits are performed to small-x data alone. Despite the fact that 
the input parton distributions are constrained to be very similar to standard parton distributions 
at large x, e.g. the MRSDO distributions, it must be remembered that these parton distributions 
provide good descriptions of the large x data when evolved at NLO in a* only when using a 
particular value of Evolving the MRSDO parton distributions, which should have A^L~^ = 

230MeV, using instead = 360MeV will produce very different results. Indeed, as some sort 

of comparison we are able to obtain a fit to the HI data and NMC data (thus constraining the 
high X parton distributions) alone to a quality rather better than in the full LO(x) fit by using 
a NLO-in-Q!s fit. However, the fit does indeed choose = 360MeV. The resulting fit to 

virtually all the other data is very poor. 

There is a clear improvement to the fit to the BCDMS data compared to the NLOi fit, and clear 
worsening of the fit to the CCFR J^ 2 (x,Q^) data. These seem to be typical systematic differences between 
all CTEQ fits and MRS type fits and are presumably due to systematic differences in the treatment of 
errors or in the numerical calculation. The comparison of the codes in [34] certainly supports the latter 
conclusion. 
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Thus, we have a number of clear demonstrations of the dangers of performing fits to only a 
restricted set of data. In order to avoid very obvious inconsistencies with some particular data it is 
necessary to perform proper global fits. One may investigate the possibility that the data in some 
region of parameter space is not really expected to be fit well using a particular expression. Indeed, 
this is what we have done in the NLO 2 fit, and believe that we have sound theoretical arguments 
for doing so, and moreover that the results of the fits back these up. 

This leads us to the question of the determination of as(M|) using global fits to structure 

function data. The complete RSC expression for structure functions only exists at leading order, 

and, as already mentioned, this leaves the renormalization scale completely undetermined. Hence, a 

determination of as{M^) does not really take place. (Note, however, the previous comments on this 

question.) It is not yet possible to extend the RSC calculation beyond the leading order due to lack 

of knowledge of NLO-in-ln(l/x) terms, but hopefully these will shortly become available [35], and 

when they do the NLO versions of (4.177)-(4.180) can be derived and put to use. Only then should 

the NLO coupling constant really be used in any global fit. Until the full renormalization-scheme- 

consistent NLO expressions become available, we believe that it is incorrect to use NLO-in-as fits 

to small-x structure function data in order to determine the NLO coupling constant (unless, of 

course, direct measurements of and other less inclusive quantities at very small x turn 

out to verify standard two-loop predictions). However, as already argued, the fixed-order-in-ag 

expressions should be accurate for CCFR, BCDMS and NMC data (except perhaps at the lowest 

X values), which after all are still much more precise than HERA data, and fits to these data alone 

will provide the best determination of the NLO as{Q‘^). However, let us be wary here. One might 

think that the above comments mean that we believe that the value of at NLO is about 

MS 

250MeV, since our fits to large x data at NLO-in-o;* support this. To a certain extent this is true, 
but we also believe that our naive treatment of the heavy quark thresholds can introduce an error 
in the determination of as{M^) which could easily be as large as 0.005. Hence, until this naive 
treatment is improved, we have nothing concrete to say about the value of as{Q‘^). 

We have given what we hope are convincing arguments for our advocated approach for cal¬ 
culating structure functions. Not only do we claim theoretical correctness, and limited predictive 
power, but we also have good quality, very comprehensive fits to data on However, we 

are well aware that only further experimental tests can prove us right (or wrong). Hence, we now 
discuss our predictions for J^l{x,Q'^). 

So far we have only probed J^l{x, Q'^) indirectly, i.e. it is simply related to the hiQ^ derivative 
of J^ 2 {x-,Q‘^) (as well as to the input !F 2 {x,Qq)). Having tied down the nonperturbative inputs 
and All and Qq from our fit to T 2 {x, Q'^), we have a prediction for J^l{x, Q^)- The result of this 
prediction for the fit with Qq = 40GeV^ is shown in hg. 4, where it is compared to the prediction 
using the NLO-in-ag approach, in particular the NLOi fit. As one can see, it is smaller than the 
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NLOl but becomes steeper at very small x. The NLO 2 fit gives a very similar form 

of J-l{x,Q'^) to the NLOl ht, while the NLO-in-as fit where momentum is not conserved leads 
to a similar form of at high x but an even larger value than the NLOi fit at x < 0.1. 

The LORSC prediction for Q^) is weakly dependent on the value of Qg chosen: the value at 

= 5GeV^ and x = 10“^ varies by ±10% within our range of Qq (increasing with Qq), and by less 
than this for higher x and The very recent results on for 0.01 ^ x ^0.1 from NMC 

[47] are matched far better by the LO(x) Tl{x,Q'^) than the NLOi Tl{x,Q‘^) (the latter being 
rather large). However, it is fair to say that any problems with the NLOi ^Lix,Q‘^) can almost 
certainly be attributed to the naive treatment of the charm quark threshold, i.e. the predicted 
J-l{x,Q‘^) from the NLO-in-a., calculation in the last of [9] matches the data well. Actually, it is 
rather uncertain whether in this region of parameter space ±l(x, would be better described by 
the LORSC calculation, which ignores some 0(a‘^) effects, or the NLO-in-ccg calculation, which 
ignores the leading-ln(l/x) effects beyond second order in 

Measurements of Tl{x, at x < 10“^ would be a better discriminant between fixed-order- 
in-a^ calculations and those involving leading-ln(l/x) terms. However, the sort of “determination” 
of J-l{x, Q'^) already performed by HI [54] is really only a consistency check for a particular NLO- 
in~as fit. It is by no means a true measurement of J-l{x, Q^). In essence, all it proves is that the 
measurements of the cross-section are consistent with a particular NLO-in-ccs fit to J^ 2 {x, Q'^) when 
the relationship between the cross-section and the values of J- 2 {x, Q'^) is determined assuming the 
correctness of the NLO-in-Og expressions for both structure functions. This is a perfectly correct 
procedure, and should be adopted for any fit to J^ 2 {x,Q‘^) data, as it has been for the LO(x) fit. 
It would cause concern about the validity of the particular fit being used if it were to fail, but in 
itself says nothing whatever about the validity of a different approach, or about the actual value 
of J^l{x, Q^)- Hence, real measurements of J^l{x, Q^) are needed in order to find the real values of 
J-l{x,Q‘^). From fig. 4 it is clear that such measurements at HERA would be an important (and 
probably essential) way of determining the validity of the approach in this paper, and the genuine 
importance of leading-ln(l/x) terms in structure functions. 

Another potentially important discriminant between different methods of calculating structure 
functions is the measurement of the charm structure function. Both the LO(x) fit and the NLOi 
fit can provide good fits to the currently available charm data, as already seen, and the value of 

The author has submitted two conference proceedings on the topic of the current paper [53] and should 
point out that these are both incomplete. In the former the expression (4.180) was not used for the input 
for ± 2 (x, Q^). Thus, the fit imposed far less constraint on Qq than the full procedure, and a value of 5GeV^ 
was used, which resulted in a prediction of that is much too small. The latter claimed that an 

input for ±l(x, Q^) a little smaller than that consistent with the full set of expressions was needed for the 
best fit, even at the optimum Qq. This was due to there being no account whatever taken of the 6-quark 
threshold. The correction makes very little difference to the phenomenological results in this latter case. 
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in the NLOi fit is rather more satisfactory than that in the LO(x) fit. However, these calculations 
involve incorrect treatments of the charm quark threshold. A correct treatment of this threshold 
at NLO in [41] shows that the value of must be at least halved in order to produce the same 
sort of value of the charm structure function as the approach used in this paper. Indeed, as seen 
in [55], correct NLO-in-cis calculations with = 2.25GeV^ undershoot the small-x data (other 
than the HI, calculation which has a form of the gluon already criticized and which would certainly 
badly undershoot the larger x EMC data). In fact, standard NLO-in-as calculations using the 
coefficient functions for charm production in [41] only produce a sufficiently large charm structure 
function if ~ IGeV^ [56], i.e. somewhat lower than expected. 

A correct treatment of the charm quark threshold within the framework advocated in this 
paper has not been completely worked out (although, as already mentioned, work is in progress), 
and has certainly not been tested. However, using the naive treatment in this paper, the quark 
mass required in the LORSC approach is somewhat higher than for the NLO-in-ag approach. If 
the same sort of conclusion is true when using a correct treatment of quark thresholds, as we might 
expect (but cannot guarantee), then this will offer further support for the LORSC calculation. More 
generally, a comparison of the theoretical calculations with the ever-improving data on the charm 
structure function seems potentially to be a very useful way of discriminating between different 
methods of calculation.^® 

In principle there are many other quantities which could be calculated within the LORSC 
framework and compared to those calculated using the NLO-in-a^ approach (or any other method) 
and to experimental data. Particularly obvious examples are the distribution of the transverse 
energy flow in the final state in lepton-hadron scattering and the the cross-section for forward 
jet production, for both of which there exists some experimental data which does not seem to be 
terribly well described by the order-by-order-in-ag approach, particularly in the latter case, see 
e.g. [57] and [58]. We have not performed a LORSC calculation of any such quantities, and it is very 
difficult to estimate the results, other than guess that they will probably lie somewhere between the 
fixed-order-in-ag predictions and those obtained using BFKL physics naively. Such calculations are 
obviously a priority, and work will begin soon. Only by comparing our theoretical predictions with 
a wide variety of experimental data can we determine unambiguously which theoretical approach 
is correct. For the moment we leave the theoretical arguments and the results of tables I and 2 as 
the evidence supporting our particular approach. 

We also note that as seen in [41] the longitudinal structure function is more sensitive to the treatment 
of the charm quark threshold than !F 2 {x,Q^), with there being a sizeable discontinuity in !Fl{x,Q^) at 
= ml when using the approach in this paper. Thus, for accurate predictions of !Fl{x,Q^) at values 
of in the vicinity of a smooth treatment of the charm quark threshold is really needed in practice. 
The problem is not as acute at the bottom quark threshold because the total squared charge of the quarks 
changes proportionally much less at this threshold than at the charm quark threshold. 
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7. Conclusion and Summary. 

In this paper we have derived expressions for the structure functions Q^) and 

in a theoretically correct, well-ordered manner. We have first done this for the particular expansions 
schemes which order the expressions strictly in orders of a*, i.e. the standard loop expansion, or in 
terms of the leading powers of ln(l/x) for given power of a^, i.e. the leading-ln(l/x) expansion. In 
both cases we have demonstrated that a correct calculation automatically leads to factorization- 
scheme-invariant results which may be expressed in terms of Catani’s physical anomalous dimen¬ 
sions [30]. Thus, we refute any suggestion that these physical anomalous dimensions are simply 
an example of anomalous dimensions in a convenient, physically motivated factorization scheme, 
but insist that they are fundamental pieces in the correct expressions for the structure functions. 
Likewise, we insist that no particular factorization scheme is fundamentally any more “physical” 
than any other. However, we also demonstrate that in the case of the leading-ln(l/x) expansion the 
correct expressions are more difficult to obtain than in the loop expansion, and a correct calculation 
requires more than just the knowledge of these physical anomalous dimensions. 

We have then argued that both the above expansion schemes are restrictive, and lead to only 
part of the correct solution at any given order. We have shown that the only calculational method 
which is truly consistent with working to a given order in ag within a given renormalization scheme 
is the renormalization-scheme-consistent expansion, in which one works to a given order in both 
as and in ln(l/x) for given power of ag- After presenting our argument that this is indeed the only 
correct expansion scheme, we have made use of Catani’s physical anomalous dimensions (though 
this is not necessary, merely very convenient) to derive the leading-order, renormalization-scheme- 
consistent expressions for the structure functions T 2 {x,Q'^) and J^l{x,Q‘^). 

As part of our overall approach we have also taken a different view of the starting scale 
from that normally taken. Rather than trying to guess some form for the input at some particular 
Qq, or simply allowing the inputs at some arbitrary Qo take any form they wish within a 
given parameterization, we have demanded firstly that our expressions should be as insensitive as 
possible to the choice of the input scale, and secondly that any deviation from the flat Regge- 
type behaviour of the structure functions must come from perturbative effects. Thus, our inputs 
take the form of nonperturbative functions, which are flat at small x, convoluted with functions of 
the physical anomalous dimensions which are evaluated at Qo determined by the requirement 
of insensitivity to the value of Qq. This leads to our inputs being determined entirely in terms 
of the flat nonperturbative inputs and one arbitrary scale All which roughly indicates the scale 
where perturbative physics should break down, i.e. All ^ IGeV^. This gives us a great deal more 
predictive power than more usual approaches. We have some idea of the form of individual structure 
functions at small x: for a full range of sensible choices of Qg (lOGeV^ — lOOGeV^) and All we 
obtain JT 2 (x, Qq) roughly oc for 0.01 > x > 0.00001, which is clearly in good qualitative 

agreement with the data. However, we also have a much stronger prediction for the relationship 
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between the small-x inputs for J^ 2 {x^Q‘^)^ !Fl{x,Q^) and for dJ^ 2 ix,Q‘^)/dlnQ‘^. Therefore, as 
well as some predictive (or at the very least, explanatory) power for individual structure functions, 
we have a consistency condition between the form of the different inputs, i.e. once we choose the 
input for J- 2 {x,Q‘^) we have determined, up to a small amount of freedom, the small-x inputs for 
dT 2 {x,Q‘^)/dl'nQ‘^ and ^l(x,( 5^). This is a feature unique to the approach in this paper. 

Not only are the features of the LORSC calculation unique and compelling, but they also work 
rather well in practice. The LORSC expressions, including this constraint on the small-x inputs, 
lead to very good fits to the data. To qualify this statement, the for the LORSC fit to 1099 data 
on !F 2 {x,Q‘^) ranging from from 0.75 >x> 0.000032 and 1.5GeV^ < < 5000GeV^ is better by 

at least 60 than any NLOGn~as fit, even though the NLO-in-ag fit is allowed arbitrary, unjustified 
powerlike behaviour at small x (i.e. the sea-quark distribution cx x~^~^'‘^), and the small-x inputs 
for J^ 2 {x,Q‘^) and d!F 2 {x,Q‘^)/dlnQ‘^ are allowed to vary with respect to each other a great deal 
more than in the LORSG fit. In fact, as we would hope, all of this superiority comes from the fit to 
the data with x < 0.1. Moreover, much the same quality of fit is obtained for a range of Qq from 
20GeV^ — SOGeV^. Let us also put the difference in the quality of the LO(x) and the NLOi fits, i.e. 
a difference of of 64, into context; if one obtains fits to the data using the NLO-in-as approach 
and allows the fit to become up to 64 worse than the (best) NLOi fit, then values of from 
190MeV to 380MeV are allowed, i.e. as(M^) = 0.109 ^ 0.122. As a caveat to the above, however, 
it is certainly true that the calculations in this paper must be improved to take account of massive 
quark thresholds in a better manner (as is the case with most NLO-in-a^ fits), and work towards 
this end is in progress. Nevertheless, with the present treatment we feel that the quality of the fit 
and the degree of explanatory (if not predictive) power, not to mention the theoretical correctness, 
give strong justification for the LORSG expressions. 

We do, however, recognize that the quality of the fit alone, although impressive, is not such a 
substantial improvement on more standard approaches that it necessarily convinces one that this 
approach has to be correct. In order to obtain verification we must compare with more and different 
experimental data. Hence, we have presented a LORSC prediction for !Fl{x,Q^), comparing it 
to that obtained using the NLO-in-ag approach. Hopefully there will be true measurements of 
J-l{x, Q^) at HERA some time in the future with which we can compare these predictions. We stress 
the importance of such measurements to the understanding of the physics which really underlies 
hadron interactions. In the near future we will also have predictions of the charm contribution 
to the structure function within the framework of a correct treatment of the massive quark, and 
comparison with the ever-improving data on the charm structure function should also be a good 
discriminant between different theoretical approaches. Calculations and measurements of other, 
less inclusive quantities, such as forward jets, are another clear goal. 

Finally, as we have noted, our calculation is at present only at leading order due to the lack of 
knowledge of next-to-leading~order-in”ln(l/x) coefficient functions and anomalous dimensions, or 
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equivalently of physical anomalous dimensions. This means that the NLO-in-Og approach is still 
in principle superior to our approach at high x, where the leading~ln(l/x) terms at third order in 
as and beyond are not important but the 0{ag) terms are. Indeed, in practice the NLO-in-a^ 
fits, particularly the NLO 2 fit, are slightly better than the LO(x) fit for data at x > 0.1. We 
might therefore expect any predictions coming from the NLO-in-a.; approach to be more accurate 
than those coming from the LORSC approach down to values of x somewhere in the region of 
0.05. The lack of the NLORSC expressions also means that a true determination of as(M^) 
from a global fit to structure function data is not yet possible, but that the best determination 
from fits to structure function data should at present come from using the NLO-in-Os approach, 
with a correct treatment of heavy quark thresholds, using only large x data. For a really fair 
comparison between the renormalization-scheme-consistent method and the conventional order- 
by-order-in-Os approach we really need the full next-to-leading-order, renormalization-scheme- 
consistent calculation. Hopefully the required NLO in ln(l/x) quantities will soon become available 
[35], and with a little work we will be able to make such a comparison. We would expect that when 
using the NLORSC expressions the fit to the large x data would become at least as good as for the 
NLO-in-Q!s approach, and that the fit to the small-x data would be of at least the same quality 
as for the LORSC fit. We wait expectantly to discover if this is indeed the case. 
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Table 1 

Comparison of quality of fits using the full leading-order (including leading-ln(l/x) terms) 
renormalization-scheme-consistent expression, LO(x), and the two-loop fits MRSRi, MRSR 2 , 
NLOi and NL02. For the LO(x) fit the HI data chooses a normalization of 1.00, the ZEUS 
data of 1.015, and the BCDMS data of 0.975. The CCFR data is fixed at a normalization of 0.95, 
and the rest is fixed at 1.00. Similarly, for the NLOi ht the HI data is fixed at a normalization 
of 0.985, the ZEUS chooses a normalization of 0.99, and the BCDMS data of 0.975. Again the 
CCFR data is fixed at a normalization of 0.95, and the rest fixed at 1.00. Also, for the NLO 2 fit 
the HI data is fixed at a normalization of 0.985, the ZEUS chooses a normalization of 0.985, and 
the BCDMS data of 0.97. Again the CCFR data is fixed at a normalization of 0.95, and the rest 
hxed at 1.00. In the Ri and R 2 fits the BCDMS data has a fixed normalization of 0.98, the CCFR 
data of 0.95 and the rest of 1.00. 


Experiment 

data 

points 

LO(x) 

NLOi 

NLO 2 

Ri 

R 2 

HI 

193 

123 

145 

145 

158 

149 

ZEUS .^2^ 

204 

253 

281 

296 

326 

308 

BCDMS 

174 

181 

218 

192 

265 

320 

NMC 

129 

122 

131 

148 

163 

135 

NMC 

129 

114 

107 

125 

134 

99 

NMC 

85 

142 

137 

138 

136 

132 

E665 

53 

63 

63 

63 

62 

63 

CCFR 

66 

59 

48 

40 

41 

56 

CCFR 

66 

48 

39 

36 

51 

47 


Table 2 

Comparison of quality of fits using the full leading-order (including leading-ln(l/x) terms) 
renormalization-scheme-consistent expression, LO(x), and the two-loop fits NLOi and NLO 2 . The 
hts are identical to above, but the data are presented in terms of whether x is less than 0.1 or not. 



data 

points 

LO(x) 

NLOi 

NLO 2 

X > 0.1 

551 

622 

615 

595 

X < 0.1 

548 

483 

554 

589 

total 

1099 

1105 

1169 

1184 
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Figure Captions 


Fig. 1. The description of the EMC and preliminary HI data for using the LO(x) fit. 

Fig. 2. The curves correspond to the value of the proton structure function obtained 

from the leading-order, renormalization-scheme-consistent (LO(x)) calculation at 12 val¬ 
ues of X appropriate for the most recent HERA data. For clarity of display we add 
0.5(12—i) to the value of .^ 2 ( 2 :, Q'^) each time the value of x is decreased, where i = 1 ^ 12. 
The data are assigned to the x value which is closest to the experimental x bin (for more 
details see the similar figure displaying the two-loop fits in [26]). E665 data are also shown 
on the curves with the five largest x values. The HI and ZEUS data are normalized by 
1.00 and 1.015 respectively in order to produce the best fit. 

Fig. 3. The description of the EMC and preliminary HI data for .Ff (x, using the NLOi fit. 

Fig. 4. Comparison of predictions for using the full leading-order, renormalization- 

scheme-consistent (LO(x)) fit and the two-loop NLOi fit. For both sets of curves 
Tl{x,Q'^) increases with increasing at the lowest x values. 
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